4th-order Runge-Kutta integrator, fully vectorized.


Updated 11 Dec 2015

View License

This is a simple implementation of a 4th-order, fixed-step, Runge-Kutta integrator. It is fully vectorized, which allows a large number of simulations to be carried out in parallel. It is designed to be used inside of optimization code.
% z = rk4(dyn,t,z0)
% This function (rk4) is used to perform a 4th-order Runge-Kutta
% integration of a dynamical system.
% dyn = handle of the form:
% dz = userFunc(t,z)
% t = [1 x nTime] vector of times, created by linspace
% z0 = [nState (x nSim)] matrix of initial states
% z = [nState (x nSim) x nTime ] matrix of trajectories

Cite As

Matthew Kelly (2023). RungeKutta4 (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2012a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes

Added a text function to show usage. Improved description.