При вычислении нелинейного уравнения методом Ньютона выдает неверное значение

#include <iostream>
#include <math.h>
using namespace std;
double Funkt(double x){
return x*x-log(x)+2*cos(x)-1;}
 
double Koren(double x0){
double x,pr,h=0.001,x1=x0;
 
do{
    x=x1;
    pr=2*x+2*sin(x)-1/x;
    x1=x-Funkt(x)/pr;
 
}
while(x1-x<h);
 
return x;}
 
 
int main(){
 
double x0,x;
cout << "Nachalnoe priblizenie x_0 = ";
cin >> x0; //x=4
 
x=Koren(x0);
cout << "x = " << x << endl;
return 0;}

while(fabs(x1-x)>h);

#include <iostream>
#include <math.h>
using namespace std;
double f(double x)
{
    return x*x-log(x)+2.*cos(x)-3.;
}
double f1(double x) // 1 производная
{
    return 2.*x-1./x-2.*sin(x);
}
double f2(double x) // 2 производная
{
    return 2.+1./(x*x)-2.*cos(x);
}
double Koren(double x0){
    double x,h=0.001,x1=x0,y;
    do{
        if(f(x1)*f2(x1)<=0)
        {
            cout<<"Can't to solve\n";
            return 0;
        }
        y=f(x1);
        x=x1;
        x1=x-y/f1(x);
    }while(fabs(x1-x)>h || fabs(f(x1))>h);
 
    return x;
}
 
 
int main(){
 
    double x0,x;
    cout << "x_0 = ";
    cin >> x0;
 
    x=Koren(x0);
    if(x!=0)
    {
        cout << "x = " << x << endl;
        cout << "f(x) = " << f(x) << endl;
    }
    return 0;
}