我想解下面的方程
我在函数的定义上做错了,但我仍然不明白
function F=myairyint(x)
F=integral(@(x)airy,x,1000)
end
functi2=@(x) myairyint;
x0=-1:-1:-15;
fsolve(functi2,x0)
为什么不起作用?
你的积分似乎没有根源[-15 -1]
并收敛于零。我已经通过绘图测试了它:
x0 = linspace(-20,20);
airyint = arrayfun(@(x0) integral(@(n) airy(n), x0, inf), x0);
plot(x0, airyint);
我也通过测试过fzero()
例如:
f = @(x) integral(@(n) airy(n), x, inf);
fzero(f, 0)
输出如预期:
ans =
115.8524
但似乎每个实数x > 115.8524看起来像这个积分方程的根。
额外的:
fzero()
是一个试图通过查看两个动态区间之间的函数值来查找根的函数。如果它捕获负值和正值,则它正在缩小捕获根的间隔。但是有一个错误率,因为你也可以从这个例子中观察到。仅当 x 趋于无穷大时,该方程才为零,这意味着该方程没有实根。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)