Восходящий интерполяционный многочлен Ньютона. Может пригодится

#include <conio.h>
#include <iostream>
 
using namespace std;
 
    double a;
    double b;
    double n;
double X[200], Y[200], A[200];
    int i, j, h;
    double  N=0,y;
 
double f(double x) {
 
    return   log10(x)-(7/(2*x+6));
}
 
double R( int i, int k) {
    if (k>1) 
         (y= R(i+1, k-1) - R(i, k-1));
    else (y= f(X[i+1]-f(X[0])));
    return y;
}
 
int main () {
 
    cout<<"vvedite a"<<endl;
    cin>>a;
    cout<<"vvedite b"<<endl;
    cin>>b;
    cout<<"vvedite n"<<endl;
    cin>>n;
 
    h=(b-a)/n;
    cout<<endl<<"h = " << h << endl<<endl;
 
    X[0]=a; Y[0]=f(a);
 
    //таблица
    for(i=1; i<n; i++)
        { 
            X[i]=X[i-1]+h;
            Y[i]=f(X[i]);
            cout<< "x = " << X[i] << endl << "y = " << Y[i]<<endl;
        }
 
    A[0]=1;
    for(i=1; i<n; i++) {
        A[i]=((h-(i-1))/i)*A[i-1];
    }
 
    for(i=0;i<n;i++) {
        N=N + (A[i]*R(0,i));
    }
 
    cout<<endl<<"N = " << N;
 
system ("pause");
    return 0;
}

#include <conio.h>
#include <iostream>
#include <math.h>
#include <cstdlib>
#include <string>
#include <stdio.h>
#include<fstream>
 
using namespace std;
 
    double a;
    double b;
    double n,m;
    double h;
    double t;
double X[200], Y[200], A[200], X2[200];
    int i, j;
    double  p=0,y;
 
double f(double x) {
 
    return   log10(x)-(7/(2*x+6));
}
 
double R( int i, int k) {
    if (k>1) 
         (y= R(i+1, k-1) - R(i, k-1));
    else (y= f(X[i+1]-f(X[i])));
    return y;
}
 
double N(double m) {
    p=Y[0]; A[0]=1;
        for(int i=1;i<m;i++) {
            
        A[i]=((t-i+1)/i)*A[i-1];    
         p = p + (A[i]*R(0,i));
         
         }
    return p;
}
 
int main () {
 
    cout<<"vvedite a"<<endl;
    cin>>a;
    cout<<"vvedite b"<<endl;
    cin>>b;
    cout<<"vvedite n"<<endl;
    cin>>n;
 
    h=(b-a)/n;
    cout<<endl<<"h = " << h << endl<<endl;
 
    X[0]=a; Y[0]=f(a);
 
    //таблица
    for(i=1; i<n; i++)
        { 
            X[i]=X[i-1]+h;
            Y[i]=f(X[i]);
            cout<< "x = " << X[i] << endl << "y = " << Y[i]<<endl;
        }
 
cout<<endl<<"________________________"<<endl<<"vvedite kol-vo"<<endl;
    cin>>n;
    double h1=(b-a)/n;
    cout <<endl<<"h="<<h<<endl;
    X2[0]=a;
        for(i=1; i<n; i++)
        { 
            X2[i]=X2[i-1]+h1;
            t=(X2[i]-a)/h;
            N(t);
    
        cout<<endl<<"X"<<i<<":" << X2[i] << "   " <<"N"<<i<<":"<<N(t)<<endl;
    }
system ("pause");
    return 0;
}