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
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