I am trying to compute the integral from 0 to infinity of the following PDF:
S=@(r) 2*pi*lambdab.*r.*exp(-pi*lambdab.*r).*exp(-2*pi*lambdab.*integral(@(u) u.*exp(-pi*lambdab.*u),0,r,'ArrayValued',true));
Logically, I need to find 1 for every >0. However, when is greater than 0.2, I surprisingly get smaller values than 1 for the integral of the PDF.
Any help please!