Register to reply

Runga-Kutta help needed (Matlab)

by Mech-Master
Tags: matlab, rungakutta
Share this thread:
Mech-Master
#1
Oct7-11, 04:03 AM
P: 13
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
Physical constant is constant even in strong gravitational fields
Montreal VR headset team turns to crowdfunding for Totem
Researchers study vital 'on/off switches' that control when bacteria turn deadly
matematikawan
#2
Oct7-11, 05:58 AM
P: 333
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