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: 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
Flapping baby birds give clues to origin of flight
Prions can trigger 'stuck' wine fermentations, researchers find
Socially-assistive robots help kids with autism learn by providing personalized prompts
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