Найти производную высшего порядка

#include<iostream>
#include<cmath>
#include<math.h>
using namespace std;
double func(double x)
{
    return pow(x,3)-x-pow(x,2)+1;
}
int main()
{
    double a,b,,f1,f2,f4,f3,h=0.25;
    setlocale(LC_ALL,".1251");
    cout<<"Введи, на каком узле искать производные: "<<endl;
    cin>>a>>b;
    double i=a;
    for(a;i<=b;a+=i)
    {
        f1=((func(a+h)-func(a))/((a+h)-a));
        f2=((func(f1)-func(i+h))/((a+h*2)-a));
        f3=((func(f2)-func(f1))/((a+h*3)-a));
        f4=((func(f3)-func(f2))/((a+h*4)-a));
    
    cout<<"Значение 1 производной в точке "<<i<<" равно:"<<f1<<endl;
    cout<<"Значение 2 производной в точке "<<i<<" равно:"<<f2<<endl;
    cout<<"Значение 3 производной в точке "<<i<<" равно:"<<f3<<endl;
    cout<<"Значение 4 производной в точке "<<i<<" равно:"<<f4<<endl;
    cout<<"Значение функции в точке "<<i<<" равно:"<<func(i)<<endl;
    i+=h;
        }
    system("pause");
}