Real-polygamma-func
evaluates an approximation to the kth derivative of the psi function by
where x is real with x≠0, -1, -2, ... and k=0,1,......6.
For more information please review the s14aec function in the NAG document.