Решение уравнения методом Ньютона

double f (double x,double a,double b,double c,int k)
{
    switch(k)
    {
    case 1: return a*pow(x,3)+b*x+c;
    case 2: return a*sin(x+b)+c;
    case 3: if (x+b!= 0) return a*log(fabs(x+b))+c;
            else cout<<"Error. Natural logarithm is not defined!"<<endl;
        return 0;
    }
    return 0;
}
 
double df (double x,double a,double b,int k)
 
{
    switch (k)
    {
    case 1: return 3*a*pow(x,2)+b;;
    case 2: return a*cos(x+b);
    case 3:
        if(x!=-b)
        {
        if(x + b > 0)
        return a / fabs(x + b);
        else
        return - a / fabs(x + b);
        }
        else
        {
        cout<<"Error. Denominator takes a zero value!"<<endl;
        return 0;
        }
    }
    return 0;
}
void newton(double a, double b, double c, double eps, int k, double x0) {
    double xk=0;
    int count=1;
    do
    {
        if (count>1)  x0 = xk;
        xk   = x0 - ( (f(x0,a,b,c,k))/( df(x0,a,b,k) )) ;
        cout << count <<" -iter, x = " << xk << endl;
        count++;
    }
    while ((fabs(f(xk,a,b,c,k)) > eps) && (count<1000));
    output(xk, f(xk,a,b,c,k));
}