Hi every one,
I'm trying to find a numeric solve for an integral, and got some weird answer.
so i tried with another function and got a different answer.
i suppose one of them is the right answer, my bet is the trapz answer is the right one. but I'm not sure.
can somebody explain where i am wrong? next thing i need to solve an 2d-integral using quad2d, and i have to know how to do it right.
10x
here is a running code:
clear Fi1=@(ro,z,theta1,k)exp(k.*z.*cos(theta1).*1i).*sqrt(cos(theta1)).* sin(theta1).^2.*besselj(1.0,k.*ro.*sin(theta1)); NA=0.95; ni=1; alpha=asind(NA/ni) alpha=asin(NA/ni); lambda=600e-9; r=0.5e-6; z=0; V=linspace(0,alpha,300); k=2*pi/lambda; I1_1=trapz( Fi1(r,z,V,k), V) I1_2=quad( @(theta1) Fi1(r,z,theta1,k), 0,alpha)
Best Answer