Программа зациклилась

[Seriy9292]

Пишу программу интерполяции кубическим сплайном, столкнулся с такой проблемой: программа зацикливается на следующем моменте:
r[d]:=3*((F[d]-F[d-1])/h[d]-(F[d-1]-F[d-2])/h[d-1]);
Я пытаюсь посчитать этот коэффициент, а на нём происходит полный вис программы всю голову сломал, не пойму что не так. Для пояснения:
F-это значение функции cos(x) в точках Xi, которые я задал ранее,h-это шаг между Xi.Походу компилятор не может посчитать значение функции F[d] или F[d-1] или же F[d-2], поэтому зацикливается. Помогите советом, что можно сделать, уже голова пухнет. Вот полный код, зацикливание в конце.

Uses Crt; 
var Xi,Fi,F,ai,bi,ci,di,h,k,l,r,s:array[1..100] of real; 
var x:real; 
var T,n,m,p,o,b,g,j,v,d:byte; 
begin 

Clrscr; 
writeln('vvedite kolichestvo elementov'); 
read(T); 
writeln; 
for n:=0 to T-1 do 
 begin 
write('vvedite X[',n,']): '); 
readln(Xi[n]); 
 end; 
 writeln; 
 For o:=1 to T-1 do 
 begin 
h[o]:=Xi[o]-Xi[o-1]; 
writeln('h[',o,']= ',h[o]:4:2); 
 end; 

writeln; 
writeln('vi vveli sleduushie znacheniya Xi: '); 

For m:=0 to T-1 do 
 begin 
writeln(Xi[m]:3:2); 
writeln; 
 end; 
writeln; 
For b:=0 to T-1 do 

 begin 
 writeln; 

 F:=cos(Xi); 
 writeln('F[',b,']= ',F:4:2); 
 ai[j]:=F[b-1]; 
 write('ai[',b,']= ',ai[j]:4:2); 

 end; 

writeln; 
For v:=1 to T do 
begin 
s[v]:=2*(h[v]+h[v-1]); 
writeln('s[',v,']= ',s[v]:4:2); 
end; 
writeln; 

For d:=1 to T do 
begin 
r[d]:=3*((F[d]-F[d-1])/h[d]-(F[d-1]-F[d-2])/h[d-1]); 
writeln('r[',d,']= ',r[d]:4:2); 
end; 
k[1]:=0;l[1]:=0; 
end.

Помогите пожалуйста




___________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.

[Seriy9292]

Извиняюсь, модератор, я новичек, не знал

[Serge_Bliznykov]

сразу не ходя далеко

var Xi,Fi,F,ai,bi,ci,di,h,k,l,r,s:array[1..100] of real; 
......
for n:=0 to T-1 do 
 begin 
write('vvedite X[',n,']): '); 
readln(Xi[n]); 
 end;

что это Вы в нулевой элемент пишете, если нумерация массива у вас от единицы...

[Seriy9292]

Дружище, спасибо большое, сразу зацикливаться перестало, но появилась ошибка в том же самом вычислении r[d] (Деление на нуль). Где я оплошался?

[Vanta11a]

Сделайте отладку.
Самый простой способ при каждой итерации открываем файл, записываем в него коэффициенты и номер итерации, закрываем.
В результате имеем конечную итерацию и значения коэффициентов перед ошибкой, что позволяет найти косяк.

[Пепел Феникса]

Самый простой способ при каждой итерации открываем файл, записываем в него коэффициенты и номер итерации, закрываем.

О, Господи...
в паскале есть отладчик, просто идем построчно и смотрим значения переменных/выражений.

[Vanta11a]

Да, отладчик хороший. При числе итераций до 100. А если программа валится на тысячной итерации или дальше?