Разработать программу "Численное дифференцирование".

double f(double x)
{
return sin(x);
}
 
double diff(double x)
{
const eps=1E-10;
return (f(x+eps)-f(x))/eps;
}

#include <iostream.h>
#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.
    clrscr();
        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);
            cout << "Mnogochlen Lagrandza " << n <<" ctepeni v tochke X RAVEN: ";
            printf ( "|%10.7f|\n", L);
            cout << "Pogreshnost ravna: ";
            printf ( "|%10.7f|\n\n", r);
 
 
        }
      getch();
        return;}

#include <iostream>
#include <cstdlib>
#include <cmath>
#define PI 3.14159
 
using namespace std;
double dx = 0.001;
 
double f(double X) // Сама Функция(в моём случае для экспоненты)
{
    return exp(X);
}
 
double P(double X) // Сама Функция ( в моём случае для синуса)
{
    return sin(X);
}
 
double diff(double (*p)(double), double x) // Первым аргументом даю возможность передачи функции
{
    double Dx;
    Dx = (p(x + dx) - p(x)) / dx; // Данная формула дифференцирования
    return Dx;
}
 
 
 
int main() {
    setlocale(LC_ALL, "Russian");
    double X, Print, func;
    cout << "  Что вы хотите найти Sin(x) или Exp(x)?(напишите 1 или 2 соответственно): "; cin >> func;
    
    if (func == 1) {
        cout << "  Введите Х для Sin(Х):  "; cin >> X; //вводим значение PI/6
        Print = diff(&P, X);// Передаю ссылку на первую функцию f, функции diff
        cout << endl << "  Для Х = " << X << "\n  sin(x) = " << Print;
    }
 
    if (func == 2) {
        cout << "  Введите Х для exp(Х):  "; cin >> X; //вводим значение 1.5
        Print = diff(&f, X);// Передаю ссылку на первую функцию f, функции diff
        cout << endl << "  Для Х = " << X << "\n  exp(x) = " << Print;
    }
    
    system("pause > null");
    return 0;
}