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

[lullaby]

Уравнение задаётся в программе в виде функции f(x)=0. Задаётся случайное начальное значение x0. Первое приближенное значение считается по методу начальных приближений.Все последующие значения считаются по формуле Вегстейна до тех пор, пока разница между нынешним и предыдущим значениями найденных x не станет меньше заданной точности.
Помогите пожалуйста найти ошибку!

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define formula(x) (x*x)-(3*x)+2

//Функция, высчитывающая функцию по методу Вегстейна
float vegstaine(float x0,float x1)
{

float f=(formula(x1)-x1)/(x1-x0);
return x1+(1/(1-f))*(formula(x1)-x1);

}

int main()
{
//x0-начальное приближение, x-промежуточный и конечный результат,e-погрешность tmp-хранение предыдущих результатов цикла

float x0,x,e,tmp;
scanf("%f %f",&x0,&e);
x=formula(x0);//Первая итерация
do
{

tmp=x;
x=vegstaine(x0,x);
x0=tmp;

}while(fabs(x-x0)>0);
printf("x=%f ",x);
getch();

}

[EUGY]

e-погрешность

Где она должна быть в вычислениях?
while ( fabs(x-x0) > e ); // тут ?

[lullaby]

Сообщение от EUGY

Где она должна быть в вычислениях?
while ( fabs(x-x0) > e ); // тут ?

Да. Нужно вычислять значения по формуле Вегстейна до тех пор, пока разница между нынешним и предыдущим значениями найденных x не станет меньше заданной точности.
Хотя вполне возможно, что можно сделать как-то подругому? Т.к. разница меньше точности может найтись, но вот значение при этом корню исходной функции не соответствует.