      SUBROUTINE SMOOTH(Y,MU,X)
************************************************************************
*                                                                      *
*     CHKS smoothing function                                          *
*                                                                      *
************************************************************************
*
      REAL*8 X,Y,MU
*
      Y=(X+DSQRT(X*X+4.0*MU*MU))/2.0
*
      RETURN
      END
*
*
      SUBROUTINE SMOOTHM(Y,MU,X)
************************************************************************
*                                                                      *
*     derivative w.r.t. epsilon                                        *
*                                                                      *
************************************************************************
*
      REAL*8 X,Y,MU
*
      Y=2.0*MU/DSQRT(X*X+4.0*MU*MU)
*
      RETURN
      END
*
*
      SUBROUTINE SMOOTHX(Y,MU,X)
************************************************************************
*                                                                      *
*     derivarive w.r.t. x                                              *
*                                                                      *
************************************************************************
*
      REAL*8 X,Y,MU
*
      Y=(1.0+X/DSQRT(X*X+4.0*MU*MU))/2.0
*
      RETURN
      END
*
*
