Euler.m

function Y = Euler

b = 6;

n = 500;

m = n + 1; % m = n+1

h = b/n;

x = 0:h:b; % x(1), x(2), ..., x(n+1)

Y = zeros(1,m);

Y(1) = 1;

for k = 1: n

Y(k+1) = Y(k) + h*(1 - x(k)*Y(k));

end

plot3(x,y,z)