SUBROUTINE FPPPTZ
*-- Author :  R. Henderson
      SUBROUTINE FPPPTZ(PERSEG,Z)
C                                                                                                             
C---  Propagates planars error matrix from z=0 to z=z                                                         
C                                                                                                             
      DOUBLE PRECISION COVAR,PERSEG                                     
C---                                                                                                          
      DIMENSION PERSEG(4,4)                                             
      DIMENSION COVAR(4,4) , DERIV(4,4) , WM(4,4)                       
C                                                                                                             
C---  Construct the derivative matrix                                                                         
C                                                                                                             
      CALL VZERO(DERIV,16)                                                                             
      DO 10 I = 1,4                                                     
      DERIV(I,I) = 1.0                                                  
   10 CONTINUE                                                          
      DERIV(1,3) = Z                                                    
      DERIV(2,4) = Z                                                    
C                                                                                                             
C---  Propagate PERSEG FROM Z=0 TO Z                                                                          
C                                                                                                             
      CALL VZERO(WM,16)                                                                                
C                                                                                                             
      DO 24 I = 1,4                                                     
        DO 25 J = 1,4                                                   
         DO 26 K = 1,4                                                  
         WM(I,J) = WM(I,J) + PERSEG(I,K) * DERIV(J,K)                   
   26    CONTINUE                                                       
   25   CONTINUE                                                        
   24 CONTINUE                                                          
C---                                                                                                          
      CALL VZERO(COVAR,32)                                                                             
      DO 54 M = 1,4                                                     
        DO 55 K = 1,4                                                   
         DO 56 I = 1,4                                                  
         COVAR(M,K) = COVAR(M,K) + DERIV(M,I) * WM(I,K)                 
   56    CONTINUE                                                       
   55   CONTINUE                                                        
   54 CONTINUE                                                          
C---                                                                                                          
      CALL UCOPY(COVAR,PERSEG,32)                                                                      
C---                                                                                                          
      RETURN                                                            
      END                                                               
*