Программа решающая кубическое уравнение

#include <cmath>
#include <iostream>
using namespace std;
 
double P(int n, double * a, double x)
{
    double res = a[0];
    for(int i = 1; i <= n; i++)
        res = a[i]*pow(x,1.0*i);
    return res;
}
 
double D(int n, double * a, double x)
{
    
    double dx  = pow(0.1,6.0);
    return (P(n, a, x + dx) - P(n, a, x))/dx;
}
 
int main()
{
    int n;//Содержит степень полинома
    cout<<"\tSIMPLE ITERATION FOR SOLVING POLYNOMS Pn(x) = 0\n";
    cout<<"Enter power of polinom : ";cin>>n;
    int i;//Переменна счётчик
    double * a = new double[n + 1];//Выделяем память под коэффициенты полинома
    for(i = 0; i <=n; i++)
    {
        cout<<"a"<<i<<" = ";cin>>a[i];
    }
    double x;//Будет содержать корень уравнения
    cout<<"Enter initial value of x   :  ";cin>>x;
    double e;
    cout<<"Enter error of calculations : ";cin>>e;
    i = 0;
    for(i = 0; (e < fabs(P(n,a,x)/D(n,a,x))); i++)
    {
        cout<<"Iteration #"<<i + 1<<endl;
        cout<<"    x = "<<x<<endl;
        cout<<"Pn(x) = "<<P(n,a,x)<<endl;
        x = x - P(n,a,x)/D(n,a,x);
    }
    cout<<"ANSWER\n";
    cout<<"    x = "<<x<<endl;
    cout<<"Pn(x) = "<<P(n,a,x)<<endl;
    system("pause");
    return 0;
}