Ошибка при использовании метода хорд

// chord method.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include "conio.h"
#include <math.h>
#include "iostream"
#define e 2.71
 
//метод хорд
using namespace std;
 
 
double func(double x)
{
    return pow(e,sqrt(x))-1/x;
}
 
double find(double infinum, double supremum, double epsilon)
{
    while (fabs(supremum - infinum) > epsilon)
    {
        infinum = supremum - (supremum - infinum) * func(supremum) / (func(
                supremum) - func(infinum));
        supremum = infinum - (infinum - supremum) * func(infinum) / (func(
                infinum) - func(supremum));
    }
 
    return supremum;
}
 
int main()
{
    double a = -5, b = 5;
    cout << find(a, b, 0.001)<<endl;
    getch();
    return 0;
}

// chord method.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include "conio.h"
#include <math.h>
#include "iostream"
#define e 2.71
 
//метод хорд
using namespace std;
 
 
double func(double x)
{
    //return pow(e,sqrt(x))-1/x;
    return pow(e,x)-1/(x*x);
}
 
double find(double infinum, double supremum, double epsilon)
{
    while (fabs(supremum - infinum) > epsilon)
    {
        infinum = supremum - (supremum - infinum) * func(supremum) / (func(
                supremum) - func(infinum));
        supremum = infinum - (infinum - supremum) * func(infinum) / (func(
                infinum) - func(supremum));
        supremum=sqrt(supremum);
    }
 
    return supremum;
}
 
int main()
{
    double a = -5, b = 5;
    cout << find(a, b, 0.001)<<endl;
    getch();
    return 0;
}

#include <iostream>
#include <math.h>
#include <conio.h>
#include <stdio.h>
#include <conio.h>
#define e 2.71
 
using namespace std;
 
void main() {
    double a, b, x, eps, fc, fx, c; //a,b - концы отрезка, x - корень, e - эпсилон, fx/fc - значения функции в точках x и c
    int count=0;
    a=1;    
    b=2;
 
    eps=0.001;
 
    x=b;
    c=a;
 
    while(1) {
        count++;
        
        fc=pow(e,sqrt(x))-1/x;
        fx=pow(e,sqrt(x))-1/x;
 
 
        x=x-fx*(x-c)/(fx-fc);
 
        cout << "x: " << x << "\n";
        
        if(abs(fx)<eps) 
            break;
 
    }
    cout << "\nSteps: " << count;
    _getch();
 
}