The analytic expression from ZF 2009 has been coded, but a factor of 0.01 is multiplied with nu*rho to overcome the negative p(E) that the code sweeps. This is just a numerical shortcut and exact prescription has to be found soon.