[frankinvini]

R это не массив, это просто расстояние вдоль координатно-временной сетки
переделал как то так:

R:=Rs;
for i:=1 to N - 1 do
begin
V[i]:=0.125*Rs*Rs*(q^(2*i))*((q*q+2*q-(2/q)-(1/(q*q)));
A:=0.5*pepr*dT*((q+1)/(q-1));
C[i]:=(V[i]+pepr*dT*((q+1)/(q-1)));
B:=0.5*pepr*dT*((q+1)/(q-1));
F[i]:=V[i]*P2[i];
R:=R+Rs*(q^(i))*(1-(1/q));
end;

здесь избавился от всех непонятных значений типа R(i+-1/2). в этом случае V[i] это просто число, которое меняется в зависимости от i. q^(2*i) это как бы q в степени 2i. C[i] и F[i] также числа. они нужны длф нахождения прогоночных коэффициентов альфа и бетта
p.s. есть полный код программы. в нем как раз таки разобраться не могу. если бы кто нибудь объяснил что там такое, был бы очень благодарен
p.p.s. программа эта моделирует давление, если что=)