clear;

h    = 0.1;
N    = 10;
x    = zeros(N+1, 1);
y    = zeros(N+1, 1);

x(1) = 0;
y(1) = 1;

for (i = 1:N)
  x(i+1) = x(i) + h;
  y(i+1) = y(i) + h * x(i) * y(i);
end

[x y]