Newton method for solving roots of polynomials

  Give me an initial value , Solving roots of polynomials by Newton method ;

 

 

#include<iostream>
#include<math.h>
#include<vector>
using namespace std;

struct   Newton {
public:
 
    double ini_para;
   
    //  The highest n Time , List of coefficients , initial value 
    double resultNewton(int n,std::vector<double> coffer,double ini_para) {
        int c = 0;
      
        double x = ini_para, x2 = 0;
        while (true) {
            c++;
            double b;
            int i = 0;
            double yuanf = 0, daof = 0;
            while( i < n+1 ){
                 yuanf += coffer[i] * pow(x,n-i); //  Primitive function 
                 if(i<n)
                     daof += (n-i)*coffer[i]*pow(x,n-i-1); // Derivative function 
                 i++;
            }
            double chu = yuanf / daof; // Primitive function / Derivative function 
            x2 = x - chu;
          
            if (abs(x2 - x) < 0.00001)
                break;
            x = x2;
        }
        int i = 0;
        double reslt = 0;
        while (i < n + 1) {
             reslt += coffer[i] * pow(x2, n - i); //  Primitive function  
            i++;
        }


        cout << " Total iterations ：" << c << "  "<< reslt;

        return x2;
    }
};


int main() {
   
    std::vector<double> coff = { 1,0 ,0,0,5,6};
    Newton newton;
    int a = newton.resultNewton(coff.size()-1, coff,-10);
    cout << "\n The result is ：\n" << a << endl;
    system("pause");
    return 0;
}