Нахождение экстремума функции

[Виктори]

Здрасте! Помогите, пожалуйста, разобраться с ошибками. Задание следующее: Найти экстремумы функции методом золотого сечения (это когда мы постепенно сужаем отрезок. на котором рассматриваем функцию). Язык программирования - Си++.
Вот что я нашла:

#include <iostream>
#define eps 1e-3
#define tay 1.618
 
void Gold(double a, double b)
{
        std::cout<<"\n\n\n\tМетод золотого сечения:\n";
        double x1, x2, _x, xf1, xf2;
        int iter(0);
        x1 = a + (b - a) / (tay * tay);
        x2 = a + (b - a) / tay;
        xf1 = Fun(x1);
        xf2 = Fun(x2);
  P:
        iter++;
        if(xf1 >= xf2)
        {
                a = x1;
                x1 = x2;
                xf1 = Fun(x2);
                x2 = a + (b - a) / tay;
                xf2 = Fun(x2);
        }
        else
        {
                b = x2;
                x2 = x1;
                xf2 = xf1;
                x1 = a + (b - a) / (tay * tay);
                xf1 = Fun(x1);
        }
        if(fabs(b - a) < eps) 
        {
                _x = (a + b) / 2;
                std::cout<<"Результат:\nx = "<<_x<<"\t\tF(x) = "<<Fun(_x)<<
                        "\nКоличество итераций: "<<iter;
        }
        else
                goto P;
}

Вопрос 1: это вообще верно? (сама проверить не могу)
Вопрос 2: Если нет...то что мне вообще надо делать?????
Заранее всем спасибо) А если поможете - будет вообще супер!

[Stilet]

По-моему у тебя дюже сложно написанно. Типичные закидоны методичек.
Я бы написал проще:

#include <math.h>

double func(double x) {return sin(x);}

int _tmain(int argc, _TCHAR* argv[])
{
	double a=0,b=90,eps=0.1,fi=0.1;
	double xa=b-(b-a)/fi,xb=a+(b-a)/fi;
	double ya=func(xa),yb=func(xb);

	for(;(a-b)<eps;){
		if(ya<=yb){
			b=xb;xb=xa;xa=b-(b-a)/fi;
			yb=ya;ya=func(xa);
		} else {
			a=xa;xa=xb;xb=a+(b-a)/fi;
			ya=yb;yb=func(xb);
		}
	}
	printf("%f",(ya<yb)?ya:yb);
	getchar();
	return 0;
}

[Виктори]

Я наверное тормоз...но можешь объяснить что тут написано???