Метод Ньютона (касательных), программа выдает неправильный результат

[Alexandr-]

Здравствуйте. Проблема в том, что с моей функцией программа выдает неверный результат(хотя если написать функцию, например x*x-2, то все нормально).
Метод Ньютона(касательных). Язык C++.

#include "stdafx.h"
#include <math.h>
#include <iostream>
using namespace std;
 
//Функция
double f(double x) {
 
    return (cos(x)*sin(6*x*x)-4*sin(x));
}
 //Производная
double f1(double x) {
 
    return (12*x*cos(x)*cos(6*x*x)-sin(x)*sin(6*x*x)-4*cos(x));
}
 //2 производная
double f2(double x) {
 
    return (-144*x*x*sin(6*x*x)*cos(x)-24*x*sin(x)*cos(6*x*x)+4*sin(x)-sin(6*x*x)*cos(x)+12*cos(x)*cos(6*x*x));
}
 
int main() {
    int n=0;
    double a,b,c,eps;
 
    cout<<"Vvedite granici otrezka\n";
    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;
}

[olej.tsil]

Сообщение от Alexandr-

с моей функцией программа выдает неверный результа

Бывает

[Alexandr-]

Код вроде правильный, так почему выдает неверный ответ... Я вставлял свою функцию в другие проги и там тоже самое.

[type_Oleg]

Это наверное недостаток метода. Он же не всегда работает. Вот - ссыла
Там еще пример приведен f(x)=x^3-2*x+2 - если начинать с 0, то все зациклится.

Тем более у вас такая функция со сложной производной: