Найдите ошибки в программе, решающей уравнения методом биссекции-Ньютона

#include <math.h> 
#include <cmath> 
#include <stdio.h> 
#include <conio.h> 
#include <iostream> 
#include <stdlib.h> 
using namespace std; 
int getch (void); //Объявление прототипов функций 
int i; 
//Функции-уравнения, решаемые в задаче 
double func1(double x) { return exp(-x) + pow(x, 2) + 2; } 
double func2(double x) { return -pow((x-2), 2); } 
double func3(double x) { return exp(x)-2*x-2; } 
double func4(double x) { return sin(x) + x - 1; } 
double func5(double x) { return exp(-pow(x,2));} 
double func6(double x) { return 2*exp(-3*x) - x +1; } 
double func7(double x) { return log(x) - x + 1; } 
double func8(double x) { return log(x) - pow((x-1), 2); } 
//Первая производная функции 
double Derivative(double x, double h, double(*func)(double x) ) 
{ 
double Derivative = ((func(x + h) - func(x - h))/(2*h)); 
return Derivative; 
} 
//Вторая производная функции 
double secondDerivative(double x, double h, double(*Derivative)(double x) ) 
{ 
double secondDerivative = ((Derivative(x + h) - Derivative(x - h))/(2*h)); 
return secondDerivative; 
} 
double secondDerivative(double, double, double(*func)(double x)); 
double ChordLocal(double, double (*func)(double x)); 
double Newton (double, double ,int,double , double (*func)(double x), double (*Derivative)(double, double, double(*func)(double x)), double (*secondDerivative)(double, double, double(*func)(double x))); 
double func1(double), func2(double), func3(double), func4(double), func5(double), func6(double), func7(double), func8(double); 
int main() 
{ 
int maxIteration, j, funcNum; 
int exit = 0;// переменные для выхода и цикла 
double x0,xn;// вычисляемые приближения для корня 
double a, b, c, h, eps, func;// границы отрезка и необходимая точность 
//Объявление массива ссылок на функции 
double (*mas[8])(double) = {func1, func2, func3, func4, func5, func6, func7, func8}; 
do 
{ 
printf ("Vvedite eps"); 
scanf ("%lf",&eps); // вводим нужную точность вычислений 
printf ("Vvedite h"); 
scanf ("%lf",&h); 
printf ("Vvedite maximum number of iteration: "); 
scanf ("%d",&maxIteration); 
for(funcNum = 0; funcNum < 8; funcNum++) 
{ 
printf("%d\n", funcNum); 
if (mas[funcNum](a)*secondDerivative(a,h,mas[funcNum](a))>0) x0 = a; // для выбора начальной точки проверяем f(x0)*p2f(x0)>0 ?
else x0 = b; 
c=(a+b)/2; 
if (mas[funcNum](c)*mas[funcNum](a)<0) 
b=c; 
if (mas[funcNum](c)*mas[funcNum](b)<0) 
a=c; 
xn=Newton (a, b, maxIteration, eps, mas[funcNum], *Derivative, *secondDerivative); 
xn = x0-mas[funcNum](x0)/Derivative(x0,h,mas[funcNum](x0)); // считаем первое приближение 
while (fabs(x0-xn)>eps) 
{ 
x0=xn; 
xn=Newton (a, b, maxIteration, eps, mas[funcNum], *Derivative, *secondDerivative); 
xn = x0-mas[funcNum](x0)/pf; // непосредственно формула Ньютона 
} 
while( xn < 0 ) // пока не достигнем необходимой точности, будет продолжать вычислять 
{ 
x0+=h; 
xn = x0-mas[funcNum](x0)/pf(x0); // непосредственно формула Ньютона 
while (fabs(x0-xn)>eps) 
{ 
x0=xn; 
xn=Newton (a, b, maxIteration, eps, mas[funcNum], *Derivative, *secondDerivative); 
xn = x0-mas[funcNum](x0)/pf(x0); // непосредственно формула Ньютона 
} 
printf("xn = %.9f\n",xn); 
printf("f(x0) = %.9f\n",mas[funcNum](x0)); 
printf("Number of iteration = %d\n\n",i); 
} 
func=exp(-xn)+xn*xn-2; 
cout<<"\nProverka kornya = exp(-x)+x*x-2 = "<<func;
cout<<"\nExit?=>";
cin>>exit;
} 
} 
while (exit!=1); // пока пользователь не ввел exit = 1 
return 0; 
}