[nb,N]=size(dvh);
V_total=3.6;
rou=10^7;
a0=0.3;
st_a=0.1;
syms a
m=0;
for i=1:nb
m=m+V_total*dvh(i,2)*exp(-a*dvh(i,1));
end
g=exp(-rou*m);
t=g*normpdf(a,a0,st_a);
s=0;
for a=a0-3*st_a:0.006*st_a:a0+3*st_a
s=s+0.006*st_a*subs(t);
end %dvh是一个多行2列的数组，这个代码如何改成不包含syms命令，比如使用quad函数。