>> f = @(a) dblquad(@(x,z) -0.004*x.*(1-(z/6.25).^2).*exp(0.05*z*sec(a)).*(cos(0.05*x.*(sec(a)))+i*sin(0.05*x.*sec(a))),0,100,-6.25,0);>> a=1;>> f(a)ans = 0.6541 -14.6963i>> a = 3;>> f(a)ans = 40.0774 -19.3193i
a