
17.06.2013, 12:19
|
Местный
|
|
Регистрация: 09.11.2010
Сообщения: 499
Репутация: 10
|
|
надо как-то так
Код:
var
R:Real;
R_i:Real;
R_i_p1:Real;
R_i_m1:Real;
R_i_p05:Real;
R_i_m05:Real;
begin
R := Rs*(1+q);
R_i := Rs;
for i:=1 to N-1 do begin
R_i_m1 := R_i;
R_i := R_i * q;
R_i_p1 := R_i * q;
R_i_m05 := 0.5*(R_i + R_i_m1);
R_i_p05 := 0.5*(R_i_p1 + R_i);
// далее заменяем
V:=0.5*((R_i_p05)*(R_i_p05)-(R_i_m-05)*(R_i_m05));
A[i]:=pepr*dT*((R_i_m05)/(R_i-R_i_m1));
C[i]:=(V+pepr*dT*(((R_i_m05)/(R_i-R_i_m1))-((R_i_p05)/(R_i_p1-R_i));
B[i]:=Pepr*dT*((R_i_p05)/(R_i_p1-R_i);
F[i]:=V*P2[i];
R:=R+(R_i-R_i_m1);
end;
|