s12rout2.m

Industrial Mathematics

% Routine 12.2   Euler's Method on Lotka-Volterra
clear all
N = 1500;
h = 0.005;
t = 0;
z = [2,2]'; 
for i=1:N                  %
a(i) = z(1);
b(i) = z(2);
 z = z + h*fn(t,z);
 t = t + h;
  end;
plot(a,b)