symsum does not work with symprod
1 回表示 (過去 30 日間)
古いコメントを表示
I'm trying to avaluate the following:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1311495/image.png)
syms j k
double(symsum((-symprod((j+0.3-1)/(j), j, 1,k))^2, k, 2, Inf))
It returns: -0.0806, which is nonsensical, since this is a sum of squared numbers.
When I evaluate the same in Wolfram Alpha using "sum( (product (j+0.3-1)/j, j=1..k)^2, k, 2, infinity)", I get 0.226456, which seems to be correct, since Matlab seems to converge to this number up to k=100000 (and then it crushes). How can I make it work? Thanks.
0 件のコメント
採用された回答
Walter Roberson
2023 年 3 月 1 日
syms j k
double(symsum((-symprod((j+0.3-1)/(j), j, 1,k))^2, k, 2, Inf))
So it works in current releases.
I am checking back to see if I can figure out when it was fixed.
3 件のコメント
Walter Roberson
2023 年 3 月 1 日
Works in R2018b.
It appears that R2018a is the last release that has this problem.
It is a problem with the symsum() part, by the way; it is returning a hypergeom() minus a constant when it should be returning an expression in gamma() minus that same constant.
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Calculus についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!