rhfeng

5 Reputation

2 Badges

13 years, 325 days

MaplePrimes Activity


These are questions asked by rhfeng

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:

Page 1 of 1