A1 Matlab

Primer funkcije: Metoda Runge-Kutta 4. reda
function [t, x] = rkutta4( F, tp, tf, x0, T )
% Metod integracije RUNGE-KUTTA 4. reda
% parametri su isti kao kod ode23,
% sem T koje predstavlja korak integracije
if nargin < 5, T = 0.1; end
t = tp : T : tf;
n = length(t);
x = zeros( n, length(x0) );
x(1,:) = x0(:)';
tl = tp;
for i = 1 : n-1
tl = tl + T;
xp = x(i,:)';
k1 = T * feval( F, tl, xp );
k2 = T * feval( F, tl+0.5*T, xp+0.5*k1 );
k3 = T * feval( F, tl+0.5*T ,xp+0.5*k2 );
k4 = T * feval( F, tl+T, xp+k3 );
x(i+1,:) = x(i,:) + (k1/6 + k2/3 + k3/3 + k4/6)';
end