登录论坛

查看完整版本 : Simpleson积分


yyc1048
2009-10-22, 16:00
请教哪位大虾帮我解决风功率密度的定积分问题!
我用Simpleson数值积分时,遇到如下问题:
>> fun = inline('v^3*(v/3.297)^(1.053-1)*exp(-(v/3.297)^1.053)');
>> q2=quadl(fun,3,25)
结果:
??? Error using ==> inlineeval at 15
Error in inline expression ==>
v^3*(v/3.297)^(1.053-1)*exp(-(v/3.297)^1.053)
Matrix must be square.

Error in ==> inline.feval at 36
INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr,
INLINE_OBJ_.expr);

Error in ==> quadl at 70
y = feval(f,x,varargin{:}); y = y(:).';
请教大虾,感谢!!!