This code defines an existing function and step size which you can change as per requirement.
P.S: This code has no new feature compared to existing codes available online. Intention behind posting this very simple code is to help students understand the concept and solve assignments.
It does not work for y(0) = 0 initial condition
Great work! What about a code for Runge Kutta method for second order ODE. Something of this nature:
d^2y/dx^2 + 0.6*dy/dx 0.8y = 0
how can i solve SIR model using RK4 method in matlab? can you write the code please
sir can you assist me ,that how we can apply 4th order Runge kutta method for 4 coupled equation?
dx/dt=−ax − eω + yz
dy/dt= by + xz
dz/dt= cz + fω − xy
dω/dt = dω – gz
a = 50, b =−16, c = 10, d = 0.2, e = 10, f = 16, g = 0.5
Step size 0.001 .
@Shahzaib Asif Very helpful program.JazakAllah
How do I run/call to this code?
for this function : f'''' - f*f''' + 4*g = 0
where i need to insert it in this code?
% Runge Kutta Method 4th Order
% function @(x,y) e.g. f=@(x,y)(x+y);
% a = the point up to which you obtain the results
% x0 = initial condition of x
% y0 = initial condition of y
% step size
x = x0:h:a;
y(1) = y0;
k1 = f(x(i),y(i));
k2 = f(x(i)+0.5*h,y(i)+0.5*h*k1);
k3 = f((x(i)+0.5*h),(y(i)+0.5*h*k2));
k4 = f((x(i)+h),(y(i)+k3*h));
y(i+1) = y(i) + (1/6)*(k1+2*k2+2*k3+k4)*h;
%Shahzaib Asif (zaibi7402)
Very good to learn. Thanks.
Just an update