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

void newton (double x01, double x02)
{
  double p1=0, p2=0, x1, x2, h11, h12, h21, h22, alfa_r=0;
  int m=2000, k=1;
  x1=x01;
  x2=x02;
  agrad(x1,x2,p1,p2);
  H (x1, x2, p1, p2, h11,h12,h21,h22);
  p1=h11*p1+h12*p2;
  p2=h21*p1+h22*p2;
  while ((sqrt(p1*p1+p2*p2)>0.01)&&(k<=m))
  {
	 p1=p1/(sqrt(pow(p1,2)+pow(p2,2)));
	 p2=p2/(sqrt(pow(p1,2)+pow(p2,2)));
	 lab3(x1,x2,p1,p2,alfa_r);
	 x1=x1+alfa_r*p1;
	 x2=x2+alfa_r*p2;
	 alfa_r=0;
	 agrad(x1,x2,p1,p2);
	 H(x1, x2, p1, p2, h11,h12,h21,h22);
	 p1=h11*p1+h12*p2;
	 p2=h21*p1+h22*p2;
	 k++;
  }
//  cout<<"sqrt(p1*p1+p2*p2)="<<fabs(p1*p1+p2*p2)<<endl;
//  cout<<"x1="<<x1<<endl;
//  cout<<"x2="<<x2<<endl;
//  cout<<"Y="<<Y(0, x1, x2, p1, p2)<<endl;
  cout<<"k="<<k-1<<endl;
}