Simpleson积分
 

请教哪位大虾帮我解决风功率密度的定积分问题！
我用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(:).';
请教大虾，感谢！！！