Грешка при използване на mupadmex, Грешка в sym/symsum (Оценка на определен интеграл)

Опитвам се да оценя определен интеграл, използвайки следния код:

m=4;
t=4;
n=5;
P=10;
reset(symengine)
syms x;
syms k;

mean=int(log(1+P*x/t)*symsum( factorial(k)*x^(n-m)*exp(-x)/factorial(k+n-m)*feval(symengine,'laguerreL',n-m,k,x)*feval(symengine,'laguerreL',n-m,k,x),k,[0,m-1]) ,x,0,inf)

Но дава следната грешка:

  Error using mupadmex
  Error in MuPAD command: The number of arguments is incorrect. [has]
  Evaluating: sum

  Error in sym/symsum (line 114)
  rSym = mupadmex('symobj::map',fsym.s,'symobj::symsum',x.s,a.s,b.s);

  Error in Untitled2 (line 32)
  mean=int(log(1+P*x/t)*symsum( factorial(k)*x^(n-m)*exp(-x)/factorial(k+n-m)*feval(symengine,'laguerreL',n-m,k,x)*feval(symengine,'laguerreL',n-m,k,x),k,[0,m-1]) ,x,0,inf)

Въпреки това, когато оценявам интеграла без някои функции:

mean=int(log(1+P*x/t)*symsum( feval(symengine,'laguerreL',n-m,k,x)*feval(symengine,'laguerreL',n-m,k,x),k,[0,m-1]) ,x,0,inf)

Получавам отговор, който е Inf. Някаква идея?


matlab mupad
person tam    schedule 25.07.2015    source източник

Отговори (1)


arrow_upward
2
arrow_downward

Намерих решението. тя се състои в използване на gamma(k+n-m+1) за изчисляване (k+n-m)!, тъй като има грешка при извършване на факториела на символна променлива.

person tam    schedule 25.08.2015