How to Determine Taylor Order for M2 in Runge Kutta Method?

MaxManus
Messages
268
Reaction score
1

Homework Statement


h is the timestep

Y' = F(Y)

Mk,1 = F(Yk)

Mk,2 = F(Yk + 0.5*Mk,1)

Mk,3 = F(Yk + h*Mk,2)

Mk,4 = F(Yk + h*Mk,3)

Yk+1 = Yk + (h/6)*(Mk,1 + 4Mk,2 + Mk,4)

Show that the local error is of fourth order


The Attempt at a Solution



I have written down the fourth order tayler series of Y

and the equation for Runge Kutta k = 0

What I need help with is what order I am supposed to write the tayler series for M4, M3 and M2
 
Physics news on Phys.org
Fourth order means that the error is on the order of h5. What does that tell you about the order of those Taylor series?
 
Thanks, my taylor series for Y has an error of order 4 so then it is of third order and not fourth.

I'm not sure about your question. Are you saying that I'm supposed to write the third order taylor series for M2,M3 and M4?
 
Last edited:
My Taylor series for Y

Y(t1) = Y(to) + Y'(t0)h + 0.5Y''(t0)h^2 + (1/6)*Y'''(t0)h^3 + (1/24)*Y''''(E)h^4

E is between t0 and t1
 
Your Taylor series should be fifth order.
 
If Y(t1) is my tayler series and Y1 is my RungeKutta expression
I thought I should get:
Y1 = Y/t0) + Y'(to)h + 0.5Y''(t0)h^2 + (1/6)Y'''h^3 + Q*h^4

Where Q is a function of Y and its derivatives

To find the local error
abs(Y(t1) - Yt1) = W*h^4

where W is a function of Y and its derivatives.

But this is wrong?
 
An nth order method means that the error is on the order of hn+1. This means your Taylor series should be fourth order plus a fifth order error term to show that the method is fourth order.
 
But I'm supposed to show that the local error is of fourth order. That means that the method is of third order?
And I should use third order Taylor series for Y(t1)?
 
Solved with two questions
My solution:
M4: Second order Taylor
M3: First order Taylor
M2: Here I used second order for the term that stands alone and zero order for the one inside M3

The trick was to do M4 first to see what Taylor order you needed for M3 to keep the error of order four

Am I allowed to use to different Taylor orders for M2. And is zero order allowed?
 
Back
Top