Методы оптимизации. Метод квадратичной интерполяции

[krasy]

Нужно написать программу, которая находит наименьшее значение функции f(x)=(x-3)*(x-3)/(x-2) на промежутке [2.1;5] методом парабол.
я сделала что-то, но думаю я где-то облажалсь. помогите найти ошибку, пожалуйста.

#include<iostream>
#include<conio.h>
#include<math.h>
using namespace std;
double a0,a1,a2,xx0,xx1,x1=2.1,x2,x3=5,e,dx;

double f(double x) {double f;
	f=(x-3)*(x-3)/(x-2);
	return f;
}

double a();
void input();

int main()
{input();
return 0;}

double a()
{a0=f(x1);
 a1=(f(x2)-f(x1))/(x2-x1);
 a2=((f(x3)-f(x1))/(x3-x1)-a1)/(x3-x2);
return a0,a1,a2;}

void input()
{
cout<<"e=";cin>>e;
cout<<"dx=";cin>>dx;
x2=x1+dx;

do
{ xx0=(x1+x2-a1/a2)/2;
  if(f(x1)>f(xx0)&&f(xx0)<f(x2)&&x1<xx0&&xx0<x2) {x1=xx0;x2=x2;x3=x3;} else
  if(f(x2)>f(xx0)&&f(xx0)<f(x3)&&x2<xx0&&xx0<x3) {x1=x2;x2=xx0;x3=x3;}
  a();
  xx1=(x1+x2-a1/a2)/2;  
}
while(fabs(xx0-xx1)<e);

cout<<"x*="<<xx1<<"; f(x*)="<<f(xx1);

getch();
}