MATLAB Runga-Kutta help needed (Matlab)

An object of mass m falls from rest at a point near the earth's surface. If the air resistance is proportional to the velocity v^2, the differential equation for the velocity as a function of time is given by:

m*dv/dt = mg - cv^2

a) Derive the exact solution

done this part, and i got v = sqrt(m*g/c)*tanh(t*sqrt(g*c/m))

b) For the given paraments g = 9.81 m/s^2. m = 68.1 kg and c = 1.5 kg/m. plot the exact solution and the numerical solution v(t) obtained from the 4th order predictor-corrector runge kutta methods using an interval of dt = 0.25 seconds in the domain of 0<t<6

i'm having trouble coding the runge kutta method with all the k1, k2, k3 and k4.. I really need help with this.

Thanks
 
Matlab already has a built-in Runge-Kutta solver, probabaly ode45. So I think you need not have to worry about those ki's.
 

Related Threads for: Runga-Kutta help needed (Matlab)

Replies
2
Views
1K
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
6
Views
4K
Replies
1
Views
3K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
0
Views
2K
Replies
7
Views
10K
Top