Вычислить интерполяционный полином Лагранжа - C++ - Ответ 14036734

#include "stdafx.h"
#include "pch.h"
#include <iostream>
 
using namespace std;
 
double f(double x)
 
{
 
    return ((3 * sin(x + 1)) / (5 - x));
}
 
 
int main()
 
{
 
    setlocale(LC_ALL, "ru");
 
    double Xi; int k;
 
    cout << "x = "; cin >> Xi;
 
    double L = 0, P, y[12], x[12] = { 0, 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4, 1.6, 1.8, 2, 2.2 };
 
    for (int i = 0; i < 12; i++)
 
        y[i] = f(x[i]);
 
    if (Xi > 2) k = 10;
 
    else
 
    {
 
        for (int i = 0; i < 12; i++)
 
            if (fabs(Xi - x[i]) <= 0.2) k = i;
 
    }
 
    for (int i = k - 1; i <= k + 1; i++)
 
    {
 
        P = 1;
 
        for (int j = k - 1; j <= k + 1; j++)
 
            if (i != j) P = P * (Xi - x[j]) / (x[i] - x[j]);
 
        L = L + y[i] * P;
 
    }
 
    cout << "L(x) = " << L << endl;
 
    system("pause");
 
    return 0;
 
}