(%i1) besselexpand:true$ (%i2)limit(bessel_k(1/6,x)/(sqrt(Pi/(2*x)) * exp(-x)) ,x,inf) (%o2)(sqrt(2)*(limit(bessel_k(1/6,x)*sqrt(x)*%e^x,x,inf)))/sqrt(Pi) which I think should be 1. so (%i3) limit(bessel_k(1/6,x)/(sqrt(Pi/(2*x)) * exp(-x)) ,x,inf) (%o3) 1 would be much nicer. Rich