Классы: найти вторую производную функции F

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
double f(double x)
{
 return exp(x) - atan(x);
}
 
int main()
{
  int iKEY;
  double x,a,b,h;
  do
        {
                cout<<"\ vedite otrezok [a;b] of tabulation\n";
                cout<<"a = ";cin>>a;x = a;
                cout<<"b = ";cin>>b;
                cout<<"\ vedite sag h of tabulation\n";
                cout<<"h = ";cin>>h;
                cout<<"\tTabulated values\n";
                cout<<"  x\t|  f(x)\n";
                while(x <= b)
                {
                        cout<<" "<<x<<"\t| "<<f(x)<<"\n";
                        x += h;
                }
                cout<<"Enter 1 for new input\n";cin>>iKEY;
        }
        while(iKEY == 1);
        return 0;
}