(* въвеждаме първоначалната информация *)
a=1.;b=2;n=10;
f[x_,y_]:=(y/x)+1
x0=a;y0=0;
yt[x_]:=x Log[x]
(*Plot[yt[x],{x,a,b}]*)
(* съставяне на мрежата *)
h=((b-a)/n);
xi=Table[a+i*h,{i,0,n}];
For[i=0;x=x0;y=y0,
 i≤n,i++,
Print["i = ",i, " (Subscript[x, i]) = ",x, " (Subscript[y, i]) = ",y, " (Subscript[f, i]) = ",f[x,y], " (SuperStar[(Subscript[y, i])]) = ",yt[x], " (Subscript[e, i]) = ",Abs[y-yt[x]]];
y=y+h*f[x,y];
x=x+h
 ]
Print["Локалната грешка на изчисленията е = ",(h^2)]
Print["Глобалната грешка на изчисленията е = ",h]

допълваме визуализация на решението:

(* въвеждаме първоначалната информация *)
a=1.;b=2;n=5;
f[x_,y_]:=(y/x)+1
x0=a;y0=0;
yt[x_]:=x Log[x]
grt=Plot[yt[x],{x,a,b}];
(* съставяне на мрежата *)
h=((b-a)/n);
xi=Table[a+i*h,{i,0,n}];
points={{x0,y0}};
For[i=0;x=x0;y=y0,
 i≤n-1,i++,
(*Print["i = ",i, " (Subscript[x, i]) = ",x, " (Subscript[y, i]) = ",y, " (Subscript[f, i]) = ",f[x,y], " (SuperStar[(Subscript[y, i])]) = ",yt[x], " (Subscript[e, i]) = ",Abs[y-yt[x]]];*)
y=y+h*f[x,y];
x=x+h;
AppendTo[points,{x,y}]
 ]
Print["Локалната грешка на изчисленията е = ",(h^2)]
Print["Глобалната грешка на изчисленията е = ",h]
grpoints=ListPlot[points,PlotStyle->{PointSize[0.03],Red}];
Show[grt,grpoints]

(* въвеждаме първоначалната информация *)
a=1.;b=2;n=50;
f[x_,y_]:=(y/x)+1
x0=a;y0=0;
yt[x_]:=x Log[x]
grt=Plot[yt[x],{x,a,b}];
(* съставяне на мрежата *)
h=((b-a)/n);
xi=Table[a+i*h,{i,0,n}];
points={{x0,y0}};
For[i=0;x=x0;y=y0,
 i≤n-1,i++,
y=y+h*f[x,y];
x=x+h;
AppendTo[points,{x,y}]
 ]
Print["Локалната грешка на изчисленията е = ",(h^2)]
Print["Глобалната грешка на изчисленията е = ",h]
grpoints=ListPlot[points,PlotStyle->{PointSize[0.03],Red}];
Show[grt,grpoints]