Fix bug in Taylor series for Bessel K
For small arg(z) and abs(z)> ~50, K(m,z) behaves very erratically, and has a large absolute and relative error compared to the implementation in scipy. Investigate the cause of this and fix the implementation. Might be an idea to look at the details of the Amos implementation (used by scipy), it is done a bit different.