积分积不出结果，画不出图

积分积不出结果，画不出图

clear;

%反射率
r1=1;
r2=2^0.5;
R=1;
%台阶高度
u1=1;
u2=-1;
%离焦量
us=0;
uz1=us-u1;
uz2=us-u2;
%偏移量
ud=1.27;
Ud=abs(ud);

a0=asin(0.4);
s0=1/(4*(sin(a0/2))^2);
syms s;
ss=s+2*s0;
syms m;
l=m;
L=abs(l);
syms v;

%第一类贝塞尔函数
J0=besselj(1,v);

f=J0.(1/2Ud.L.(ss-l^2/4).^0.5).*r1.*sin(v.m)./m.exp(1i(uz1-ud/2).s);
result_1=int(f,s,l.^2/4,l-L.(1-L/2));
result_2=int(int(f,s,l.^2/4,l-L.(1-L/2)),m,-2,2);
% result_3=subs(result_2,v,1);
% result_4=vpa(result_3,5);
fplot(@(v)result_2,[-10,10]);

警告: 函数处理数组输入时行为异常。要改善性能，请将您的函数正确向量化，以返回大小和形状与输入参数相同的输出。

位置：matlab.graphics.function.FunctionLine>getFunction
位置: matlab.graphics.function/FunctionLine/updateFunction
位置: matlab.graphics.function/FunctionLine/set.Function_I
位置: matlab.graphics.function/FunctionLine/set.Function
位置: matlab.graphics.function.FunctionLine
位置: fplot>singleFplot (第 245 行)
位置: fplot>@(f)singleFplot(cax,{f},limits,extraOpts,args) (第 200 行)
位置: fplot>vectorizeFplot (第 200 行)
位置: fplot (第 166 行)
警告: 更新 FunctionLine 时出错。

计算 FunctionLine 更新中的函数时报告了以下错误: Unable to convert expression containing symbolic variables into double array. Apply 'subs' function first
to substitute values for variables.

单步调试或加日志