Знайти розв’язок диференціального рівняння y’=x2 y2 -1 на інтервалі а=0 b=1 методом Ейлера-Коші

При початковій умові y0(0)=1

Рішення методом Ейлера-Коші в середовищі Turbo Pascal

Program Е_К;

uses crt ;

var a,b,h,y0,y,x,k1,k2,k3,k4:real;

i,n:integer;

label @;

Begin

clrscr;

a:=0;

b:=1;

n:=10;

y0:=1;

x0:=0;

x:=0;

i:=1;

h:=(b-a)/n;

@:

x:=x+h;

y:= y0+(h/2)*(((sqr(x0)*sqr(y0))-1)+(sqr(x)*sqr(y0+h*(sqr(x)*sqr(y))-1)-1));

writeln(‘x=’,x:5:5, ‘y=’,y:5:5);

if i<n then

Begin

i:=i+1;

y0:=y;

goto @;

end;

end.

Реалізація в Mathcad