c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc       parameters of fit3 for Si                      ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 )  =  0.03318d0
        a( 1 )  = -0.01218d0
        a( 2 )  = -0.00145d0
        a( 3 )  = -0.00197d0
        a( 4 )  = -0.00099d0
        b( 1 )  = -0.00940d0
        b( 2 )  = -0.00281d0
        b( 3 )  = -0.00093d0
        b( 4 )  =  0.d0
        c       =  0.00123d0
        d       =  0.13959d0
        e( 0 )  =  0.06964d0
        e( 1 )  = -0.02211d0
        e( 2 )  = -0.00656d0
        e( 3 )  = -0.00379d0
        e( 4 )  = -0.00236d0
        f( 1 )  = -0.00863d0
        f( 2 )  = -0.00450d0
        f( 3 )  = -0.00162d0
        f( 4 )  =  0.d0
        g       = -0.01325d0
        h       =  0.29081d0
        i( 0 )  =  0.00573d0
        i( 1 )  = -0.00370d0
        i( 2 )  =  0.00076d0
        i( 3 )  =  0.00000d0
        i( 4 )  =  0.00007d0
        j( 1 )  =  0.00826d0
        j( 2 )  = -0.00279d0
        j( 3 )  =  0.00065d0
        j( 4 )  =  0.d0
        k       = -0.00448d0
        l       =  0.02820d0
        p( 0 )  =  0.00636d0
        p( 1 )  = -0.00389d0
        p( 2 )  =  0.00050d0
        p( 3 )  =  0.00013d0
        p( 4 )  =  0.00008d0
        q( 1 )  =  0.01587d0
        q( 2 )  = -0.00465d0
        q( 3 )  =  0.00100d0
        q( 4 )  =  0.d0
        r       = -0.00931d0
        s       =  0.05857d0
        ap( 0 ) =  0.01767d0
        ap( 1 ) = -0.00455d0
        ap( 2 ) = -0.00177d0
        ap( 3 ) = -0.00084d0
        ap( 4 ) = -0.00041d0
        bp( 1 ) = -0.00244d0
        bp( 2 ) = -0.00152d0
        bp( 3 ) = -0.00044d0
        bp( 4 ) =  0.d0
        cp      = -0.00239d0
        dp      =  0.09466d0
        ip( 0 ) =  0.00213d0
        ip( 1 ) = -0.00094d0
        ip( 2 ) = -0.00003d0
        ip( 3 ) =  0.00010d0
        ip( 4 ) = -0.00003d0
        jp( 1 ) =  0.00520d0
        jp( 2 ) = -0.00146d0
        jp( 3 ) =  0.00030d0
        jp( 4 ) =  0.d0
        kp      = -0.00303d0
        lp      =  0.01906d0
        alpha0  =    0.5502d0
        alpha1  =   15.4934d0
        alpha2  = -148.0760d0
        alpha3  =  250.3260d0
        beta0   =    0.6380d0
        beta1   =   13.3510d0
        beta2   = -136.2650d0
        beta3   =  235.7860d0
        alphap0 =    0.2448d0
        alphap1 =  -12.1275d0
        alphap2 =  173.0720d0
        alphap3 = -457.7020d0
        betap0  =    0.2809d0
        betap1  =  -12.7225d0
        betap2  =  177.8420d0
        betap3  = -472.0400d0