Не могу осилить метод Ньютона

#include <conio.h>
#include <math.h>
#include <iostream.h>
#define pi 3.14
 
double f(double x) {
 
    return   x*x-(cos(pi*x));
}
 
double f1(double x) {
 
    return   2*x+(1/x);
}
 
double f2(double x) {
 
    return   2+(-1/(x*x));
}
 
int main() {
    int n=0;
    double a,b,c,eps;
    cout<<"a="; cin>>a;
    cout<<"b="; cin>>b;
    cout<<"eps="; cin>>eps;
    if(f(a)*f2(a)>0) c=a;
    else c=b;
    do {
        c=c-f(c)/f1(c);
        n+=1;
    }
    while (fabs(f(c))>=eps);
        cout<<"c="<<c<<"\n";
        cout<<"n="<<n<<"\n";
        getch();
    return 0;
}

#include <conio.h>
#include <math.h>
#include <iostream>
#define pi 3.14
using namespace std;
 
double f(double x) {
 
    return   acos(x) - sqrt(1.0 - 0.3 * x * x * x);
}
 
double f1(double x) {
 
    return  0.45 * x * x / sqrt(1 - 0.3 * x * x * x) - 1 / sqrt(1 - x * x);
}
 
double f2(double x) {
 
    return   0.9 / (1 - 0.3*pow(x, 3)) - 3 / 2 * (pow(x, 2) + 1) + 0.3*pow(x, 4) / 1 - 0.3*pow(x, 3);
}
 
int main() {
    int n = 0;
    double a, b, c, eps;
    cout << "a="; cin >> a;
    cout << "b="; cin >> b;
    cout << "eps="; cin >> eps;
    if (f(a)*f2(a)>0) c = a;
    else c = b;
    do {
        c = c - f(c) / f1(c);
        n += 1;
    } while (fabs(f(c)) >= eps);
    cout << "c=" << c << "\n";
    cout << "n=" << n << "\n";
    
    system("pause");
    return 0;
}