机械优化设计惩罚函数法多变量约束优化.doc

C=PROGRAM COMPLEC=DIMENSION X(25),GX(50),XCOM(1250)COMMON /ONE/ ITE,IXE,ILI,NPE,NFX,NGRREAD(*,*) N,KG,KWRITE(*,10001) N,KG,K10001FORMAT(25X,= PRIMARY DATA =/5X, 1 N=,I4,5X,KG=,I4,5X,K=,I4)CALL MAISUB(N,K,KG,X,GX,XCOM)STOPENDC=SUBROUTINE MAISUB(N,K,KG,X,GX,XCOM)C=DIMENSION X(N),GX(KG),XCOM(N,K),FXK(50),XR(25)DIMENSION