Вместо числа выводится строка "-1.#IND"

// X*x+4sin(x)=0
#include <iostream>
#include <math.h>
#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
#include <iomanip>
using namespace std;
#define  db double
const db c=1, d=2, e=0.001;
/*------------------------------------*///вычисление значения функции
db f(db x, db p)
{
    return pow((db)x,2)-10-pow(p,pow((db)x-2,0.5));
}
/*------------------------------------*/
db f1 (db x,db p)   // 1-z производная 
{
    return 2*(db)x-pow(p,(db)x)*log(p);
}
/*------------------------------------*/
db f2 (db x,db p)     // 2-я производная 
{
    
    return 2-pow(p,(db)x)*pow(log(p),2);
}
/*------------------------------------*/
db the_chord_method(db left_bound, db right_bound, db p)//метод хорд
{
 
    return (left_bound-(f(left_bound,p))/(f(right_bound,p)-f(left_bound,p))*(right_bound-left_bound));
}
/*------------------------------------*/
db method_of_tangents(db left_bound, db right_bound, db p)//метод касательных
{
    return right_bound-(f(right_bound,p))/(f1(right_bound,p));
}
/*------------------------------------*/
void P_return(double c, double d, int m, db*P)//вычисление коээфициентов Р
{
    db h=(d-c)/m;
 
    for(int i=0;i<m+1;i++)
    {   P[i]=0;
        P[i]=(c+i*h);       
    }   
}
/*------------------------------------*/
void print(db *X,db *P, db *left_bound, db* right_bound, int m)
{
     cout<<"\n\t\t__________________________________________"<<endl;
     cout<<"\t\t"<<" N"<<setw(7)<<"a"<<setw(10)<<"b"<<setw(10)<<"P"<<setw(10)<<"X"<<endl;
     cout<<"\t\t__________________________________________"<<endl;
     for(int i=0;i<m+1;i++)
     {
        cout<<"\t"<<setw(10)<<i+1<<setw(10)<<setprecision(5)<<left_bound[i]<<setw(10)<<setprecision(5)<<right_bound[i];
        cout<<setw(10)<<setprecision(5)<<P[i]<<setw(10)<<setprecision(5)<<X[i]<<endl;   
     }
        cout<<"\t\t__________________________________________"<<endl;
}
/*------------------------------------*/
void find_values_X(db *X,db *P, db *left_bound, db* right_bound, int m)//комбинированный метод
{
     for(int i=0;i<m+1;i++)
     {
            if ((f1(c,P[i])*f1(d,P[i])<0) || (f2(c,P[i])*f2(d,P[i])<0)) 
                exit(0);//если производные имеют разные знаки
                    if (f(c,P[i])*f2(c,P[i])>0)
                    {
                        right_bound[i]=c; 
                        left_bound[i]=d;
                    }    
                    else 
                        if (f(d,P[i])*f2(d,P[i])>0)
                        {
                            right_bound[i]=d; 
                            left_bound[i]=c;
                        }
 
                    bool OK=false;
                    int n=0; 
         
                        while ((!OK) && (n<m))
                        {
                            right_bound[i]=method_of_tangents(left_bound[i], right_bound[i],P[i]);
                            left_bound[i]=the_chord_method(left_bound[i], right_bound[i],P[i]);
                            X[i]=(right_bound[i]+left_bound[i])/2;
                                if (fabs(right_bound[i]-left_bound[i])<e)
                                    OK=true;
                            n++;
                        }
     }
}
/*------------------------------------*/
void pre_main()
{    
    system("cls");
    setlocale(LC_ALL,"");
    cout<<"\t\t\tРезультаты работы программы:"<<endl;
    int m=15;
    db *left_bound, *right_bound, *X;
    X=new db[m+1];
    left_bound=new db[m+1];
    right_bound=new db[m+1];
    db *P= new db[m+1];
    for(int i=0;i<m;i++)
    {
            X[i]=0;
            left_bound[i]=0;
            right_bound[i]=0;
    }    
    P_return(c,d,m,P);  
    find_values_X(X,P,left_bound,right_bound,m);
    print(X,P,left_bound,right_bound,m);
}
/*-------------------------------------------------*/
int main ()
{
     pre_main();
     _getch();
     return 0;
}