C-8K MODV 11-219 01.05 *;A J;*;F I=1,99;*;A Y(I);* 01.10 S ND=99;S NV=5;S SU=0;S LA=1;S N=NV; S IN=1 01.12 F I=1,IN,ND; S SU=SU+Y(I); IF (Y(I)-Y(LA)) 1.14,1.14; S LA=I 01.14 C FOR INVERSE LINE (DECREASE AT CENTER), 'IF (Y(LA)-Y(I))' IN 1.12 01.16 S V(1)=LA; S V(2)=.1; S V(4)=SU*IN/ND; S V(3)=Y(LA)-V(4); S V(5)=00 01.18 S SU=0;S Q=1; F I=1,IN,ND;D 3;S SU=SU+(SI-Y(I))^2 01.21 T !!!!, " CENTER WIDTH AMPLI BASE DISPER SUM",!! 01.26 F I=1,NV;T %7.04,V(I) 01.28 T %12.04,SU,! 01.30 F QQ=1,6;D 2;D 1.26;D 1.28 01.32 F J=1,NV;S DV(J)=FSQT(A((NV+1)*(J-2)+(NV+2))*SU/(ND-NV)) 01.34 T !!, "UNCERTAINTIES IN LAST ITERATION",!;F I=1,NV;T %7.04,DV(I) 01.36 S D=V(3)/FSQT(SU/(ND/IN)) 01.38 T !!,"SIGNAL TO NOISE RATIO",D,!!!!!!!!!! 01.40 GOTO 1.05 02.02 F J=1,NV; S C(J)=0.0; S DV(J)=.01*V(J); S DV(5)=.01 02.03 F J=1,NV; F K=1,NV; S B(NV*(J-1)+K)=0.0 02.04 F I=1,IN,ND; S Q=1; D 3; D 2.08; F J=1,NV; F K=1,NV;D 2.14 02.05 G 2.18 02.08 F J=1,NV; S V(J)=V(J)+DV(J);S Q=2; D 3;D 2.12 02.12 S V(J)=V(J)-DV(J); S C(J)=C(J)+( (G(J)-SI)/DV(J))*(SI-Y(I)) 02.14 S H=NV*(J-1)+K; D 2.16 02.16 S B(H)=B(H)+((G(J)SI)/DV(J))*((G(K)-SI)/DV(K)); S A(H)=B(H) 02.18 D 9 02.22 F J=1,NV; F K=1,NV; S V(J)=V(J)-A(NV*(J-1)+K)*C(K) 02.24 S SU=0.0;S Q=1; F I=1,IN,ND;D 3;S SU=SU+(SI-Y(I))^2 02.32 S N=N 03.01 S Z=(I-V(1))*V(2) 03.50 S D=V(3)*(1+V(5)*Z)/(1+Z*Z)+V(4) 03.70 IF (Q-1.5)3.72;S G(J)=D;G 3.9 03.72 S SI=D 03.90 S N=N 09.01 C REPLACE A SQUARE MATRIX "A(N^2)" BY ITS INVERSE 09.10 S D=1;S NK=-N; F K=1,N;DO 10 09.11 S K=N 09.12 S K=K-1; IF (-K)9.13;GOTO 9.25 09.13 S I=L(K); IF (I-K) 9.17,9.17,9.14 09.14 S JQ=N*(K-1); S JR=N*(I-1); F J=1,N ; S JK=JQ+J; DO 9.16 09.15 GOTO 9.17 09.16 S HOLD=A(JK); S JI=JR+J; S A(JK)=-A(JI); S A(JI)=HOLD 09.17 S J=M(K); IF (J-K) 9.12,9.12,9.18 09.18 S KI=K-N; F I=1,N; S KI=KI+N; S HOLD=A(KI); DO 9.20 09.19 GOTO 9.12 09.20 S JI=KI-K+J; S A(KI)=-A(JI); S A(JI)=HOLD 09.25 S N=N 10.10 S NK=NK+N; S L(K)=K; S M(K)=K; S KK=NK+K; S BIGA=A(KK) 10.11 F J=K,N; S IZ=N*(J-1); F I=K,N;S IJ=IZ+I; DO 10.13 10.12 GOTO 10.16 10.13 IF (FABS(BIGA)-FABS(A(IJ))) 10.14;C 10.14 S BIGA=A(IJ); S L(K)=I; S M(K)=J 10.16 S J=L(K); IF (J-K) 10.20,10.20;S KI=K-N; F I=1,N;DO 10.18 10.17 GOTO 10.20 10.18 S KI=KI+N; S HOLD=-A(KI); S JI=KI-K+J; DO 10.19 10.19 S A(KI)=A(JI); S A(JI)=HOLD 10.20 S I=M(K); IF (I-K) 10.24,10.24; S JP=N*(I-1); F J=1,N; DO 10.22 10.21 GOTO 10.24 10.22 S JK=NK+J; S JI=JP+J; S HOLD=-A(JK); S A(JK)=A(JI); S A(JI)=HOLD 10.24 IF (BIGA) 10.26,10.25,10.26 10.25 S D=0.0; R 10.26 F I=1,N; IF (-FABS(I-K)) 10.28; C 10.27 GOTO 10.29 10.28 S IK=NK+I; S A(IK)=A(IK)/(-BIGA) 10.29 F I=1,N; S IK=NK+I; S IJ=I-N; DO 10.31 10.30 GOTO 10.35 10.31 F J=1,N; S IJ=IJ+N; IF (-FABS(I-K)) 10.32; C 10.32 IF (-FABS(J-K)) 10.33; C 10.33 S KJ=IJ-I+K; S A(IJ)=A(IK)*A(KJ)+A(IJ) 10.35 S KJ=K-N; F J=1,N; S KJ=KJ+N; IF (-FABS(J-K)) 10.37; C 10.36 GOTO 10.38 10.37 S A(KJ)=A(KJ)/BIGA 10.38 S D=D*BIGA; S A(KK)=1.0/BIGA *