Не могу найти ошибку

[3K@BH]

Нужно написать на С++ решения линейных алгебраических уравнений комбинированным методом хорд и касательных.
Вот мой код программы:
#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;

double v( double p ) {
return 3*pow(p,2)-4*p-9;
}//первая производная функции
double z( double w ) {
return 6*w-4;
}//вторая производная функции
double f( double y ) {
return pow(y,3)+2*pow(y,2)-9*y-18;
}//уравнение

int main() {
int n=0;
double a,b,eps,c(0);
cout<<"a="; cin>>a;
cout<<"b="; cin>>b;
cout<<"eps="; cin>>eps;
do {

a=a-(f(a)/(f(b)-f(a)))*(b-a);
b=b-(f(b)/v(b));
n++;
}
while(fabs(a-b)>=eps);

c = ((a+b)/2);

cout<< "n=" << n <<endl;
cout<< "x=" << c <<endl;

return 0;
}

проблема заключается в том, что в некоторых случаях результат выводится неправильно...например взяв отрезок (2;4) мне выдает корень "-3". Может я конечно не совсем понял суть этого метода, но проблема кажется всё-таки в коде.
Заранее, спасибо.

[mrChester]

Не вникал в суть метода, но сразу бросается в глаза, что вы не используете функцию нахождения второй производной, хотя описали её, в остальном похоже на правду, в коде ошибок не заметно