Метод Ньютонов в переполнении стека

Привет, я проверил множество других вопросов, но не могу понять, почему моя реализация не сходится. Я продолжаю получать часть «Ошибка, нет конвергенции» в моей программе, даже когда я вхожу в корень.

Функция у = х ^ 2 — 1

Вот код:

// Newton sqaure root finder function

#include <iostream>
#include <cmath>

int main()
{
using namespace std;

// Enter an initial guess x
cout << "Enter an initial guess: ";
double x;
cin >> x;// Define & initialize the error, tolerance and iteration variables
double tol = 1e-12;
cout << 1e-12;
double error = tol + 1;
int it = 0;
int max_it = 100;

// Define the x1 variable to hold the latest result (root approximation)
double x1;

// Start while loop with guess x to find the root
while (error > tol && it < max_it)
{
x1 = x + (x*x-1) / (2*x);
error = fabs(x1 - x);
x = x1;
it++;
cout << error << endl;
}

if (error <= tol)
{
cout << "The root is " << x << endl;
}
else
{
cout << "Error, no convergence" << endl;
}

cin.get();
cin.get();

return 0;
}