先建立了一个函数:
function y= f(x)
y= x^(1/2)*log(x);
又建立函数是这样的:
function s= trapz( a, b, n)
h= ( b- a)/n; s= 0;
for k= 1: ( n- 1)
x= a+ h*k;
s= s+ f( x) ;
end
s=h*(f(a) + f(b))/2+ h*s;
然后输入trapz( 1, 0, 3)时]
trapz( 1, 0, 3)
ans =
NaN
哎,谁能告诉我该怎么用,才能打印出我建立的那个函数所求出来的值啊。。。。额。。我建立那个函数本意是建立用复合梯形方法求积分的。。
------解决方案--------------------------------------------------------
最后一步计算f(b)的时候,f(x)中的log(x)项造成的,你输入的参数b=0。