第一步
将
function y = myfun(x)
y = 1./(x.*(x.^2 - 3*x + 2 ).^(1/3));
保存为M文件.(文件→新建→M文件)
第二步
将
warning off all
Q = quad(@myfun,2,100000)
Q = quad(@myfun,2,10000000)
Q = quad(@myfun,2,1000000000)
warning on
粘贴至(命令窗口)
结果显示,
当上限为100000,10000000,1000000000时,
定积分的值为x=1.4389,1.4396,1.4396.
因此,可以将1.4396作为此定积分的值.