CALL AXES(XL,-2.0,5HSL0PE,-5,2.0,90.,1.0,-4.0,2.0,-l) 

 CALL PL0T(0.0,0.0,3) 

 859 AVSLP=0.0 

 KKNT=0 



SLWID=SLWID2-SLWID1 

 ENWID=ENWI02-ENWID1 

 INWID1=IFIX(ENWID1) 

 LSWID1=IFIX(SLWID1) 

 INWID2=IFIX(ENWID2) 

 LSWID2=IFIX(SLWID2) 



DO 877 K=(NTWICE+l)+(NFIT/2),XT-((NTWICE+l)+(NFIT/2)) 

 KKNT=KKNT+1 

 C=(K-1)*.02 

 B=(TZ(K,ll))/2. 



AVSLP=AVSLP+(TZ(K,11)/(JCT-(4*NUMFIL+(NFIT+1)))) 

 IF(IPLOT.NE.l) GO TO 876 

 IF(KKNT.EQ.l) CALL PLOT (C,B,3) 

 CALL PL0T(C,B,2) 



876 IF(KKNT.LE.LSWID2) GO TO 877 

 IF(K.GE.(JCT-(NTWICE+NFIT/2+l+LSWID2))) GO TO 877 

 PREAV=0.0 



POSTAV=0.0 

 C DO 8761 KNB=LSWID1,LSWID2 

 C PREAV=PREAV+(TZ(K-KNB,11)/SLWID) 

 C8761 P0STAV=P0STAV+(TZ(K+KNB,11)/SLWID) 

 C IF(ABS(POSTAV-PREAV).LT.SLWIN) GO TO 877 

 C CALL SYMB0L(C,-2., .2,16,0.0,-1) 

 C CALL PL0T(C,B,3) 



877 CONTINUE 



C PLOT ARRAY OF RMS ENERGY 



IF(IPLOT.EQ.l) CALL NEWPEN(l) 



KKNT=0 



AVENER=0.0 



DO 878 L=(NTWICE+l)+(NFIT/2),JCT-((NTWICE+l)+(NFIT/2)) 

 C CALCULATE BAND LIMITED RMS-FIRST SQUARE FUNCTION AND INTEGRATE 

 C INTEGRAL((A*F**B)**2)=(A**2/(2*B+1)*F**(2*B+1) 



XINTL0=(TZ(L,12)**2/(2.*TZ(L,11)+1.))*FR1**(2.*TZ(L,11)+1.) 



XINTHI=(TZ(L,12)**2/(2.*TZ(L,11)+1.))*FR2**(2.*TZ(L,11)+1.) 

 C SINCE SPECTRUM IS SYMMETRIC ABOUT ZERO, CAN EVALUATE INTEGRAL 

 C BY MULTIPLYING EVALUATED RESULT BY TWO 



ENERGY(L)=2.*(XINTHI-XINTL0) 

 C TO CALCULATE MEAN SQUARE, DIVIDE BY WIDTH 



ENERGY(L)=ENERGY(L)/(2.*(FR2-FR1)) 

 C NOW TAKE SQUARE ROOT TO DETERMINE RMS 



ENERGY(L)=SQRT(ENERGY(L)) 



878 AVENER=AVENER+(ENERGY(L)/(JCT-(4*NUMFIL+(NFIT+1)))) 

 ENSCAL=4.*AVENER 



IF(IPLOT.NE.l) GO TO 8782 



CALL AXES(0.0,0.0,10HRMS ENERGY,10,2. ,90. ,0.5,0.0.AVENER/2. , 

 *-l) 

 CALL PL0T(0.0,0.0,3) 

 8782 DO 8781 L = (NTWICE+l)+(NFIT/2) ,XT-((NTWICE+l)+(NFIT/2) ) 



154 



