|Oct7-11, 04:03 AM||#1|
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.
|Oct7-11, 05:58 AM||#2|
Matlab already has a built-in Runge-Kutta solver, probabaly ode45. So I think you need not have to worry about those ki's.
|Similar Threads for: Runga-Kutta help needed (Matlab)|
|Modelling the spread of disease (in MatLab, Runge-Kutta)||Engineering, Comp Sci, & Technology Homework||0|
|MATLAB Help needed!||Engineering, Comp Sci, & Technology Homework||1|
|Eulers & Runge-Kutta using Matlab||Calculus & Beyond Homework||0|
|Matlab Problem: Runge-Kutta Fehlberg||Calculus & Beyond Homework||0|