Пользуясь интерполяционный полином Лагранжа, вычислить приближенное значение функции

#include <iostream>
#include "stdafx.h"
 
using namespace std;
 
 void main()
{
 
float x[7] = {1,1.02,1.05,1.08,1.16,1.19,1.2};
float y[7] = {0,0.019,0.049,0.086,0.122,0.157,0.182}; 
int i,j;
int n=8;
float f1,f2,L;
float X=1.137;
 
for (i=0;i<n;i++)
{
    f1=1;
    f2=1;
    L=0;
        for(j=0;j<n;j++)
            if(i!=j)
            {
                f1*=X-x[j];
                f2*=x[i]-x[j];
            }
            
L+=y[i]*f1/f2;
cout<<L<<endl;
return L;
}
 
 }

#include "stdafx.h"
#include <iostream>
 
using namespace std;
 
 int main()
{
 
float x[] = {1,1.02,1.05,1.08,1.16,1.19,1.2};
float y[] = {0,0.019,0.049,0.086,0.122,0.157,0.182}; 
int i,j;
int n=1;
float f1,f2,L;
 
float X=1.137;
 
for (i=0;i<n+1;i++)
{
        f1=1;
        f2=1;
        L=0;
        for(j=0;j<n+1;j++)
            if(i!=j)
            {
                f1*=X-x[j];
                f2*=x[i]-x[j];
            }
            
L+=y[i]*f1/f2;
cout<<L<<endl;
 
}
return L;
 }

cout<<L<<endl;
return L;

 void main()
{
 
    float x[7] = {1,1.02,1.05,1.08,1.16,1.19,1.2};
    float y[7] = {0,0.019,0.049,0.086,0.122,0.157,0.182};
    int i,j;
    int n=7;
    float f1,f2,L;
    float X=1.137;
 
    for (i=0;i<n;i++)
    {
        f1=1;
        f2=1;
        L=0;
            for(j=0;j<n;j++)
                if(i!=j)
                {
                    f1*=X-x[j];
                    f2*=x[i]-x[j];
                }
 
    L+=y[i]*f1/f2;
 
    }
    cout<<L<<endl;
    system("pause");
}