matlab程序仿真画图报错，求解
syms x; %指曲率半径r/r0
syms m %指半径r
w=0:0.01:10; %指束腰b/r0
p=632.8*10^(-9); %指波长
z=10; %指传输距离
c=3*10^8; %指光速
u=0.185*((p^2)/(z*2.5*10^(-17)))^(3/5); %指参数r0
x=m/u;
l=(1/2)*w*u; %指w0
r=(1/2)*((2*pi)/p).*(l.^2);%指瑞利长度Z（R）
d=l.*sqrt((1+(z./r).^2));%指w（z）
T0=0.35; %指盐度
t=10;
m0=1;
syms v
q1=T0/(2*sqrt(pi)).*int(exp(-(v.^2.*T0.^2)/4+1i.*v.*t)*m.^m0./d.^(m0+1).*exp((1i*(v+l)*z)./c).*exp(-(1-1i*(z./r)).*m.^2./d.^2).*exp(-1i*(m0+1).*atan(z./r)),v);
v=1;
z1=0;
syms y
a=exp(-6.88*2.^(2/3)*x.^(5/3)*abs(sin(y/2)).^(5/3));
b1=exp(-1i*z1*y);
s1 = int(a*b1,0,2*pi);
c1=(1/(2*pi))*s1;
n1 = integral( matlabFunction( (abs(q1)).^2*m*c1 ), 0, inf, 'ArrayValued', true );
g1=n1/0.16;
plot(w,g1)
hold on
说一下积分变量，这个q1是对v进行不定积分，这个s1是对y进行定积分，最后这个n1是对m进行定积分，v，y，m都是有自己意义的参数，不能混淆，而且q1和s1积分后都含有m，所以会带入n1中对m进行最后的积分
