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)