惩罚函数法

#include #include #define m 10 //f函数 float f(float x[],float p); //进退法 void mjtf(int n,float x0[],float h,float s[],float a[],float b[],float p); void mhjfgf(int n,float a[],float b[],float flag,float x[],float p); //鲍威尔法 void mbwef(int n,float x0[],float h,float flag,float a[],float b[],float x[],float p); float max(float l) { float resu; if(l<=0) resu=0; else resu=l*l;