Question: Numerical integration problem

I have a problem with the numerical integration in Maple. Any help or suggestion will be greatly appreciated!

 

K[1] := 1.0;
K[2] := 1;
K[3] := 0.16;
nu:=-5.92;
alpha:=8;
y:=0.01;

dF:=proc(z,y)
local k,q,h;
q:=BesselI(alpha,z)/BesselI(0.0,z);
h:=2*(K[1]*q-y)/(K[3]+(K[3]^2+4*K[2]*z^2*K[1]*q-4*K[2]*z^2*y)^(1/2));
k:=-(K[3]^2+4*K[2]*z^2*(K[1]*q-y))^(-1/2);
z^nu*h^(nu-1)*exp(-1/2/h-z^2*h/2)*k;
end proc:

evalf(int(1/z*dF(z,y)*BesselI(lambda,z),z=8..100));

For some reason, Maple does not seem to be able to return any value for this integral. I waited for very long time without getting anything. But if you run the plot the integrand over this interval. The function looks perfectly smooth.

plot(1/z*dF(z,y)*BesselI(lambda,z),z=8..100);

Does anyone have an idea what's wrong with this work? Many many thanks for your help!

Please Wait...