Поиск экстремума функции двух переменных

#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;
 
float f(float l, float i)
{
    return pow((i-pow(l,2)),2)+pow((1-l),2);//Исходная функция
}
 
float pf1(float l, float i)
{
    return -4*l*(i-pow(l,2))-2*(1-l);//производная функции по х1
}
 
float pf2(float l, float i)
{
    return 2*(i-pow(l,2))+2*(1-l);//производная функции по х2
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(LC_ALL, "RUSSIAN");
    float n;
    float a = 1;//коэффициент пропорциональности
    float x1, x2, fu,df1,df2;
    //ввод начальные значения
    cout << "Введите начальное значение x1= ";
    cin>>x1;
    cout << "Введите начальное значение x2= ";
    cin>>x2;
    fu = f(x1, x2);
    cout << "fu= " << fu << endl;
    cout << "pf1= " << pf1(x1,x2) << endl;
    cout << "pf2= " << pf2(x1,x2) << endl;
    cout << "____________" << endl;
    do
    {
        a = 0.5*a;
        n = fu;
        //значения на следующем шаге
        x1 = x1 - a*pf1(x1,x2);
        x2 = x2 - a*pf2(x1,x2);
        fu = f(x1, x2);
        cout << "x1= " << x1 << endl;
        cout << "x2= " << x2 << endl;
        cout << "fu= " << fu << endl;
        cout << "pf1= " << pf1(x1,x2) << endl;
        cout << "pf2= " << pf2(x1,x2) << endl;
        cout << "____________" << endl;
        df1 = pf1(x1,x2);
        df2 = pf2(x1,x2);
        
    } 
    while((df1!=0) || (df2!=0));//цикл работает, пока производные не обнулятся
    float y = fu;
    cout << "Minimum =  " << y << " v tochke x = " << x1 << ","<< x2 << endl;
    system("pause");
    return 0;
}