Программа для поиска минимального значения функции

#include <stdio.h>
 
int main(void)
{
    double Y,x1,x2,x3,x4,x5,Ynim,Xnim1,Xnim2,Xnim3,Xnim4,Xnim5;
    int i1,i2,i3,i4,i5;
    Ynim=1000;
    for(i1=0;i1<=10;i1++)
    {
     x1=0.2*i1-1;
     for(i2=0;i2<=10;i2++)
     {
      x2=0.2*i2-1;
      for(i3=0;i3<=10;i3++)
      {
       x3=0.2*i3-1;
       for(i4=0;i4<=10;i4++)
       {
        x4=0.2*i4-1;
        for(i5=0;i5<=10;i5++)
        {
         x5=0.2*i5-1;
        
        
        Y = 226.17 + 57.72*x1 + 17.51*x4 + 22.92*x5 - 11.43*x1*x3 + 15.22*x2*x3 + 29.83*x2*x5 + 6.86*x4*x4;
        if (Y<Ynim)
        { 
         Ynim=Y;
         Xnim1=x1;
         Xnim2=x2;
         Xnim3=x3;
         Xnim4=x4;
         Xnim5=x5;
        }
        }
       }
      }
     }
    }
    
    printf("x1 = %f; x2 = %f; x3 = %f; x4 = %f; x5 = %f; Ymin = %f\n",Xnim1,Xnim2,Xnim3,Xnim4,Xnim5,Ynim);
        
    return 0;
}

#include <iostream>
#include <iomanip>
using namespace std;
 
int main() 
{   double x2 = -22.92/29.83;
    double x3 = 57.72/11.43;
    double x4 = -17.51/13.72;
    double x1 = 15.22*x2/11.43;
    double x5 = -15.22*x3/29.83;
 
    double y = 226.17+57.72*x1+17.51*x4+22.92*x5-11.43*x1*x3+15.22*x2*x3+29.83*x2*x5+6.86*x4*x4;
 
    cout<<fixed<<setprecision(4)<<"x1 = "<<x1<< "\tx2 = "<<x2<<"\tx3 = "<<x3
        <<"\tx4 = "<<x4<<"\tx5 = "<<x5<<"\n";
 
    cout<<"Ymin = "<<y<<"\n";
        
    cin.sync(); cin.get();
    return 0;
}