Поиск точек экстремума методом золотого сечения

#include <iostream>
#include <conio.h>
#include <fstream>
#include <string>
#include <Windows.h>
#include <locale.h>
#include <vector>
#include <cmath>
 
using namespace std;
 
double F_Parabola(const double &X, const double &A, const double &B, const double &C)
{ return A * X*X + B * X + C; }
 
double Extremum(double &a, double &b, const double &tochnosty, const double &A, const double &B, const double &C)
{
    double X_M = (a + b) / 2;
    double X_L, X_R;
    do
    {
        X_L = X_M - tochnosty;
        X_R = X_M + tochnosty;
 
        if (F_Parabola(X_L, A, B, C) < F_Parabola(X_M, A, B, C) &&
            F_Parabola(X_M, A, B, C) < F_Parabola(X_R, A, B, C))
        {
            a = X_M;
            X_M = (a + b) / 2;          
        }
        else if (F_Parabola(X_L, A, B, C) > F_Parabola(X_M, A, B, C) &&
                 F_Parabola(X_M, A, B, C) > F_Parabola(X_R, A, B, C))
        {
            b = X_M;
            X_M = (a + b) / 2;      
        }
        else if (F_Parabola(X_L, A, B, C) < F_Parabola(X_M, A, B, C) &&
                 F_Parabola(X_M, A, B, C) > F_Parabola(X_R, A, B, C))   
            return X_M;
 
    }while(true);
}
 
int main()
{
    setlocale (LC_ALL, "rus");
    cout << "Введите исходные данные для уравнения вида A*x^2 + B*x + C." << endl; 
    double A, B, C;
    cout << "A: ";
    cin >> A;
    cout << "B: ";
    cin >> B;
    cout << "C: ";
    cin >> C;
    
    double a, b;
    cout << "Границы интервала (a и b): ";
    cin >> a >> b;
    
    double tochnosty;
    cout << "Точность: ";
    cin >> tochnosty;   
    double res = Extremum(a, b, tochnosty, A, B, C);
 
    int aa = fabs(log10(tochnosty)); 
    char s[100];
 
    sprintf(s, "%dlf", aa);
    for(int i = strlen(s) + 2; i >= 2; i--)
        s[i] = s[i -2];
 
    s[0] = '%';
    s[1] = '.';
 
    printf("X = ");
    printf(s, res);
    printf(", Y = %lf", F_Parabola(res, A, B, C));  
 
    _getch();
    return 0;
}

    int aa = fabs(log10(tochnosty)); 
    char s[100];
 
    sprintf(s, "%dlf", aa);
    for(int i = strlen(s) + 2; i >= 2; i--)
        s[i] = s[i -2];
 
    s[0] = '%';
    s[1] = '.';