function [ y ] = Euler1( h )
%Euler1 Approximate the solution of y' = y, x \in [a,b],
%    y_0 = 1 using Euler method.
%   Input: h = mesh size. Output = y_N.

x = 0:h:1;
N = length(x);

y = 1;
for i=1:N-1
    y = y + h*y;
end

end