Задание с курсовой.программа для решения дифференциальных уравнений произвольного порядка.

[saydmc]

люди,помогите с заданием...делаю курсовую,вылезает ошибка в этой строчке:

program kurs;
const jm = 2;
type mas = array[1..jm] of real;
var
y0, y, f: mas;
t0, t, h, tk: real;
j: 1..jm;
i, l: integer;

procedure RP(var t0: real; var y0, f: mas);
begin
f[1] := y[2];
f[2] := y[2] * (y[2] + y[1] * ln * y[1]) * (1 / y[1]);
end;

procedure EU(var t0, h, t: real; var y0, y, f: mas);
var
j: 1..jm;
begin
RP(t0, y0, f);
for j := 1 to jm do
y[j] := y0[j] + h * f[j];
t := t0 + h;
end;

begin
t0 := 0.5; tk := 1.0;
y0[2] := 1.0;
h := 0.001; l := 0;
i := 100;
while t <= tk do
begin
EU(t0, h, t, y0, y, f);
if (l mod i) = 0 then writeln('t =', t0, ' ', 'y1 =', y[1]:5:5, ' ', 'y2= ', y[2]:5:5);
l := l + 1;
for j := 1 to jm do
y0[j] := y[j];
t0 := t;
end;
end.


Program1.pas(18) : Нельзя преобразовать тип function(x: real): real к real

Что именно нужно поменять?

[saydmc]

Ну что)))Кто-нибудь помочь может?

[GetMax]

ln * y[1]

От чего логарифм натуральный берется? Если от y[1] то записать надо так

Ln(y[1])

[saydmc]

Ну в общем я переделал всё это уравнение...вышло так:

program kurs;
const jm = 2;
type mas = array[1..jm] of real;
var
y0,y, f: mas;
t0, t, h, tk: real;
j: 1..jm;
i, l: integer;

procedure RP(var t0: real; var y, f: mas);
begin
f[1] := y[2];
f[2]:=(y0[2]*y0[2])/y0[1]-y0[1]*y0[2]*ln(y0[1]);
end;

procedure EU(var t0, h, t: real; var y0, y, f: mas);
var
j: 1..jm;
begin
RP(t0, y0, f);
for j := 1 to jm do
y[j] := y0[j] + h * f[j];
t := t0 + h;
end;

begin
t0 := 0.5; tk := 1.0;
y[2] := 1.0;
h := 0.001; l := 0;
i := 100;
while t <= tk do
begin
EU(t0, h, t,y0, y, f);
if (l mod i) = 0 then writeln('t =', t0, ' ', 'y1 =', y[1]:5:5, ' ', 'y2= ', y[2]:5:5);
l := l + 1;
for j := 1 to jm do
y0[j] := y[j];
t0 := t;
end;
end.

но она не считает точки!!!
t =0.5 y1 =0.00000 y2= NaN
t =0.6 y1 = NaN y2= NaN
t =0.7 y1 = NaN y2= NaN
t =0.8 y1 = NaN y2= NaN
t =0.9 y1 = NaN y2= NaN

кто-нить может саму программу проглядеть,что в теле не так?

[saydmc]

Никто не глянет?

[saydmc]

Видно никто....(((
Прост сижу с ней весь вечер,сдавать завтра,а в голову уже не лезет ничё,она уже кругом идёт...надеялся на вашу помощь...