1. If air resistance is proportional to the square of the instantaneous velocity, then the velocity v of a mass m dropped from a given height is determined from:
m*dv/dt = mg - kv^, k>0.

Let: v(0) = 0, k = 0.125, m = 5 slugs, and g = 32 ft/s^2

a) Use a fourth-order Runge-kutta (RK4) method with h = 1 to approximate the velocity v(5)
b) Use excel to graph the solution of the IVP on the interval [0,6].
c) Use separation of variables to solve the IVP and then find the actual value.

2. Find the analytic solution of the initial-value problem (IVP)

y' = -y +10sin(3x), y(0) = 0, on the interval [0,2].

a) Graph the solution and find its positive roots.
b) Use Runge-Kutta (RK4) method with h = 0,1 to approximate a solution of the initial-value problem (IVP)

Hello and thank you for posting your question to Brainmass.
The solution is attached below in two files. the files are identical in ...

Solution Summary

The solution shows in a step by step manner how to solve a first order differential equation using order 4 Runge-Kutta method and comparing it to the analytic solution

$2.19

See Also This Related BrainMass Solution

Runge-Kutta Method Errors

Question:
How does this compare with the theoretical rate of convergence of 0(H2). Explain your result as best as you can?

Data:
We used the Runge-Kutta method to solve
y'(x) = - y(x) + x^0.1[1.1 + x], y(0) = 0
whose solution is y(x) = x^1.1. We solved the equation on [0,5] and we printed the errors at x = 1,2,3,4,5. We used stepsize h-0.1, 0.05, 0.0024, 0.0125, 0.00625. We calculated the ratios by which the errors decrease when h is halved.

The ratios by which the errors decrease when h is halved (see attached).

... Data: We used the Runge-Kutta method to solve y'(x) = - y(x) + x^0.1[1.1 + x], y(0) = 0 ... the equation on [0,5] and we printed the errors at x = 1,2,3,4,5. We ...

Question Use Runge-Kutta method of order four to approximate the solution to the given initial value problem and compare the results to the actual values. ...

... 1. Detailed explanation of methods used (Euler and runge kutta) 2. Matlab code ... strengths of both methods) 3. similarities/differences of methods 4. Plots with ...

... In a Matlab script, we demonstrate the application of the Runge-Kutta numerical method for a Lorenz attractor, the Butterfly Effect caused by a small ... PROBLEM 4. ...

I am posting formula called Runge Kutta Fourth Order. ... rk1 = increment_value*f(t, w); %% To keep all t_i, w_i, make all the t and w in these 4 %% steps as t(i ...

... us write out what the order conditions are for a general r-stage Runge-Kutta method ... 1 3 2 To get third order accuracy, conditions (4) and (5) must be satisfied ...

... 3D plots of Runge-Kutta starting from the two versions of initial conditions. Fig 3 compares plots versus time of the two versions of RK. Fig 4 plots versus ...

... at equilibrium and from equations (4), (1), and (3), we have. ... x(t) and y(t) for initial condition x(0)=0 and y(0) = 0, using the 4th order Runge-Kutta routine. ...

... π π π I= [f (0) + 2f ( ) + f ( )] (6) 8 4 2 This will give 0.7854. We can see that there is not much diﬀerence in the answer. 1 2 Runge-Kutta Method The ...