Matlab程序中求和表达式symsum求和边界出错

Question

问题遇到的现象和发生背景

原式要求解的是下限为0，上限为正负pai的求和，如图所示，三个方程组中其一的表达式

问题相关代码，请勿粘贴截图

其中自己写的求解这部分的表达式代码(MATLAB语言）如下：

h1=(((((((sind(a0)+(x1/A)+Ri*(x3/A)*cos(v))^2)+((cosd(a0)+(x2/A)*cos(v))^2))^0.5)-1)^1.5)*(sind(a0)+(x1/A)+Ri*(x3/A)*cos(v)))/((((sind(a0)+(x1/A)+Ri*(x3/A)*cos(v))^2)+((cosd(a0)+(x2/A)*cos(v))^2))^0.5);
ZS1=symsum(h1,v,-pi,+pi);
f1=Fa-Kn*(A^1.5)*ZS1;

运行结果及报错内容

运行后这部分出错代码：
错误使用 symengine
Invalid bounds for summation index.

出错 sym/symsum (line 70)
rSym = mupadmex('symobj::map',fsym.s,'symobj::symsum',x.s,a.s,b.s);

出错 fun (line 36)
ZS1=symsum(h1,v,-pi,+pi);%v是滚珠方位角

我的解答思路和尝试过的方法

我一开始想symsum直接求解，但是尝试解决后，边界上下限必须是整数，不能解决我这个求和上下限。

我想要达到的结果

将这方程式中的求和部分表达出来

Answer 1

如果是轴承问题，那么这个求和代表的意义应该是沿着轴承一周分布的每个滚珠的作用力（或其它物理量）进行求和，那么虽然公式上写的是对角度进行求和，实际上是对每个滚珠进行求和，将每个滚珠的方位角作为变量带入求出求和符号内的函数值后，直接将各个滚珠的函数值相加即可。