mwdfhs.c

运筹学非线性规划的c语言编写的外点惩罚函数法子程序

#include "powell.c"
void main()
{int i;
 double p[]={7,8};
 double fom,fxo,c,x[2];
 c=10;
 r0=0.1;
 fom=100;
 do 
 {fxo=powell(p,0.1,0.0001,0.0001,2,x);
  if(fabs(fom-fxo)>0.001)
    {fom=fxo;
     r0=c*r0;
     for(i=0;i<2;i++)
      *(p+i)=x[i];
    }
  else
   {printf("输出最优点及其目标函数值：\n");
    printf("x[0]=%f,x[1]=%f,ff=%f",x[0],x[1],fxo);
    return;
   }
 }while(1);
}