4th Order Runge-Kutta method and over/under estimates

  • Thread starter Thread starter dwdoyle8854
  • Start date Start date
  • Tags Tags
    Method Runge-kutta
Click For Summary
SUMMARY

The forum discussion centers on the application of the 4th Order Runge-Kutta (RK4) method to approximate the solution of the differential equation dF/dt = -0.1F + 70, with initial condition F(0) = 0. The exact solution is identified as F(t) = 700 - 700 * exp(-0.1t). Despite the expectation that RK4 would yield an overestimate due to its linear characteristics, the user found that the method produced an under-approximation. The conversation highlights the importance of understanding the polynomial nature of RK4 and its implications for error analysis.

PREREQUISITES
  • Understanding of differential equations and initial value problems
  • Familiarity with the 4th Order Runge-Kutta method
  • Basic knowledge of polynomial functions and their properties
  • Proficiency in using Excel for numerical computations
NEXT STEPS
  • Study the derivation and application of the 4th Order Runge-Kutta method
  • Learn about error analysis in numerical methods, particularly for RK4
  • Explore the concept of concavity in relation to solution curves of differential equations
  • Investigate the impact of varying step sizes on the accuracy of numerical approximations
USEFUL FOR

Students and educators in mathematics or engineering fields, particularly those focusing on numerical methods for solving differential equations.

dwdoyle8854
Messages
16
Reaction score
0

Homework Statement



"Use Excel to approximate dF/dt=-0.1F+70, F(0)=0 to generate approximations for F at t=1,2 and 4 using step size 0.1. Explain whether these approximation are greater than or less than the exact values. Determine whether the shape of the solution curve is increasing, decreasing, concave up or concave down based on the data alone. Explain."

Homework Equations



ynext = ynow + (1/6)(k1 +2k2 + 2k3 + k4)
k1= Δx*f'(xnow, ynow)
k2= Δx*f'(xnow+.5Δx,ynow +.5k1)
k3= Δx*f'(xnow+.5Δx,ynow +.5k2)
k4= Δx*f'(xnow+Δx,ynow+k3)

I found the exact solution to be F(t)=700-700*exp(-.1t)

The Attempt at a Solution



i've attached my excel file.

Since Runge-Kutta is al inear technique and I observed all the slopes in RK to be decreasing I predicted that the method would give an over estimate. However, comparing my estimate with the exact solution shows that I infact get an under-approximation. I am completely lost as to why this occurs and am looking for some explanation. I do not feel i am able to say anything about the concavity given my results that the decreasing rates (f'1, f'2 et cetera on my excel) gave an underestimation.

futher, how exactly does putting weight on the k2 and k3 terms affect the approximation?




[
 

Attachments

Physics news on Phys.org
RK4 is not a linear technique. Where did you get the idea that it is?
 
doesnt it assume that rate of change is constant over an interval of time? that to me says linear. but nonetheless still don't know what to say about the concavity or the reason why I got an underestimation.
 
dwdoyle8854 said:
doesnt it assume that rate of change is constant over an interval of time? that to me says linear.
No. That's the whole point of evaluating the derivative four times during one step. The "4" in RK4 is short for fourth order. RK4 essentially comes up with a fourth order polynomial for each step.
 
okay, so since its a polynomial, is there anything that can be said about whether we expect it to over/under approximate?
 
It's not some random 4th order polynomial. It's a very specific 4th order polynomial. Which is it? What does that tell you about the error? This is homework, so at this stage I'm leaving the rest up to you.
 
to be honest, i have no idea what you are hinting at with the whole 4th order polynomial thing. I don't know what that has to do with error, or what precisely you mean by a 4th order polynomial.

this is more of a discovery project than a homework. Its graded sure, but its purpose is to explore.
 

Similar threads

  • · Replies 14 ·
Replies
14
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 6 ·
Replies
6
Views
51K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 3 ·
Replies
3
Views
7K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
9
Views
7K