Показать сообщение отдельно
  #3  
Старый 17.06.2013, 12:19
icWasya icWasya вне форума
Местный
 
Регистрация: 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;
Ответить с цитированием