Runge Kutta 4th order ode
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.
Cite As
Judah S (2021). Runge Kutta 4th order ode (https://www.mathworks.com/matlabcentral/fileexchange/29851-runge-kutta-4th-order-ode), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (21)
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Acknowledgements
Inspired: ODE4 gives more accurate results than ODE45, ODE23, ODE23s
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
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
Thank you
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 .
regards
faiz
@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
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
clear coding
Very good to learn. Thanks.
excellent work
Excellent program,
very helpful.