Интерполяция Лагранжа

[Smitt&Wesson]

Второй день бьюсь, никак не могу понять, что я делаю не так.
Есть формулы и программа на Бейсике. Есть в описании контрольный пример. Но, беда втом, что результаты выдаваемые программами на Бейсике и С++, не совпадают с контрольным примером.
Вот условие: Я знаю, ошибка в представлении данных, но не погу понять, где именно.

[WinCoder]

Так совпадают же. Просто во втором случае идёт округление до 6 символов после запятой, а в первом до 7.

[Smitt&Wesson]

Сообщение от WinCoder

Так совпадают же. Просто во втором случае идёт округление до 6 символов после запятой, а в первом до 7.

Совпадают на консоли у разных программ. Не совпадение с контрольными примерами (подчёркнуто красным).

[Аватар]

Так может везде двойную точность для констант поставить? Все таки 7 значащих

[Smitt&Wesson]

Да при чём тут точность? Тьфу на неё. 0,6087558 много больше чем 0,4204928.
В кубической интерполяции, округление результата, на сам расчёт не влияет.
И ещё не понятно. В программах выводится один результат - Y, а в примерах даны результаты двух расчётов Y(0,9) и Y(1,1) с какого марса они упали? Формулы на этот счёт вооще молчат.
Догадался. Y(0,9) это X = 0.9 но результат всё равно не совпадает с примером.

[Аватар]

Называли бы переменные нормально, то и проблем бы не было
Например 1-ое присвоение игрека заменить на:
y:=(1-p)*c+p*d;

[Smitt&Wesson]

С переменными у меня всё ОК :

// Интерполяция при H=const по формулам Лагранжа
// Степень полинома от 1 до 5
bool InterpolationLagranj(double *Y, double H, double X,
                          double X0, uint degree, double &res)
{
  bool ret = true;
  double p = (X - X0) / H, i, j, k, l, m, y;
  switch(degree)
  {
    case 1:{res = (1 - p) * Y[0] + p * Y[1]; break;}
    case 2:{
             res = p * (p - 1) * Y[0] / 2 + (1 - p * p) * Y[1] + p * (p + 1) * Y[2] / 2;
             break;
           }
    case 3:{
             m = (p - 2) / 2;
             y = -p * (p - 1) * m * Y[0] / 3 + (p * p - 1) * m * Y[1];
             res = y - p * (p + 1) * m * Y[2] + p * (p * p - 1) * Y[3] / 6;
             break;
           }
    case 4:{
             m = (p * p - 1) / 2; k = (p * p - 4) / 2;
             y = m * p * (p - 2) * Y[0] / 12 - (p - 1) * p * k * Y[1] / 3 + m * k * Y[2];
             res = y - (p + 1) * p * k * Y[3] / 3 + m * p * (p + 2) * Y[4] / 12;
             break;
           }
    case 5:{
             m = p * p - 1;
             i = m / 24;
             j = -p * i * (p - 3);
             k = m - 3;
             l = k * (p - 3) / 12;
             y = j * (p - 2) * Y[0] / 5 + l * p * (p - 1) * Y[1] / 2 - m * l * Y[2];
             res = y + p * (p + 1) * l * Y[3] + j * (p + 2) * Y[4] + p * i * k * Y[5] / 5;
             break;
           }
    default : ret = false;
  }
  return ret;
}

Формулы тоже проверены - перепроверены 100 раз. Где-то у меня с теорией пробел, а именно с правильностью задания коэффициентов.
Исходная прога, на бейсике:

[Аватар]

Нет, не порядок. Почему у меня y:=(1-p)*c+p*d; дает 0,6087558. Дальше формулы не смотрел
смотрим что такое c и d, и что такое a и b как в вашей формуле. И что такое y0 и y1. В исходной проге чепуха

[Smitt&Wesson]

Сообщение от Аватар

Почему у меня y:=(1-p)*c+p*d; дает 0,6087558.

Может я Р не правильно вычисляю т.к. ошибка во всех формулах? И почему у вас c+p*d если даже в формуле Y0 и Y1, а вводится A , B?
P = (X - X0) / H = (1 - 1.1) / 0.2 = 0.5
Y = (1 - P) * Y0 + P * Y1 = (1 - 0.5) * 0.36699114 + 0.5 * 0.4739944 = 0,42049277. Посчитал вручную на калькуляторе

[Аватар]

В примере y(0)=0.5686566
С какого перепугу вместо него a подставлять? а=y(-2); b=y(-1); c=y(0) и т.д.
Я про код в #1