5 * Revision 1.1.1.1 1996/04/01 15:02:20 mclareni
10 SUBROUTINE MINSQ (M,N,F,X,E,IPRINT,NFUN ,NW,W,COV,XSTEP)
11 DIMENSION F(M),X(N),E(N),W(NW)
20 CALL FCN1(M,N,F,X,IFLAG,NW,W,KEND,KEND2N)
32 CALL FCN1(M,N,F,X,IFLAG,NW,W,KEND,KEND2N)
48 7 FORMAT(5X,'E(',I3,') UNREASONABLY SMALL')
59 W(II)=W(II)+W(KK)*W(K)
76 W(IIP)=W(IIP)+W(KK)*W(II)
78 IF (JL-ILESS) 18,18,19
82 W(IIP)=W(IIP)+W(KK)*W(JJ)
83 W(JJP)=W(JJP)+W(KK)*W(II)
100 15 GO TO (27,24),IINV
114 IPP=IPRINT*(IPRINT-1)
120 28 WRITE(6,30)ITC,MC,FF
121 30 FORMAT(//5X,'ITERATION',I4,I9,' CALLS OF FCN ',5X,'F=',E24.14)
122 WRITE(6,31)(X(I),I=1,N)
123 31 FORMAT(5X,'VARIABLES',/(5E24.14))
124 WRITE(6,32)(F(I),I=1,M)
125 32 FORMAT(5X,'FUNCTIONS',/(5E24.14))
128 29 GO TO (34,35),ICONT
129 35 IF (CHANGE-1.) 10,10,36
133 38 FORMAT(//5X,' FINAL VALUES OF FUNCTIONS AND VARIABLES')
137 CALL FCN1(M,N,F,X,IFLAG,NW,W,KEND,KEND2N)
148 90 W(IWC)=W(IWC)+W(ICOV)*W(JCOV)
151 1000 FORMAT(///40X,'VARIANCE-COVARIANCE MATRIX'///)
156 92 WRITE(6,1001)(W(J),J=IA,IB)
157 1001 FORMAT(10X,5G20.6//)
177 W(IIP)=W(IIP)+W(K)*W(II)
183 W(IIP)=W(IIP)+W(K)*W(JJ)
184 W(JJP)=W(JJP)+W(K)*W(II)
187 43 IF (DM-ABS (W(II)*W(IIP))) 45,41,41
188 45 DM=ABS (W(II)*W(IIP))
203 IF (ABS (E(I)*CHANGE)-ABS (W(I))) 48,48,46
204 48 CHANGE=ABS (W(I)/E(I))
218 IF (CHANGE-1.) 50,50,51
220 51 CALL VD01A (IT,XC,FC,20,ACC,0.1,XSTEP)
221 GO TO (52,53,53,53),IT
223 IF (MC-MAXFUN) 54,54,55
225 56 FORMAT(5X,I6,' CALLS OF FCN')
233 CALL FCN1(M,N,F,X,IFLAG,NW,W,KEND,KEND2N)
249 59 IF (FC-FSEC) 64,51,51
253 74 IF (FC-FMIN) 65,51,66
324 IF (FMIN-FF) 82,83,83
328 CHANGE=ABS (XC)*CHANGE
330 SUM=1./SQRT (SUM+DM*XL)
339 W(K)=SUM*(W(K)+XL*F(I))