решение на системата от задача 3  стр.34

(* въвеждаме първоначалната информация *)
a=1.;b=1.5;
(*n=5;h=((b-a)/n);(* в случай, че е дадено n и търсим h *)*)
h=0.1;n = ((b-a)/h);
f[x_,y_,z_]:=x -2z
g[x_,y_,z_]:=z+((3y)/(x+z))
x0=a;y0=-1;z0=2;

For[i=0;x=x0;y=y0;z = z0,
 i≤n,i++,
Print["i = ",i,
 " (Subscript[x, i]) = ",x, 
" (Subscript[y, i]) = ",y,
" (Subscript[z, i]) = ",z, 
" (Subscript[f, i]) = ",f[x,y,z], 
" (Subscript[g, i]) = ",g[x,y,z]];
y1=y+h*f[x,y,z];
z1=z+h*g[x,y,z];
x=x+h; y=y1;z=z1
 ]
Print["Локалната грешка на изчисленията е = ",(h^2)]
Print["Глобалната грешка на изчисленията е = ",h]