Скачиваний:
26
Добавлен:
09.12.2013
Размер:
909 б
Скачать
program Miln;
uses crt;
var
i,j:integer;
x,h,f:real;
x1,y,y1: array [0..20] of real;
begin
clrscr;
y[0]:=0.2;h:=0.1;X:=0;

for i:=0 to 3 do begin {Њ®¤ЁдЁжЁа®ў ­­л© ¬Ґв®¤ ќ©«Ґа }
y[i+1]:=y[i]+h*(0.3*(x+h/2)+sqr(y[i]+(h/2)*(0.3*(x)+sqr(y[i]))));
x:=x+h;
end;x:=0;

for i:=0 to 10 do begin {‚лзЁб«пҐ¬ §­ зҐ­Ёп 㧫®ў}
x1[i]:=x;
x:=x+h;
end;

for i:=4 to 10 do begin
y[i]:=y[i-4]+(4*h/3)*(2*(0.3*(x1[i-3])+ {”®а¬г«  Їа®Ј­®§ }
sqr(y[i-3]))-(0.3*(x1[i-2])+sqr(y[i-2]))+
2*(0.3*(x1[i-1])+sqr(y[i-1])));

f:= (sqr(x1[i])+sqr(y[i]));

y1[i]:=y[i-2]+(h/3)*(f+4*(0.3*(x1[i-1])+ {”®а¬г«  Є®а४樨}
sqr(y[i-1]))-(0.3*(x1[i-2])+sqr(y[i-2])));
end;
writeln(' Їа®Ј­®§ Є®а४жЁп');
for i:=0 to 10 do
if i<>10 then writeln(i,' ',x1[i]:2:1,' ',y[i]:8:7,' ',y1[i]:8:7) else
writeln(i,' ',x1[i]:2:1,' ',y[i]:8:7,' ',y1[i]:8:7);
readln;
end.
Соседние файлы в папке Выч методы