Hello All:
I tried two methods as follows, but did not work correctly. Your help will be greatly appreciated.
R=2;
photonnumber=100000;
r=5;
meff=0.9;
F=0.4;
NA=pi/2
D=1.2;
1.
syms beta ceta z
bdetectpower=3*photonnumber/(4*pi*R^3)*int(int(int(exp(-sqrt(r^2+z^2-2*r*z*cos...
