s12rout5.m

来自「Industrial Mathematics」· M 代码 · 共 9 行

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