>> a=0.001;
>> M=1032000;
>> x=0.0025;
>> z=-0.0008:0.00005:0.0008;
>> f=@(n,z0)(a.*cos(n).*(z-z0))./(sqrt((x-a.*cos(n)).^2+(0-a.*sin(n)).^2+(z-z0).^2)).^3;
>> M=1032000;
>> x=0.0025;
>> z=-0.0008:0.00005:0.0008;
>> f=@(n,z0)(a.*cos(n).*(z-z0))./(sqrt((x-a.*cos(n)).^2+(0-a.*sin(n)).^2+(z-z0).^2)).^3;
>> q=quadv(@(n)quadv(@(z0)f(2*pi,z0),-1,1),0,2*pi)
对f函数进行二重积分,积分对象为z0,n,而且z也为一个矩阵,调用quadv这个函数哪里写错了呢,有其他的表示方法么?本问题的公式在图中