内点惩罚函数法子程序6页.doc

#include stdio.h#include stdlib.h#include math.hconst int kkg=3;double r0;double f(double x)double ff; ff=pow(x0-8),2)+pow(x1-8),2); return(ff);/*约束条件子程序*/void strain(double x,double g)g0=x0-1; g1=x1-1; g2=11-x0-x1;/*惩罚函数子程序*/double objf(double p)int i; double ff,sg,*g;g=(double *)malloc(kkg*sizeof(double);sg=0;strain(p,g);for(i=0;i0) sg=sg+r0/(*(g+i); elsesg=