Интерполяционный полином ньютона-где то ошибка в коде

[DiGris]

ввожу isk_x=2 x0=1 y0=1 x1=3 x2=2 выводит 1.5, а по формуле 0.5


#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <iostream>
using namespace std;


void main(void){
int n;
double isk_X;
float y;
float x;
cout<<"enter kol tochek"<<endl;
//ввод с клавиатуры числа точек
cin>>n;
//объявление массивов для хранения значений координат точек
float* Xi=new float[n];
float* Yi=new float[n];
//заполнение массивов с клавиатуры

for (int i=0;i<n;i++){
cout<<"enter koord toch po X"<<i<<endl;
cin>>x;
Xi[i]=x;
cout<<"enter koord toch po Y"<<i<<endl;
cin>>y;
Yi[i]=y;
}
//ввод координат искомой точки
cout<<"vvesti absc isk toch"<<endl;
cin>>isk_X;
double f, LN, XXX, XX=1.;
int j, k;

for (i=1, LN=Yi[0]; i<n; i++)//LN=Yi[0]-начало вычисления функции Pn(x)
{
XX *= (isk_X-Xi[i-1]);//вычисляем разницу между искомой Х и последующей
for (j=0,f=0; j<=i; j++)//накапливаем сумму полиномов ньютона
{
for (k=0,XXX=1.; k<=i; k++)
{
if (k!=j)
XXX *= Xi[j]-Xi[k];//знаменатель в полиноме
}
f += Yi[j]/XXX;//накапливаем сумму полиномов ньютона
}
LN += XX * f;//накапливаем значение Pn(x)
}
printf("Y[X] = %g",LN);
getch();

}

[DiGris]

Код рабочий, все в порядке