laba3

program laba3;

const h=0.1;


var x,y,w,t,yp,yk: array[0..40] of real;
i:integer;
q:text;
begin
x[0]:=0; y[0]:=0; w[0]:=0;
assign(q,'rez.txt');
rewrite(q);

{prostoy eyler}

writeln(q,'prostoy');

for i:=0 to 39 do
begin
x[i+1]:=x[i]+h;
w[i+1]:=w[i]+h*(exp(x[i])+5*w[i]-6*y[i]);
y[i+1]:=y[i]+h*w[i];
t[i+1]:=-exp(2*x[i])+0.5*exp(3*x[i])+0.5*exp(x[i]);
writeln(q,w[i+1],y[i+1],t[i+1]);
end;

writeln(q);


{ispr}
writeln(q,'ispr');
for i:=0 to 39 do
begin
w[i+1]:=w[i]+h/2*(5*w[i]-6*y[i]+exp(x[i])+5*(w[i]+h*(5*w[i]-6*y[i]+exp(x[i])))-6*(y[i]+h*w[i])+exp(x[i]+h));
y[i+1]:=y[i]+h/2*(w[i]+w[i]+h*(5*w[i]-6*y[i]+exp(x[i])));
writeln(q,w[i+1],y[i+1],t[i+1]);
end;

writeln(q);


{mod}
writeln(q,'mod');
for i:=0 to 39 do
begin
w[i+1]:=w[i]+h*(exp(x[i]+h/2)+5*(w[i]+h/2*(exp(x[i])+5*w[i]-6*y[i]))-6*(y[i]+h/2*w[i]));
y[i+1]:=y[i]+h*(w[i]+h/2*(5*w[i]-6*y[i]+exp(x[i])));
writeln(q,w[i+1],y[i+1],t[i+1]);
end;


writeln(q);





{milna}

writeln(q,'milna');
for i:=0 to 3 do
yk[i]:=y[i];



for i:=3 to 39 do

begin

yp[i+1]:=yk[i-3]+4*h/3*(2*yk[i-2]-yk[i-1]+2*yk[i]);
yk[i+1]:=yk[i-1]+h/3*(yp[i+1]+4*yk[i]+yk[i-1]);

writeln(q,yp[i+1],yk[i+1],t[i+1]);
end;








close(q);



end.