Ok, I am not a very experienced user, but I am mystified by this.  The essence of the problem is that when entering a range of values in an integration, I get a different answer if the range is given as 0..1 to if it is given as 0.0..1 or 0..1.0. Here is the program

### dgig
dgig := proc(psi, chi, lambda, t)
    (psi/chi)^(lambda/2)*t^(lambda - 1)*exp(-(psi*t+chi/t)/2)/\
    (2*BesselK(lambda,sqrt(chi*psi))):
end: ## End of dgig