机械优化设计惩罚函数内点法(共6页).doc

精选优质文档-倾情为你奉上#include #include #define m 12double f(double x,double r);void jintuifa(double abmm,int n,double x0,double h,int ij,double a,double b,double r0);void hongjinfa(int n,double a,double b,double flag,double x,double r0);void baoweifa(int n,double x0,double h,double flag,double a,double b,double x,double r0);double fahansu(double tt) double ty; if(tt0) ty=-tt;else ty=0; return ty;double yuanhansu(dou