假设y的范围为[-10,10],则获取F与y关系的代码及图像如下所示:
>> fun = @(t,y,z)-20*sin(t).^2./((20*sin(t)).^2+(y-20*sin(t).^2).^2+(0.7-z).^2).^1.5; % 包含参数的被积函数
>> y = -10:0.2:10; % 假设参数y的范围[-10,10]
>> F = arrayfun(@(y) integral2(@(t,z)fun(t,y,z),0,2*pi,-4,0),y); % 计算每一个y对应的F值
>> plot(y,F)