c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc       parameters of fit3 for Fe                      ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 )  =  0.04192d0
        a( 1 )  = -0.01768d0
        a( 2 )  = -0.00007d0
        a( 3 )  = -0.00241d0
        a( 4 )  = -0.00080d0
        b( 1 )  = -0.01705d0
        b( 2 )  = -0.00268d0
        b( 3 )  = -0.00141d0
        b( 4 )  =  0.d0
        c       =  0.00818d0
        d       =  0.13629d0
        e( 0 )  =  0.09256d0
        e( 1 )  = -0.03290d0
        e( 2 )  = -0.00523d0
        e( 3 )  = -0.00539d0
        e( 4 )  = -0.00276d0
        f( 1 )  = -0.02574d0
        f( 2 )  = -0.00630d0
        f( 3 )  = -0.00285d0
        f( 4 )  =  0.d0
        g       =  0.00022d0
        h       =  0.31871d0
        i( 0 )  =  0.00977d0
        i( 1 )  = -0.00653d0
        i( 2 )  =  0.00171d0
        i( 3 )  = -0.00024d0
        i( 4 )  =  0.00017d0
        j( 1 )  =  0.00869d0
        j( 2 )  = -0.00323d0
        j( 3 )  =  0.00075d0
        j( 4 )  =  0.d0
        k       = -0.00439d0
        l       =  0.02766d0
        p( 0 )  =  0.01464d0
        p( 1 )  = -0.00957d0
        p( 2 )  =  0.00222d0
        p( 3 )  = -0.00022d0
        p( 4 )  =  0.00025d0
        q( 1 )  =  0.01867d0
        q( 2 )  = -0.00615d0
        q( 3 )  =  0.00133d0
        q( 4 )  =  0.d0
        r       = -0.01023d0
        s       =  0.06442d0
        ap( 0 ) =  0.02584d0
        ap( 1 ) = -0.00894d0
        ap( 2 ) = -0.00097d0
        ap( 3 ) = -0.00125d0
        ap( 4 ) = -0.00045d0
        bp( 1 ) = -0.00739d0
        bp( 2 ) = -0.00190d0
        bp( 3 ) = -0.00062d0
        bp( 4 ) =  0.d0
        cp      =  0.00167d0
        dp      =  0.09950d0
        ip( 0 ) =  0.00480d0
        ip( 1 ) = -0.00271d0
        ip( 2 ) =  0.00047d0
        ip( 3 ) =  0.00001d0
        ip( 4 ) = -0.00001d0
        jp( 1 ) =  0.00604d0
        jp( 2 ) = -0.00197d0
        jp( 3 ) =  0.00044d0
        jp( 4 ) =  0.d0
        kp      = -0.00320d0
        lp      =  0.02012d0
        alpha0  =    0.6798d0
        alpha1  =   12.7527d0
        alpha2  = -140.1800d0
        alpha3  =  268.8290d0
        beta0   =    0.7783d0
        beta1   =   10.2315d0
        beta2   = -124.2640d0
        beta3   =  241.3060d0
        alphap0 =    0.2847d0
        alphap1 =  -13.0828d0
        alphap2 =  192.5030d0
        alphap3 = -543.1480d0
        betap0  =    0.3221d0
        betap1  =  -13.8640d0
        betap2  =  201.0700d0
        betap3  = -573.1160d0