24 Reputation

15 years, 79 days

Promlem soved...

Thank you gentlemen - problem solved.

plot3d - orientation...

Thanks Alex, I tried this but do not see the values. I should have mentioned that I am using Maple version 8. Regards

How to simplify taylor series...

I am trying to write a script that generates the truncation error of a discretization. For example, a 4-point (equaly spaced), finite difference formula for the second derivative of y(x) is ` y''(x) = (2*y1-5*y2+4*y3-y4)/(h^2) (1)` My approach has been to expand y(x) as a Taylor series and substitute into (1) evaluated at y1=y(x=0), y2=y(x=h) and y3=y(x=2h). However, whilst I obtain the correct Taylor series expansions, I cannot simplify to get the truncation error. The code I am attempting to use is as folows: ``` > restart;dy :=taylor(y(x),x=0,5): > y1:=subs(x=0,dy): > y2:=subs(x=h, dy): > y3:=subs(x=2*h, dy): > y4:=subs(x=3*h, dy): > dy1:=simplify((2*y1-5*y2+4*y3-y4)/(h^2)); ``` This produces the correct output but it is not simplified. I would like to be able to simplify it to get the truncation error. Any help appreciated. Graham
 Page 1 of 1
﻿