SUBROUTINE LAGRAN(X0,X1,X2,Y0,Y1,Y2,X,Y)                         
C     =========================================                               
C
C ***                                                                   
C *** Lagrange interpolation for three points                           
C ***
      INCLUDE 'IMPLIC.FOR'
C                                                                        
      XL0=(X-X1)*(X-X2)/(X0-X1)/(X0-X2)                                 
      XL1=(X-X0)*(X-X2)/(X1-X0)/(X1-X2)                                 
      XL2=(X-X0)*(X-X1)/(X2-X0)/(X2-X1)                                 
                                                                        
      Y=Y0*XL0+Y1*XL1+Y2*XL2                                            
                                                                        
      RETURN                                                            
      END