
08.06.2010, 17:27
|
Прохожий
|
|
Регистрация: 08.06.2010
Сообщения: 5
Репутация: 10
|
|
Что не так?
Это метод Эйлера с заданной точностью. написанно в Turbo Delphi.
Код:
program Eiler;
{$APPTYPE CONSOLE}
uses
SysUtils,
InOut,
Math;
const
E = 0.01;
var
x,y,g,m,y0,x0 :Extended;
n,i: integer;
function F(x,y:extended):Extended;
begin
Result := 0.2*(power(y,0.5))*x;
end;
begin
Write(RusWD('Введите значение Хо: '));
readln(x);
Writeln;
Write(RusWD('Введите значение Yо: '));
readln(y);
Writeln;
Write(RusWD('Введите значение Хк: '));
readln(m);
Writeln;
Write(RusWD('Введите шаг: '));
readln(g);
Writeln;
Writeln(' X ',' ',' Y ');
y:= y0;
while abs(y0-y) <= E do
begin
y:= y0;
n:= (x+m)/g;
for i := 1 to n do
begin
y := y+g*f(x,y);
x := x+g;
end;
end;
Writeln(x:7:5,' ',y:7:5);
Readln;
end.
|