SUBROUTINE MDL0(K,L,MATRIX,A,HEADER,YRONE)
C
      IMPLICIT DOUBLE PRECISION (A-G,O-Z)
      INTEGER YRONE,BANDS
      DOUBLE PRECISION MATRIX,MAT,N, LOWS,HIS,LOWFX,
     1 HIFX,LOWF,HIF,MLS,LOWAFX,LOWAF,LOWMLS
      CHARACTER HEADER*80, FORMT*30
      DIMENSION MATRIX(21,21),A(20)
      COMMON /MDL0CM/ N(20),R(20),C(20),T(20),F(20)
     1 , FX(20),S(20),E(20,20),VARF(20),VARFX(20),VARS(20)
     1 , SDF(20),SDFX(20),SDS(20),IMATT(20)
     1 , CVFFX(20),CVFXS(20),CVFXPS(20),CVFS(20),CVSSP(20)
     1 , CVFPS(20),ISAVE(20),IMAT(20,20),MAT(20,20),GAMMA(20)
     1 , CHISQ(20,20),CRFFX(20),CRFXS(20),CRFXPS(20),CRFS(20)
     1 , CRSSP(20),CRFPS(20),RHO(20),SFPROD(20),JYEAR(20),IN(20)
     1 , SML(20),FML(20), TMSPAC(479)
C
      DIMENSION THCHI(100)
      DATA THCHI/3.84D0,5.99D0,7.81D0,9.49D0,11.10D0,12.60D0,14.10D0,
     1     15.50D0,16.90D0,18.30D0,19.70D0,21.00D0,22.40D0,23.70D0,
     2     25.00D0,26.30D0,27.60D0,28.90D0,30.10D0,31.40D0,32.70D0,
     3     33.90D0,35.20D0,36.40D0,37.70D0,38.90D0,40.10D0,41.30D0,
     4     42.60D0,43.80D0,44.98D0,46.19D0,47.40D0,48.60D0,49.80D0,
     5     50.99D0,52.19D0,53.38D0,54.57D0,55.76D0,56.94D0,58.12D0,
     6     59.30D0,60.48D0,61.65D0,62.83D0,63.99D0,65.16D0,66.34D0,
     7     67.50D0,68.67D0,69.83D0,70.99D0,72.15D0,73.30D0,74.46D0,
     8     75.62D0,76.78D0,77.93D0,79.08D0,80.23D0,81.38D0,82.53D0,
     9     83.67D0,84.81D0,85.96D0,87.10D0,88.25D0,89.38D0,90.53D0,
     $     91.66D0,92.80D0,93.94D0,95.08D0,96.21D0,97.34D0,98.48D0,
     $ 99.61D0,100.74D0,101.88D0,103.00D0,104.14D0,105.26D0,106.39D0,
     $107.52D0,108.64D0,109.77D0,110.89D0,112.02D0,113.14D0,114.26D0,
     $115.38D0,116.51D0,117.63D0,118.75D0,119.87D0,120.98D0,122.10D0,
     $123.22D0,124.34D0/
C
C      INITIALIZE VARIABLES
      MDIFF=L-K
      MDIFFP=MDIFF+1
      LL=L
      KPLUS=K+1
      KN=K+1
      KK=K-1
      AVEF=0.0
      AVEFX=0.0
      AVES=0.0
      CHISUM=0.0
      SUMS=0.0
      SUMVAR=0.0
      DF=0.0
      DO 5 J=1,L
      IN(J)=A(J)
      N(J)=A(J)
      R(J)=0.0
      C(J)=0.0
      T(J)=0.0
      S(J)=0.0
      SML(J)=0.0
      FML(J)=0.0
      VARS(J)=0.0
      VARF(J)=0.0
      VARFX(J)=0.0
      SDS(J)=0.0
      FX(J)=0.0
      ISAVE(J)=0
      GAMMA(J)=0.0
      CVFPS(J)=0.0
      CRFFX(J)=0.0
      CRFXS(J)=0.0
      CRFXPS(J)=0.0
      CRFS(J)=0.0
      CRSSP(J)=0.0
      CRFPS(J)=0.0
      CVFXS(J)=0.0
      CVFXPS(J)=0.0
      CVFS(J)=0.0
      CVSSP(J)=0.0
      SDF(J)=0.0
      SDFX(J)=0.0
      RHO(J)=0.0
      F(J)=0.0
      CVFFX(J)=0.0
      DO 5 I=1,K
      IMAT(I,J)=MATRIX(I,J)
      E(I,J)=0.0
      CHISQ(I,J)=0.0
      MAT(I,J)=MATRIX(I,J)
    5 CONTINUE
C
      CALL MOD0W1
      WRITE (6,80) HEADER
   80 FORMAT (' ',20X,A///)
      WRITE (6,85)
   85 FORMAT (' ','BANDING AND RECOVERY INPUT DATA'/)
      WRITE (6,90)
   90 FORMAT (/' ','YEAR NUMBER',35X,'RECOVERY MATRIX')
      WRITE (6,95)
   95 FORMAT (' ',5X,'BANDED')
      WRITE (6,100)
  100 FORMAT ('+','____ ______',35X,'_______________'/)
      IYEAR=YRONE-1
      DO 110 I=1,K
      IYEAR=IYEAR+1
      IS=(I-1)*6+1
      WRITE(FORMT,105) IS
      WRITE (6,FORMT) IYEAR,IN(I),(IMAT(I,J),J=I,L)
  105 FORMAT ('(1X,I4,1X,I5,',I2,'X,20I6)')
  110 CONTINUE
C     CALCULATE INTERMEDIATE STATISTICS
      DO 115 I=1,K
      DO 115 J=1,L
      R(I)=R(I)+MATRIX(I,J)
  115 CONTINUE
      SUMA=0.0
      SUMR=0.0
      DO 120 I=1,K
      SUMA=SUMA+A(I)
      SUMR=SUMR+R(I)
  120 CONTINUE
C
      DO 125 J=1,L
      DO 125 I=1,K
      C(J)=C(J)+MATRIX(I,J)
  125 CONTINUE
C
      T(1)=R(1)
      DO 130 I=2,K
      T(I)=T(I-1)-C(I-1)+R(I)
  130 CONTINUE
C
      DO 145 I=1,K
      RHO(I)=R(I)/A(I)
      FX(I)=MATRIX(I,I)/A(I)
      IF (I.EQ.1) GO TO 135
      IF (MDIFF.EQ.0.AND.I.EQ.K) GO TO 135
      IF (R(I)-MATRIX(I,I).EQ.0) GO TO 140
      F(I)=((R(I)-MATRIX(I,I))/A(I))*((C(I)-MATRIX(I,I))/(T(I)-R(I)-C(I)
     1+MATRIX(I,I)+1.0))
      FML(I)=((R(I)-MATRIX(I,I))/A(I))*((C(I)-MATRIX(I,I))/(T(I)-R(I)-C(
     1I)+MATRIX(I,I)))
  135 IF (I.EQ.K) GO TO 140
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 140
      S(I)=((R(I)-MATRIX(I,I))/A(I))*((A(I+1)+1.0)/(R(I+1)-MATRIX(I+1,I+
     11)+1.0))*(1.0-((C(I+1)-MATRIX(I+1,I+1))/(T(I+1)-R(I+1))))
      SML(I)=((R(I)-MATRIX(I,I))/A(I))*((A(I+1))/(R(I+1)-MATRIX(I+1,I+
     11)))*(1.0-((C(I+1)-MATRIX(I+1,I+1))/(T(I+1)-R(I+1))))
  140 CONTINUE
      VARFX(I)=(FX(I)*(1.0-FX(I)))/A(I)
      IF (VARFX(I).GT.0.0) SDFX(I)=DSQRT(VARFX(I))
      IF (I.EQ.1) GO TO 145
      IF (MDIFF.EQ.0.AND.I.EQ.K) GO TO 145
      VARF(I)=F(I)**2*((1.0/(R(I)-MATRIX(I,I)))-(1.0/A(I))+(1.0/(T(I)-R(
     1I)-C(I)+MATRIX(I,I)))+(1.0/(C(I)-MATRIX(I,I))))
      IF (VARF(I).GT.0.0) SDF(I)=DSQRT(VARF(I))
      IF (MDIFF.EQ.0) SDF(K)=99999.
  145 CONTINUE
C
      DO 150 I=KN,L
      GAMMA(I)=(C(I)*F(K))/C(K)
  150 CONTINUE
      DO 155 I=1,KK
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 155
      VARS(I)=S(I)**2*((1.0/(R(I)-MATRIX(I,I)))-(1.0/A(I))+(1.0/(R(I+1)-
     1MATRIX(I+1,I+1)))-(1.0/A(I+1))+(1.0/(T(I+1)-R(I+1)-C(I+1)+MATRIX(I
     2+1,I+1)))-(1.0/(T(I+1)-R(I+1))))
      IF (VARS(I).GT.0.0) SDS(I)=DSQRT(VARS(I))
      IF (VARS(I).LT.0) SDS(I)=99999.
      IF (MDIFF.EQ.0) SDS(K-1)=99999.
      SUMVAR=SUMVAR+VARS(I)
      SUMS=SUMS+S(I)
  155 CONTINUE
      IF (MDIFF.EQ.0) AVES=SUMS/(K-2)
      IF (MDIFF.GT.0) AVES=SUMS/(K-1)
C
      IF (K.LT.4) GO TO 190
      DO 185 I=1,K
      IF (I.EQ.1) GO TO 160
      IF (MDIFF.EQ.0.AND.I.EQ.K) GO TO 160
      CVFFX(I)=-1.0*((FX(I)*F(I))/A(I))
  160 IF (I.EQ.K) GO TO 165
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 165
      CVFXS(I)=-1.0*((FX(I)*S(I))/A(I))
  165 IF (I.EQ.K) GO TO 170
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 170
      CVFXPS(I)=(FX(I+1)*S(I))/A(I+1)
  170 IF (I.EQ.1.OR.I.EQ.K) GO TO 175
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 175
      CVFS(I)=F(I)*S(I)*((1.0/(R(I)-MATRIX(I,I)))-(1.0/A(I)))
  175 IF (I.EQ.K) GO TO 180
      IF (MDIFF.EQ.0.AND.I.EQ.KK) GO TO 180
      CVFPS(I)=-1.0*(S(I)*F(I+1))*((1.0/(R(I+1)-MATRIX(I+1,I+1)))-(1.0/A
     1(I+1)))
  180 IF (I.GT.(K-2)) GO TO 185
      IF (MDIFF.EQ.0.AND.I.EQ.(K-2)) GO TO 185
      CVSSP(I)=-1.0*(S(I)*S(I+1))*((1.0/(R(I+1)-MATRIX(I+1,I+1)))-(1.0/A
     1(I+1)))
  185 CONTINUE
  190 CONTINUE
      WRITE (6,195)
  195 FORMAT (///1X,28X,'RECOVERY RATE    F(I)  (%)',33X,'FIRST-YEAR RE'
     1,'COVERY RATE    F*(I)  (%)')
      WRITE (6,200)
  200 FORMAT ('+',15X,'_______________________________________________',
     1'____________',9X,'_____________________________________________',
     2'__')
      WRITE (6,205)
  205 FORMAT (1X,6X,'I',8X,'ESTIMATE',14X,'STANDARD',14X,'95% CONFIDENC'
     1,'E',10X,'ESTIMATE',9X,'STANDARD',6X,'95% CONFIDENCE')
      WRITE (6,210)
  210 FORMAT (1X,37X,'  ERROR  ',16X,'INTERVAL',30X,'  ERROR  ',8X,'INT'
     1,'ERVAL')
      WRITE (6,215)
  215 FORMAT ('+',15X,'________',14X,'_________',13X,'______________',10
     1X,'________',9X,'_________',5X,'______________'/)
      IYEAR=YRONE-1
      SUMONE=0.0
      SUMTWO=0.0
      DO 225 I=1,K
      IYEAR=IYEAR+1
      IF (I.GT.1) SUMONE=SUMONE+VARF(I)
      SUMTWO=SUMTWO+VARFX(I)
      PF=F(I)*100.0
      PSDF=SDF(I)*100.0
      IF (I.EQ.1) PF=99999.
      IF (I.EQ.1) PSDF=99999.
      IF (I.EQ.K.AND.MDIFF.EQ.0) PF=99999.
      LOWF=PF-1.96*PSDF
      HIF=PF+1.96*PSDF
C
      PFX=FX(I)*100.0
      PSDFX=SDFX(I)*100.0
      LOWFX=PFX-1.96*PSDFX
      HIFX=PFX+1.96*PSDFX
C
      IF(I.NE.K.OR.MDIFF.NE.0) GO TO 217
      WRITE (6,216) IYEAR,PSDF,LOWF,HIF,PFX,PSDFX,LOWFX,HIFX
  216 FORMAT(' ',4X,I4,7X,'*******',16X,F7.3,11X,F7.3,' - ',F7.3,10X,
     1F7.3,10X,F7.3,4X,F7.3,' - ',F7.3/)
C
      GO TO 225
  217 IF(I.NE.1) GO TO 219
      WRITE (6,218) IYEAR,PFX,PSDFX,LOWFX,HIFX
  218 FORMAT (' ',4X,I4,7X,'*******',16X,'*******',11X,'*******',
     1' - ','*******',10X,F7.3,10X,F7.3,4X,F7.3,' - ',F7.3/)
C
      GO TO 225
  219 WRITE (6,220) IYEAR,PF,PSDF,LOWF,HIF,PFX,PSDFX,LOWFX,HIFX
  220 FORMAT (' ',4X,I4,7X,F7.3,16X,F7.3,11X,F7.3,' - ',F7.3,10X,F7.3,10
     1X,F7.3,4X,F7.3,' - ',F7.3/)
  225 CONTINUE
C
      WRITE (6,240)
      SUMF=0.0
      SUMFX=0.0
      AVEF=0.0
      AVEFX=0.0
      DO 230 I=1,K
      SUMF=SUMF+F(I)
      SUMFX=SUMFX+FX(I)
  230 CONTINUE
      IF (MDIFF.EQ.0.AND.(K-2).GT.0) AVEF=SUMF/(K-2)
      IF (MDIFF.EQ.0.AND.(K-2).EQ.0) AVEF=SUMF
      IF (MDIFF.GT.0) AVEF=SUMF/(K-1)
      AVEFX=SUMFX/K
      PAVEF=AVEF*100.0
      PAVEFX=AVEFX*100.0
      VARAFX=SUMTWO/(K*K)
      IF (MDIFF.EQ.0.AND.(K-2).GT.0) VARAF=SUMONE/((K-2)**2)
      IF (MDIFF.EQ.0.AND.(K-2).EQ.0) VARAF=VARF(2)
      IF (MDIFF.GT.0.AND.(K-2).GT.0) VARAF=SUMONE/((K-1)**2)
      IF (MDIFF.GT.0.AND.(K-2).EQ.0) VARAF=VARF(2)
      IF (VARAFX.GT.0.0) SDAFX=DSQRT(VARAFX)
      IF (VARAF.GT.0.0) SDAF=DSQRT(VARAF)
      PSDAF=SDAF*100.0
      PSDAFX=SDAFX*100.0
      LOWAF=PAVEF-1.96*PSDAF
      LOWAFX=PAVEFX-1.96*PSDAFX
      HIAF=PAVEF+1.96*PSDAF
      HIAFX=PAVEFX+1.96*PSDAFX
      WRITE (6,235) PAVEF,PSDAF,LOWAF,HIAF,PAVEFX,PSDAFX,LOWAFX,HIAFX
  235 FORMAT (/' ',4X,'MEAN',7X,F7.3,16X,F7.3,11X,F7.3,' - ',F7.3,10X,F7
     1.3,10X,F7.3,4X,F7.3,' - ',F7.3/)
C
      WRITE (6,240)
  240 FORMAT (' ',130('_'))
      WRITE (6,245)
  245 FORMAT (1X,'(QUANTITIES SHOWN AS ******* ARE NOT ESTIMABLE UNDER'
     1 ,' MODEL 0)'//)
      IF (K.GT.3) CALL FANDFP(K,L,R,C,T,MATRIX,A,YRONE)
      WRITE (6,250) HEADER
  250 FORMAT (/' ',20X,A/)
      WRITE (6,255)
  255 FORMAT (' ',34X,'SURVIVAL S(I) (%)')
      WRITE (6,290)
      WRITE (6,260)
  260 FORMAT (1X,3X,'YEAR',8X,'SURVIVAL',7X,'STANDARD',6X,
     1 'COEFFICIENTS',4X,'95% CONFIDENCE')
      WRITE (6,265)
  265 FORMAT (' ',30X,'  ERROR  ',5X,'OF VARIATION',7X,'INTERVAL')
      WRITE (6,270)
  270 FORMAT ('+',15X,'________',7X,'_________',5X,'____________',4X,'__
     1______________')
      IYEAR=YRONE-1
C
      AVES=0.0
      DO 280 I=1,KK
      IF (I.EQ.(K-1).AND.MDIFF.EQ.0) GO TO 280
      IYEAR=IYEAR+1
      CV=(SDS(I)/S(I))*100.0
      PS=S(I)*100.0
      PSDS=SDS(I)*100.0
      LOWS=PS-1.96*PSDS
      HIS=PS+1.96*PSDS
      AVES=AVES+PS
C
      WRITE (6,275) IYEAR,PS,PSDS,CV,LOWS,HIS
  275 FORMAT (' ',3X,I4,3(9X,F6.2),7X,F6.2,' - ',F6.2/)
  280 CONTINUE
C
      IF (MDIFF.EQ.0) AVES=AVES/(K-2)
      IF (MDIFF.GT.0) AVES=AVES/(K-1)
      PAVES=AVES
      SUMCOV=0.0
      DO 285 I=1,KK
      SUMCOV=SUMCOV+CVSSP(I)
  285 CONTINUE
      IF (MDIFF.EQ.0.AND.(K-3).GT.0) VARAVS=(SUMVAR+2.*SUMCOV)/((K-2)**2
     1)
      IF (MDIFF.EQ.0.AND.(K-3).EQ.0) VARAVS=VARS(1)
      IF (MDIFF.GT.0.AND.(K-2).GT.0) VARAVS=(SUMVAR+2.*SUMCOV)/((K-1)**2
     1)
      IF (MDIFF.GT.0.AND.(K-2).EQ.0) VARAVS=VARS(1)
      IF (VARAVS.GT.0.0) SDAVS=DSQRT(VARAVS)
      PSDAVS=SDAVS*100.0
      CVAVES=(PSDAVS/PAVES)*100.0
      LOWS=PAVES-1.96*PSDAVS
      HIS=PAVES+1.96*PSDAVS
      WRITE (6,290)
  290 FORMAT (' ','_____________________________________________________
     1_______________________')
      WRITE (6,295) PAVES,PSDAVS,CVAVES,LOWS,HIS
  295 FORMAT (/' ',3X,'MEAN',9X,F6.2,9X,F6.2,9X,F6.2,7X,F6.2,' - ',F6.2/
     1)
      WRITE (6,290)
C
      AVES=AVES/100.0
      MLS=-1.0/DLOG(AVES)
      VARMLS=(1.0/AVES**2)*(1.0/(DLOG(AVES)**4))*VARAVS
      SDMLS=(MLS**2*SDAVS)/AVES
      IF (LOWS/100..GT.0.0) LOWMLS=-1.0/DLOG(LOWS/100.)
      IF (HIS/100..LT.1.0) HIMLS=-1.0/DLOG(HIS/100.)
      WRITE (6,300) MLS
  300 FORMAT (//' ',2X,'MEAN LIFE SPAN AS AN ADULT = ',F7.2)
      WRITE (6,305) SDMLS
  305 FORMAT (' ',2X,'STANDARD ERROR OF THE MEAN LIFE SPAN = ',F7.2)
      IF(LOWS/100..GT.0.AND.HIS/100..LT.1.0) GO TO 311
      IF (LOWS/100..GT.0.AND.HIS/100..GT.1.0) GO TO 307
      IF (LOWS/100..LT.0.AND.HIS/100..LT.1.0) GO TO 309
      WRITE (6,306)
306   FORMAT(' ',2X,'95% CONFIDENCE INTERVAL OF LIFE ',
     1'SPAN ******* - *******'//)
      GO TO 313
307   WRITE (6,308) LOWMLS
308   FORMAT(' ',2X,'95% CONFIDENCE INTERVAL OF LIFE ',
     1'SPAN ',F7.2,' - *******'//)
      GO TO 313
309   WRITE (6,310) HIMLS
310    FORMAT(' ',2X,'95% CONFIDENCE INTERVAL OF LIFE ',
     1'SPAN ******* - ',F7.2//)
      GO TO 313
311   WRITE (6,312) LOWMLS,HIMLS
312   FORMAT (' ',2X,'95% CONFIDENCE INTERVAL OF LIFE SPAN ',F7.2,' - ',
     1F7.2//)
313   CONTINUE
      DO 315 I=1,K
      E(I,I)=A(I)*FX(I)
315   CONTINUE
      SFPROD(1)=(R(K-1)-MATRIX(K-1,K-1))/N(K-1)
      DO 325 J=2,K
      II=J-1
      DO 325 I=1,II
      PROD=1.0
      DO 320 KM=I,II
      IF (SML(KM).NE.0.0) PROD=PROD*SML(KM)
  320 CONTINUE
      IF (J.LT.L) E(I,J)=A(I)*PROD*FML(J)
      IF (J.EQ.L) E(I,J)=A(I)*PROD*SFPROD(1)
  325 CONTINUE
      IF (MDIFF.EQ.0) GO TO 336
      DO 335 I=1,MDIFF
      SUMC=0.0
      DO 330 J=1,MDIFF
      SUMC=SUMC+C(K+J)
  330 CONTINUE
      SFPROD(I+1)=(R(K)-MATRIX(K,K))/N(K)*(C(K+I)/SUMC)
  335 CONTINUE
  336 IF (L-K) 355,355,340
  340 SML(K)=1.0
      DO 350 I=1,K
      DO 350 J=KPLUS,L
      PROD=1.0
      DO 345 KM=I,K
      IF (SML(KM).NE.0.0) PROD=PROD*SML(KM)
  345 CONTINUE
      E(I,J)=A(I)*PROD*SFPROD(J-K+1)
  350 CONTINUE
C
      SML(K)=0.0
  355 CONTINUE
C
      IF (MDIFF.EQ.0) DF=(K*(K+1))/2-(3*K-3)
      IF (MDIFF.GT.0) DF=((K*(K+1))/2)+(K*MDIFF)-(3*K-2+MDIFF)
      DO 375 I=1,K
      DO 360 M=1,LL
      IF (I.EQ.K.AND.M.EQ.LL.AND.MDIFF.EQ.0) GO TO 375
      IF (E(I,LL-M+1).GT.2.0) GO TO 365
      E(I,LL)=E(I,LL)+E(I,LL-M)
      E(I,LL-M)=0.0
      IF (E(I,LL).GT.2.0) GO TO 370
  360 CONTINUE
  365 ISAVE(I)=M-1
      GO TO 375
  370 ISAVE(I)=M
C     ISAVE(I) = NUMBER OF DEGREES OF FREEDOM LOST IN ROW I
  375 CONTINUE
      DO 385 I=1,K
      MM=ISAVE(I)
      DO 380 M=1,MM
      IF (ISAVE(I).NE.0) MAT(I,L)=MAT(I,L)+MAT(I,L-M)
      IF (ISAVE(I).NE.0) MAT(I,L-M)=0.0
  380 CONTINUE
  385 CONTINUE
      WRITE (6,390)
  390 FORMAT (///)
      WRITE (6,90)
      WRITE (6,95)
      WRITE (6,100)
      IYEAR=YRONE-1
      DO 415 I=1,K
      BANDS=A(I)
      IYEAR=IYEAR+1
      DO 405 J=1,L
      IMATT(J)=MAT(I,J)
  405 CONTINUE
      IS=(I-1)*6+1
      WRITE(FORMT,105) IS
      WRITE (6,FORMT) IYEAR,BANDS,(IMATT(J),J=I,L)
  415 CONTINUE
      WRITE (6,420)
  420 FORMAT (//' ','MATRIX OF EXPECTED VALUES UNDER MODEL 0'/)
      DO 430 I=1,K
      IS=(I-1)*6+13
      WRITE(FORMT,425) IS
      WRITE (6,FORMT) (E(I,J),J=I,L)
  425 FORMAT ('(',I2,'X,20F6.1)')
  430 CONTINUE
C
      CHISUM=0.0
      DO 435 I=1,K
      DO 435 J=1,L
      IF (E(I,J).GT.0.0 .AND.E(I,J).LT.N(I)) THEN
         CHISQ(I,J)=((MAT(I,J)-E(I,J))**2)/E(I,J)
      ELSE
         CHISQ(I,J)=0.
      ENDIF
      CHISUM=CHISUM+CHISQ(I,J)
  435 CONTINUE
      WRITE (6,440)
  440 FORMAT (//' ','MATRIX OF CHI-SQUARE VALUES UNDER MODEL 0'/)
C
      DO 450 I=1,K
      IS=(I-1)*6+13
      WRITE(FORMT,445) IS
  445 FORMAT ('(',I2,'X,20F6.2)')
      WRITE (6,FORMT) (CHISQ(I,J),J=I,L)
  450 CONTINUE
      SUMDF=0.0
      DO 455 I=1,K
      SUMDF=SUMDF+ISAVE(I)
  455 CONTINUE
      DF=DF-SUMDF
      IF (SUMDF.GT.0.0) WRITE (6,460)
  460 FORMAT (/1X,'(FREQUENCIES WERE COMBINED WHERE EXPECTED VALUES WER'
     1,'E SMALL)'/)
C
      WRITE (6,465)
  465 FORMAT (//1X,'TEST OF THE NULL HYPOTHESIS THAT THE DATA FIT MODEL'
     1,' 0'/)
      WRITE (6,470) CHISUM
  470 FORMAT (' ','   CHI-SQUARED VALUE (SAMPLE) = ',F8.2)
      WRITE (6,475)
  475 FORMAT (' ','   THEORETICAL CHI-SQUARE VALUE AT THE 5% LEVEL = ')
      IF (DF-100.0) 480,480,490
  480 J=DF
      WRITE (6,485) THCHI(J)
  485 FORMAT ('+',50X,F8.2)
      GO TO 500
  490 THECHI=.5*(1.64+DSQRT((2.0*DF)-1))**2
      WRITE (6,485) THECHI
      IDF=DF
  500 WRITE (6,505) J
  505 FORMAT (' ','   DEGREES OF FREEDOM = ',I6)
      CALL CHI1(J,CHISUM)
C
      IF (K.LT.4) GO TO 571
      WRITE (6,510)
  510 FORMAT ('0',' ESTIMATED COVARIANCE AND CORRELATION COEFFICIENTS:
     1S(I)  = SURVIVAL RATE IN YEAR I')
      WRITE (6,515)
  515 FORMAT (' ',53X,'F*(I) = FIRST-YEAR RECOVERY RATE IN YEAR I')
      WRITE (6,520)
  520 FORMAT (' ',53X,'F(I)  = RECOVERY RATE IN YEAR I'/)
C
      DO 525 I=1,K
      IF (I.NE.1.AND.VARF(I).NE.0.0) CRFFX(I)=CVFFX(I)/(SDF(I)*SDFX(I))
      IF (VARS(I).NE.0.0) CRFXS(I)=CVFXS(I)/(SDFX(I)*SDS(I))
      IF (I.NE.1.AND.VARS(I).NE.0.0) CRFXPS(I)=CVFXPS(I)/(SDFX(I+1)*SDS(
     1I))
      IF (I.NE.1.AND.VARS(I).NE.0.0.AND.VARF(I).NE.0.0) CRFS(I)=CVFS(I)/
     1(SDF(I)*SDS(I))
      IF ((SDS(I)*SDS(I+1)).NE.0.0) CRSSP(I)=CVSSP(I)/(SDS(I)*SDS(I+1))
      IF ((SDF(I+1)*SDS(I)).NE.0.0) CRFPS(I)=CVFPS(I)/(SDF(I+1)*SDS(I))
      JYEAR(I)=YRONE+I-1
  525 CONTINUE
C
      WRITE (6,530)
  530 FORMAT (/' ',' I',14X,'COVAR(F*(I),F(I))',18X,'COVAR(F*(I),S(I))',
     118X,'COVAR(F*(I+1),S(I))')

      DO 540 I=1,KK
      IF (I.EQ.1) GO TO 532
      IF(I.EQ.KK.AND.MDIFF.EQ.0) GO TO 534
      WRITE (6,531)JYEAR(I),CVFFX(I),CVFXS(I),CVFXPS(I)
531   FORMAT(' ',I5,13X,F12.9,24X,F12.9,24X,F12.9)
      GO TO 540
532   WRITE (6,533) JYEAR(I),CVFXS(I),CVFXPS(I)
533   FORMAT(' ',I5,13X,'************',24X,F12.9,24X,F12.9)
      GO TO 540
534   WRITE (6,535)JYEAR(I),CVFFX(I)
535   FORMAT(' ',I5,13X,F12.9,24X,'************',24X,'************')
540   CONTINUE
C
      WRITE (6,545)
  545 FORMAT (/' ',' I',14X,'CORR(F*(I),F(I))',19X,'CORR(F*(I),S(I))',20
     1X,'CORR(F*(I+1),S(I))')
      DO 550 I=1,KK
      IF(I.EQ.1) GO TO 546
      IF(I.EQ.KK.AND.MDIFF.EQ.0) GO TO 547
      WRITE (6,531)  JYEAR(I),CRFFX(I),CRFXS(I),CRFXPS(I)
      GO TO 550
546   WRITE (6,533) JYEAR(I),CRFXS(I),CRFXPS(I)
      GO TO 550
547   WRITE (6,535) JYEAR(I),CRFFX(I)
550   CONTINUE
C
      WRITE (6,551)
  551 FORMAT (/' ',' I',14X,'COVAR(F(I),S(I)) ',18X,'COVAR(F(I+1),S(I))'
     1,17X,'COVAR(S(I),S(I+1))')
      DO 560 I=1,KK
      IF(I.EQ.1) GO TO 552
      IF(I.EQ.KK-1.AND.DIFF.EQ.0) GO TO 553
      IF (I.EQ.KK.AND.MDIFF.EQ.0) GO TO 555
      WRITE (6,531) JYEAR(I),CVFS(I),CVFPS(I),CVSSP(I)
      GO TO 560
552   WRITE (6,533) JYEAR(I),CVFPS(I),CVSSP(I)
      GO TO 560
553   WRITE (6,554) JYEAR(I),CVFS(I),CVFPS(I)
554   FORMAT(' ',I5,13X,F12.9,24X,F12.9,24X,'************')
      GO TO 560
555   WRITE (6,556) JYEAR(I)
556   FORMAT(' ',I5,13X,'************',24X,'************',24X,
     1'************')
560   CONTINUE
C
      WRITE (6,565)
  565 FORMAT (/' ',' I',14X,'CORR(F(I),S(I))  ',18X,'CORR(F(I+1),S(I)) '
     1,17X,'CORR(S(I),S(I+1))')
      DO 570 I=1,KK
      IF(I.EQ.1) GO TO 566
      IF(I.EQ.KK-1.AND.MDIFF.EQ.0) GO TO 567
      IF(I.EQ.KK.AND.MDIFF.EQ.0) GO TO 568
      WRITE (6,531) JYEAR(I),CRFS(I),CRFPS(I),CRSSP(I)
      GO TO 570
566   WRITE (6,533) JYEAR(I),CRFPS(I),CRSSP(I)
      GO TO 570
567   WRITE (6,554) JYEAR(I),CRFS(I),CRFPS(I)
      GO TO 570
568   WRITE (6,556) JYEAR(I)
570   CONTINUE
  571 WRITE (6,245)
      RETURN
      END
      SUBROUTINE MOD0W1
      WRITE (6,10)
   10 FORMAT ('0','MODEL 0')
      WRITE (6,15)
   15 FORMAT (' ','********'/)
      WRITE (6,20)
   20 FORMAT (1X,'ANALYSIS UNDER THE ASSUMPTIONS OF TIME-SPECIFIC SURVI'
     1,'VAL AND RECOVERY RATES WHERE THE FIRST-YEAR RECOVERY RATES ARE')
      WRITE (6,25)
   25 FORMAT (1X,'DIFFERENT FROM RECOVERY RATES OF PREVIOUSLY-BANDED CO'
     1,'HORTS.'/)
      WRITE (6,30)
   30 FORMAT (1X,'THIS MODEL IS AN EXTENSION OF THE SEBER-ROBSON-YOUNGS'
     1,' MODEL (MODEL 1) IN THAT THE FIRST-YEAR RECOVERY RATES ARE ',
     2 'ALLOWED TO')
      WRITE (6,35)
   35 FORMAT (1X,'DIFFER.  THIS MODEL IS USEFUL IF THE BAND REPORTING '
     1,'RATE IS QUITE DIFFERENT THE FIRST YEAR AFTER BANDING (E.G.,',
     2 ' BAND')
      WRITE (6,40)
   40 FORMAT (1X,'COLLECTING ACTIVITIES CONDUCTED BY CONSERVATION AGENC'
     1,'Y PERSONNEL OFTEN AFFECT PRIMARILY THE FIRST-YEAR RECOVERY',
     2' RATES).'/)
      WRITE (6,45)
   45 FORMAT (1X,'(MODEL DEVELOPED BY BROWNIE (1974. CORNELL BIOMETRICS'
     1,' UNIT PAPER 535-M)).  REFER TO BROWNIE (1974) FOR A DISCUSSION',
     2' OF')
      WRITE (6,50)
   50 FORMAT (1X,'THE USE OF THIS MODEL VS. MODEL 1.'//)
      WRITE (6,55)
   55 FORMAT (' ','SPECIFICALLY, THE MODEL STRUCTURE IS: '/)
      WRITE (6,60)
   60 FORMAT (1X,'    F*(1)     S(1)F(2)     S(1)S(2)F(3)     S(1)S(2)S'
     1,'(3)F(4)     S(1)S(2)S(3)S(4)F(5)')
      WRITE (6,65)
   65 FORMAT (1X,'                 F*(2)         S(2)F(3)         S(2)S'
     1,'(3)F(4)         S(2)S(3)S(4)F(5)')
      WRITE (6,70)
   70 FORMAT (1X,'                                  F*(3)             S'
     1,'(3)F(4)             S(3)S(4)F(5)')
      WRITE (6,75)
   75 FORMAT (1X,'                                                     '
     1,'  F*(4)                 S(4)F(5)'///)
      RETURN
      END