我想对一个积分进行数值求解,该积分取决于其先前的路径。积分 7.14 从这里 http://www.mathworks.com/matlabcentral/answers/uploaded_files/8998/1.pdf.
使用下面的代码我收到错误。怎么了?我用这段代码求解这个积分吗?
“使用积分计算/最终输入检查时出错(第 516 行)函数的输出必须与输入的大小相同。如果 FUN 是数组值被积数,请将 'ArrayValued' 选项设置为 true。”
谢谢你。
z = 0:1/1000:10^-2
dz = z(2)-z(1); %integration step
sigma = 1;
q=0; %Integral value at z=0
for rr = 1:length(z)
fun = @(z) sigma*((10^4-exp(q))/((10^4+exp(q))));
q = integral(fun,0,z(rr))*dz;
end
q
更换线路:
q = integral(fun,0,z(rr))*dz;
与:
q = integral(fun,0,z(rr),'ArrayValued',true)*dz;