5,стр43

uses crt;
var i:integer;
x,y,xs,ys,h:real;
yp,yi,ym:array[0..20]of real;
function f(X,Y:REAL):REAL;
begin
f:=(exp(y)-1)/x;
end;
function fprx(x,y:real):real;
begin
fprx:=-(exp(y)-1)/(x*x);
end;
function fpry(x,y:real):real;
begin
fpry:=exp(y)/x;
end;

begin
clrscr;
xs:=2;
ys:=1;
h:=0.1;
yp[0]:=1;
yi[0]:=1;
ym[0]:=1;
write('x ',xs:2:3,' ');
for i:=1 to 10 do begin
x:=2+i*0.1;
write(x:2:3,' ');
end;
writeln;
writeln('Їа®бв®© ¬Ґв®¤ ќ©«Ґа ');
write('y ',ys:2:3,' ');
for i:=1 to 10 do begin
x:=2+i*0.1;
y:=ys+h*f(xs,ys);
write(y:2:3,' ');
ys:=y;
yp[i]:=y;
xs:=x;
end;
writeln;
writeln('ЁбЇа ў«Ґ­­л© ¬Ґв®¤ ќ©«Ґа ');
write('y ');
ys:=1;
xs:=2;
write(ys:2:3,' ');
for i:=1 to 10 do begin
x:=2+i*0.1;
y:=ys+(2*f(xs,ys)*h+h*h*(fprx(xs,ys)+f(xs,ys)*fpry(xs,ys)))/2;
write(y:2:3,' ');
ys:=y;
xs:=x;
yi[i]:=y;
end;
writeln;
writeln('¬®¤ЁдЁжЁа®ў ­­л© ¬Ґв®¤ ќ©«Ґа ');
write('y ');
ys:=1;
xs:=2;
write(ys:2:3,' ');
for i:=1 to 10 do begin
x:=2+i*0.1;
y:=ys+h*f(xs+h/2,ys+h*f(xs,ys));
write(y:2:3,' ');
ys:=y;
xs:=x;
ym[i]:=y;
end;
writeln;
writeln('в®з­®Ґ аҐиҐ­ЁҐ');
write('y ');
ys:=1;
write(ys:2:3,' ');
for i:=1 to 10 do begin
x:=2+i*0.1;
y:=-ln(1-0.317*x);
write(y:2:3,' ');
end;

writeln;
writeln('¬Ё«­  Їа®бв®©');
write('y ');
for i:=3 to 10 do
yp[i+1]:=yp[i-3]+
(4/3)*h*(2*f(2+(i-2)*0.1,yp[i-2])-f(2+(i-1)*0.1,yp[i-1])+2*f(2+i*0.1,yp[i]));

for i:=1 to 10 do
yp[i+1]:=yp[i-1]+
h/3*(f(2+(i+1)*0.1,yp[i+1])+4*f(2+i*0.1,yp[i])+f(2+(i-1)*0.1,yp[i-1]));
for i:=0 to 10 do write(yp[i]:2:3,' ');

writeln;
writeln('¬Ё«­  ЁбЇа ў«Ґ­­л©');
write('y ');
for i:=3 to 10 do
yi[i+1]:=yi[i-3]+
(4/3)*h*(2*f(2+(i-2)*0.1,yi[i-2])-f(2+(i-1)*0.1,yi[i-1])+2*f(2+i*0.1,yi[i]));

for i:=1 to 10 do
yi[i+1]:=yi[i-1]+
h/3*(f(2+(i+1)*0.1,yi[i+1])+4*f(2+i*0.1,yi[i])+f(2+(i-1)*0.1,yi[i-1]));
for i:=0 to 10 do write(yi[i]:2:3,' ');

writeln;
writeln('¬Ё«­  ¬®¤ЁдЁжЁа®ў ­­л©');
write('y ');
for i:=3 to 10 do
ym[i+1]:=ym[i-3]+
(4/3)*h*(2*f(2+(i-2)*0.1,ym[i-2])-f(2+(i-1)*0.1,ym[i-1])+2*f(2+i*0.1,ym[i]));

for i:=1 to 10 do
ym[i+1]:=ym[i-1]+
h/3*(f(2+(i+1)*0.1,ym[i+1])+4*f(2+i*0.1,ym[i])+f(2+(i-1)*0.1,ym[i-1]));
for i:=0 to 10 do write(ym[i]:2:3,' ');

readln;
end.