mldfhs1.c

来自「运筹学非线性规划C语言编写的内点惩罚函数法子程序」· C语言 代码 · 共 23 行

#include "dfpopt.c"
main()
{int i;
 double p[]={3,4};
 double fom,fxo,x[2];
 double c=0.1;
 double  r0=120;
 fom=100;
 do
  { fxo=dfpopt(p,0.2,0.0001,0.000001,2);
    if(fabs(fom-fxo)>0.001)
    {fom=fxo;
     r0=c*r0;
     for(i=0;i<2;i++)
     *(p+i)=x[i];
    }
   else
    {printf("\nx[0]=%f,x[1]=%f,ff=%f",x[0],x[1],fxo);
     return; }
    }while(1);
    getch();
}