История изменений
Исправление cvs-255, (текущая версия) :
Аппроксимация кривой не параболой, а полиномом 4 порядка.
Коэффициенты для 3 точек: -1/(2*H), 0, 1/(2*H)
Коэффициенты для 5 точек: 1/(12*H), -2/(3*H), 0, 2/(3*H), -1/(12*H)
3-х точечной производной хватает для 3-х оборотов вокруг черной дыры по фотонной сфере. 5-точечной для 9 оборотов (тут еще long double понадобился). Использовал рунге-кутта 4 порядка
Исправление cvs-255, :
Аппроксимация кривой не параболой, а полиномом 4 порядка.
Коэффициенты для 3 точек: -1/(2*H), 0, 1/(2*H)
Коэффициенты для 5 точек: 1/(12*H), -2/(3*H), 0, 2/(3*H), -1/(12*H)
3-х точечной производной хватает для 3-х оборотов вокруг черной дыры по фотонной сфере. 5-точечной для 9 оборотов (тут еще long double понадобился)
Исходная версия cvs-255, :
Аппроксимация кривой не параболой, а полиномом 4 порядка.
Коэффициенты для 3 точек: -1/(2*H), 0, 1/(2*H)
Коэффициенты для 5 точек: 1/(12*H), -2/(3*H), 0, 2/(3*H), -1/(12*H)
3-х точечной производной хватает для 3-х оборотов вокруг черной дыры по фотонной сфере. 5-точечной для 9 оборотов