s12rout5.m

Industrial Mathematics

% Routine 12.5  Runge-Kutta  
T = [0,7]; 
z0=[2,2]; 
  [t,z] = ode45('fn',T,z0); 
x = z(:,1);
y = z(:,2);
plot(x,y)