Полином Лагранжа - Сделать интерполяцию по точкам - C++ - Ответ 7567288

//Lagrange
double L(double xc, double x[],double y[]) {//хс точка в которой щитаеца всё
double Ch; 
double Zn;
 int k;
 double R=0;
 for (int i = 0; i < n; i++) { 
 
 
Ch = 1; Zn = 1;
 
for (k = 0; k < n; k++ ) {
 
if ( k == i ) continue; 
 
Ch *= xc - x[k];
}
for(k= 0; k < n;k++) {
 
if (x[i] == x[k]) continue; 
Zn *= x[i] - x[k];
}
 R += y[i]*Ch/Zn;
}
 
return R; 
 }
void main(void){
double Xi[n],Yi[n];
FILE *fp = fopen("test.txt","r");
FILE *fp1 = fopen("testy.dat","w");
 
int nf,i=0;
 
fscanf(fp,"%d",&nf);
    
    while (!feof(fp))
    {
        fscanf(fp,"%lf%lf",&Xi[i],&Yi[i]);
        i++;
        if (i==n+1) 
        break;
    }
    for(int i=0;i<n;i++)
 {
     printf("%f  %f\n",Xi[i],Yi[i]);}   
 
 
double h=-0.55;
 for(int i=0;i<21;i++)
 {h+=0.05;
     printf("=============%f  %f\n",L(h,Xi,Yi),h);
fprintf(fp1,"%f %f %f %f\n",L(h,Xi,Yi));}
 
fclose(fp);
fclose(fp1);
 getch();
 
}