function [X,Y]=myeuler_1(x,y,x1,n)

%................................... 
%set initial condition and compute h

h=(x1-x)/n;
X=x;
Y=y;


%................................... 
%  Print table:

fid = fopen('table.txt','w');
fprintf(fid,'              h=%8.4f\n \n',h);
fprintf(fid, '    x         y     \n \n'); 
fprintf(fid,'%8.4f  %8.4f \n',x,y);

%................................... 
% Euler's method 

for i=1:n
  y=y+h*f(x,y);
  x=x+h;
  X=[X;x];
  Y=[Y;y];
  fprintf(fid,'%8.4f  %8.4f \n',x,y);
end

fclose(fid);

%...................................
% Plot graph 

clf
plot(X,Y)
xlabel('x')
ylabel('y(x)')
Ha_ax=gca;set(Ha_ax,'Fontsize',12);
hold off

%...................................