обратная трассировка лучей и тор _ ошибка в коде или рассуждениях?

[Anna221]

Помогите, пожалуйста, никак не могу найти ошибку. Вот кусок кода, из-за которого, как я предполагаю, вылетает ошибка:

tmp5 := (a4*(sqr(t)*sqr(t))+ a3*(t*sqr(t))+ a2*(sqr(t)) + a1*t + a0) ;
tmp4 := trunc(tmp5);
if ( tmp4 < 2 ) then result := t
else result := -1;

если ставить if tmp4=0 then то программа работает, но не находит лучей, пересекающих тор, т.е., как я понимаю, нужно учитывать погрешность вычислений, но тогда вылетает ошибка..

tmp4 - integer, tmp5 - double

Может, все элементарно, но никак не могу придумать, как исправить..

и еще, ход мыслей при нахождении нормали к тору и пересечение луча с тором описаны верно? (с комментариями в torusUnit) или нужен другой подход?
Заранее благодарю
Код прилагается.