fixed error in bessel_k.c : if fabs(val_new)<eps
use taylor approximation in bessel_k.c for nu=0, x and y small
included upper bound for modified bessel function
Merge branch 'master' of github.com:scafacos/scafacos