Различные алгоритмы минимизации, реализованные в С++ / all / LAB3_1_1

#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream.h>

double func(double x1, double x2)
	{
	 return (2*x1*x1*x1+4*x1*x2*x2*x2-10*x1*x2+x2*x2);
	}
double df(double x1, double x2, double p1, double p2)
	{
	 return ((6*x1*x1+4*x2*x2*x2-10*x2)*p1+(12*x1*x2*x2-10*x1+2*x2)*p2);
	}

void main()
{
 clrscr();
 double x10,x20,x1,x2,p1,p2,a1,a2,b1,b2,l1,l2,m1,m2,xmin1,xmin2,eps=0.01,alfa;
 int k1=1,k2=0;
 double delta=0.1*eps;

 cout<<"‚ўҐ¤ЁвҐ X:\n";
 cin>>x10;
 cin>>x20;

 cout<<"\n‚ўҐ¤ЁвҐ ­ Їа ў«Ґ­ЁҐ P:\n";
 cin>>p1;
 cin>>p2;

 double n=fabs((func(x10,x20)-df(x10,x20,p1,p2))/df(x10,x20,p1,p2));
 if(n>1) {alfa=1;}
 else {alfa=n;}

 x1=x10+alfa*p1;
 x2=x20+alfa*p2;

	if (func(x1,x2)>func(x10,x20))
	{
	 p1=-p1;
	 p2=-p2;
	}

 x1=x10+alfa*p1;
 x2=x20+alfa*p2;

 while(func(x10,x20)>func(x1,x2))
	{
	 x10=x1;
	 x20=x2;
	 alfa=alfa*2;
	 x1=x10+alfa*p1;
	 x2=x20+alfa*p2;
	 k1++;
	}

 a1=x1-(3*alfa*p1)/2;
 a2=x2-(3*alfa*p2)/2;

 while ((a1>x1)||(a2>x2))
	{
	 if (a1>x1)
		{
		 b1=a1;
		 a1=x1;
		 b2=x2;
		}
	 else
		{
		 b1=x1;
		 b2=x2;
		};

	 if (a2>x2)
		{
		 b2=a2;
		 a2=x2;
		 b1=x1;
		}
	 else
		{
		 b1=x1;
		 b2=x2;
		}
	}

 cout<<"Ќ з «м­л© Ё­вҐаў « = "<<a1<<"     "<<b1<<"\n";
 cout<<"                     "<<a2<<"     "<<b2<<"\n";



 while (sqrt((b1-a1)*(b1-a1)+(b2-a2)*(b2-a2))>eps)
	{
	 l1=(a1+b1-delta)/2;
	 l2=(a2+b2-delta)/2;
	 m1=(a1+b1+delta)/2;
	 m2=(a2+b2+delta)/2;

	 if (func(l1,l2)<func(m1,m2))
		{
		 b1=m1;
		 b2=m2;
		}
	 else
		{
		 a1=l1;
		 a2=l2;
		}
	 k2++;
	}
 xmin1=(a1+b1)/2;
 xmin2=(a2+b2)/2;

 cout<<"ЊЁ­Ё¬г¬            = "<<xmin1<<"\n";
 cout<<"                     "<<xmin2<<"\n";
 cout<<"—Ёб«® ЁвҐа жЁ© = "<<k2;

 getch();
}