Runga-Kutta help needed (Matlab)


by Mech-Master
Tags: matlab, rungakutta
Mech-Master
Mech-Master is offline
#1
Oct7-11, 04:03 AM
P: 10
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
Phys.Org News Partner Science news on Phys.org
NASA's space station Robonaut finally getting legs
Free the seed: OSSI nurtures growing plants without patent barriers
Going nuts? Turkey looks to pistachios to heat new eco-city
matematikawan
matematikawan is offline
#2
Oct7-11, 05:58 AM
P: 330
Matlab already has a built-in Runge-Kutta solver, probabaly ode45. So I think you need not have to worry about those ki's.


Register to reply

Related Discussions
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