C CHNSNG SOURCE CHAT 05/01/12 22:00:42 5004 * impose par penalisation une relation dans am * subroutine chnsng(am,nbnn,i,j,k) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) dimension am(nbnn,nbnn) * xpenl=am(i,i)*1d-3 am(i,i)=am(i,i)-2*xpenl am(i,j)=am(i,j)+xpenl am(i,k)=am(i,k)+xpenl * end