akb(i,j,bi,bj)=exp((-8966.90 - 2890.53*sqrts - 77.942*s +
& (148.0248 + 137.1942*sqrts + 1.62142*s) +
& (-24.4344 - 25.085*sqrts - 0.2474*s) *
& dlogtk + 0.053105*sqrts*tk)
& (-106.736*invtk + 0.69171)*sqrts +
& (-160.340*invtk + 1.3566) * sqrts +
& sqrts + (-44.99486*invtk - 0.09984) * s)
& sqrts - 0.01615 * s)
tmpa2 = +(-0.77712 + (0.0028426*tk) + (178.34/tk) )*sqrts
_RL sqrts
sqrts=sqrt(s)
akb(i,j,bi,bj)=exp((-8966.90 _d 0- 2890.53 _d 0*sqrts -
& (148.0248 _d 0 + 137.1942 _d 0*sqrts + 1.62142 _d 0*s) +
& (-24.4344 _d 0 - 25.085 _d 0*sqrts - 0.2474 _d 0*s) *
& dlogtk + 0.053105 _d 0*sqrts*tk)
& (-106.736 _d 0*invtk + 0.69171 _d 0)*sqrts +
& (-160.340 _d 0*invtk + 1.3566 _d 0) * sqrts +
& sqrts + (-44.99486 _d 0*invtk - 0.09984 _d 0) * s)
& * sqrts - 0.01615 _d 0 * s)