The analytic expression from ZF 2009 is coded, but a factor of 0.01 is multiplied with, to overcome negative p(E) sweeped by the code. This is just a shortcut, but exact prescription will be found soon.