Численное дифференцирование (Лагранж)

#include <iostream>
#include <conio.h>
#include <stdio.h>
using namespace std;
 int main ()
{
system("COLOR 0A");
double ob, x1 , x2 , y1 , y2, p1 , p2 , pi , skolko , status; 
 
    system("echo Интерполяция X1 - X2 ");
    system("echo Ввести число:");
    cin >> ob; 
 
    system("echo Например 62, C1 = 60, L1 = 1.31 ,C2 = 80 , L2 = 1.29");
    cout << "C1 = ";
           cin >> x1;
    cout << "L1 = " ;
           cin >> x2;
    cout << "C2 = " ;
           cin >> y1;
    cout << "L2 = " ;  
           cin >> y2;
 
      p1 = y1 - x1;
      p2 = y2 - x2;
 
    pi = p2 / p1;
      skolko = ob - x1;
      status = x2 + (pi * skolko);
 
    cout << "_______________________________________________________________";
    cout << endl <<  ob << " = " << status;
    _getch ();
    system("CLS");
 
    return 0;
}

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <conio.h>
 
void main()
    {
    int n;
    double l, L, r ;
    int j , k ;
    float X=1.33;
    float x[6]={0.5, 0.70, 0.90, 1.10, 1.30, 1.50};
    double F[6]={0.8775826, 0.7648422, 0.6216100, 0.4535961, 0.2674988, 0.0707372} ;     // zna4eni9 funkcii iz tabl.
   int clrscr();
        std::cout<<"Istinnoe znachenie RAVNO: ";
        printf ( "%11.7f\n", cos(X)) ;
 
        for (n=2 ; n<6 ; n++)
        {
             L = 0 ;
            for (k = 0 ; k<n ; k++)
            {
              l=1;
                for (j = 0 ; j<n  ; j++)
                {
                    if (k!=j) { l = l *   (X - x[j]) / (x[k] - x[j]);}
                }
                L = L + (l * F[k]);
            }
 
 
            r= fabs (cos(X) - L);
            std::cout << "Mnogochlen Lagrandza " << n <<" ctepeni v tochke X RAVEN: ";
            printf ( "|%10.7f|\n", L);
            std::cout << "Pogreshnost ravna: ";
            printf ( "|%10.7f|\n\n", r);
 
 
        }
      getch();
        return;}