Лабораторки по информатике / SAM / sam / AQWAR

const n=100;
type dim=array[0..n] of real;
f=function(x,y:real):real;
var od1,od2: dim;
c,x,x0,xn:real;
a:integer;
t:text;
function yr(x,c:real):real;
begin
yr:=x/(sqrt(2)-1);

end;

procedure ejler (n:integer;h,y0,x0,xn:real;var y:dim);
var a:integer;
x:dim;
function yg(x,y:real):real;
begin
yg:=(y-x)/(x+y);
end;
begin
x[0]:=x0;
y[0]:=y0;
for a:=1 to n do begin
x[a]:=x0+h*a;
y[a]:=y[a-1]+h*yg(x[a-1],y[a-1]);
end;
end;

procedure output;
begin
assign(t,'e:\2tm-25d\data.dat');
rewrite(t);
writeln (t,'*************************************************************');
writeln (t,'* x * y(x) * в®з­®бвм 0.1 * в®з­®бвм0.01 *');
writeln (t,'*************************************************************');
x:=x0;
for a:=0 to 9 do begin
writeln(t,'*' ,x:3:1, '*' ,yr(x,c):8:4, '*' ,od1[a]:8:4, '*' ,od2[a*10]:8:4, '*');
x:=x+0.1;
end;
writeln(t,'***************************************************************');
writeln('‚лў®¤ ў д ©« data.dat-OK');
close(t);
end;

begin
writeln('‚ўҐ¤ЁвҐ ­ з «м­лҐ Ё Є®­Ґз­лҐ §­ зҐ­Ёп x');
readln(x0);
readln(xn);
write('‚ўҐ¤ЁвҐ C');
readln(c);
ejler(10,0.1,x0,yr(x0,c),xn,od1);
ejler(100,0.01,x0,yr(x0,c),xn,od2);
output;
end.