function [T,Y]=myeuler_2(t,x,t1,n)

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

h=(t1-t)/n;
T=t;
X=[];
X=x


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

fid = fopen('table.txt','w');
fprintf(fid,'              h=%8.4f\n \n',h);
fprintf(fid, '    t         x_1        x_2\n \n'); 
fprintf(fid,'%8.4f  %8.4f %8.4f\n',t,x(1),x(2));

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

for i=1:n
  x=x+h*f2(t,x);
  t=t+h;
  T=[T;t];
  X=[X;x];
  fprintf(fid,'%8.4f %8.4f  %8.4f \n',t,x(1),x(2));
end
fclose(fid);
%...................................
% Plot graph 

figure(1)
clf
plot(T,X(:,1))
xlabel('t')
ylabel('x_1(t)')
Ha_ax=gca;set(Ha_ax,'Fontsize',12);
figure(2)
clf
plot(T,X(:,2))
xlabel('t')
ylabel('x_2(t)')
Ha_ax=gca;set(Ha_ax,'Fontsize',12);

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