s12rout1.m

Industrial Mathematics

% Routine 12.1   Euler's Method
clear
x=0;
y = 1;
h=0.1;
N = 11;
                   %
for i = 1:N
 a(i) = x;
 b(i) = y;
 y = y + h*fn(x,y);
 x = x + h;
 end;
plot(a,b)