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

#include "pch.h"
#include <iostream>
#include <cmath>
using namespace std;
double y, f, fa, fb, k, x, e;
int main()
{
    setlocale(LC_ALL, "rus");
    cout << "Веедите x\n" << "x= ";
    cin >> x;
    cout << "Введите e\n" << "e= ";
    cin >> e;
    f = cos(x) - 3 * x + 1;
    fa = -sin(x) - 3;
    fb = cos(x);
    for (k = 0;abs(f) < e;k = k++)
    {
        y = f - f / fa;
    }
    cout << "y= " << y <<"\n";
    system("pause");
    return 0;
}

#include <iostream>
#include <cmath>
 
using namespace std;
 
const int N = 10;
 
//Функция для полинома:
double F(double x)
{
 return cos(x) - 3 * x + 1;
 }
//Функция поиска корня:
 
double FindRoot(double (*f)(double),double x0,int n)
{
 double x = x0,df,h = 0.00001;
 
 df=(f(x + h) - f(x)) / h;
 
 for(int i = 1; i <= n;i++)
 {
  x = x - f(x) / df;
 }
 
return x;
}
int main()
{
//Начальное приближение и корень:
double x0,x;
cout << "initial x0 = ";
cin >> x0;
//Поиск решения:
x = FindRoot(F,x0,N);
cout << "x = " << x << endl;
return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
double  f, fa, x, e;
int main()
{
setlocale(LC_ALL, "rus");
cout << "Веедите x\n" << "x= ";
cin >> x;
cout << "Введите e\n" << "e= ";
cin >> e;
//x=1;
//e=0.00001;
f = cos(x) - 3 * x + 1;
while(abs(f)>e)
{
    f = cos(x) - 3 * x + 1;
    fa = -sin(x) - 3;
    x = x - f / fa;
cout << "x= " << x <<"\n";
}
cout << "Конечное решение х= " << x <<"\n";
system("pause");
return 0;
}