Квадратичная интерполяция

#include <iostream>
#include <cmath>
#include <math.h>
using namespace std;
double f(double x)
{
    double ret = 2* pow(x, 2) + 16/x;
    return ret;
}
int main()
{
    double  help, p,p1, x0, x1, x2, h, e1, e2,  index, tmp, minintpo,  help1;
    double mas[2];
    bool sled = true;
    bool proverka = true;
    double x[2];
    cout << "Enter x1 : "; cin >> x[0];
    cout << "Enter h : "; cin >> h;
    cout << "Enter e1 : "; cin >> e1;
    cout << "Enter e2: "; cin >> e2;
    do
    {
        if (sled == true)
        {
            if (proverka == true)
            {
                x[1] = x[0] + h;
            }
            else
            {
                x[1] = minintpo;
 
            }
            x0 = x[0];
            x1 = x[1];
            if (f(x0) > f(x1))
            {
                cout << "f(x0)>f(x1" << endl;
                x[2] = x0 + 2 * h;
            }
            if (f(x0) <= f(x1))
            {
                x[2] = x0 - h;
            }
            x2 = x[2];
        }
            mas[0] = f(x0);
            mas[1] = f(x1);
            mas[2] = f(x2);
            tmp = mas[0];
            for (int i = 0; i<2; i++)
            {
                if (tmp>mas[i+1])
                {
                    tmp = mas[i + 1];
                    index = x[i + 1];
                    
                }
            }
        cout << "xmin" << index << endl;
        cout << "Fmin" << tmp << endl;
        minintpo = 0.5*(((x1*x1) - (x2*x2))*mas[0] + ((x2*x2) - (x0*x0))*mas[1] + ((x0*x0) - (x1*x1))*mas[2]) / ((x1 - x2)*mas[0] + (x2 - x0)*mas[1] + (x0 - x1)*mas[2]);
        help = fabs((tmp - f(minintpo)) / (f(minintpo)));
        help1 = fabs((index - minintpo) / (minintpo));
        cout << "" << help << endl;
        cout << "" << help1 << endl;
        if (minintpo >= x0 && minintpo <= x2)
        {
            sled = false;
            cout << "yes" << endl;
            if (index > minintpo)
            {
 
                x[1] = minintpo;
            }
            else
            {
                x[1] = index;
            }
            p = modf(x[1], &p1);
            x[0] = x[1] - p;
            x[2] = x[1] + (1 - p);
            cout << "x0" << x[0] << endl;
            cout << "x2" << x[2] << endl;
 
        }
            if (minintpo <= x0 && minintpo >= x2)
            {
                sled = true;
                proverka = false;
                x[1] = minintpo;
 
            }
        } while ((help > e1) && (help1 > e2));
        cout << "tmp" << tmp << endl;
        cout << "index" << index << endl;
        cout << "minintpo" << minintpo << endl;
        cout << "f(interpolinom)" << f(minintpo);
    }