File Exchange

image thumbnail

Runge Kutta 4th order ode

version 1.4 (1.32 KB) by

solves ode using 4th order Runge Kutta method

4.66667
7 Ratings

102 Downloads

Updated

View License

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.

Comments and Ratings (13)

how can i solve SIR model using RK4 method in matlab? can you write the code please

faiz islam

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 .
regards
faiz

Ali Abbas

@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?
thank you

Shahzaib Asif

function RK4(f,a,x0,y0,h)

% 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;

for i=1:(length(x)-1)
   
    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;

end

y(:)

%Shahzaib Asif (zaibi7402)
%shahzaib.7402@gmail.com

Chris FUNG

clear coding

Christoph

Ying

Ying (view profile)

Ying

Ying (view profile)

Very good to learn. Thanks.

Arun

Arun (view profile)

Pi Ting

excellent work

Ido

Ido (view profile)

Excellent program,
very helpful.

Updates

1.4

Just an update

1.1

Nothing much

MATLAB Release
MATLAB 7.10 (R2010a)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video