非常慢的无限系列与matlab中的symsum

Question

我有以下无限系列我需要评估。目前，根据输入参数（kappa，sigmav，thetav），此函数需要很长时间才能计算。有没有可能加快速度？耗时的表达式是symsum（...，0，Inf）

 j = sym('j', 'positive');
 c =  1/((1/2*kappav)*(sigmav^2)*(1-exp(-kappav*deltat)));
 lambda = 2*c*Vt*exp(-kappav*deltat);
 v = 4*thetav/sigmav^2;
 ez = 2*c*exp(Vtt);
 % 
 p = ((exp(-1/2*(ez+lambda))*ez^(1/2*v)) / ( 2^(0.5*v)))* ... 
   vpa( symsum(   ((1/4*ez*lambda)^j)/ (gamma(1/2*v+j)*factorial(j)) ,j, 0, Inf  ));

PS。这是non-central chi-square distribution的pdf，稍作修改，因此我无法使用ncx2pdf。我已经尝试https://stackoverflow.com/a/15966126/321749，这只会导致无穷大。