次を使用中のエラー: symengine Division by zero.

以下の式のε（コード内ではe）を求めたいのですが，

下のようなコードを組むと，

syms tdvs s 
m=4;x=1;p1=5;p2=5;p3=5;bw=0.58/2;n=40;a=0.3387;
Rtdvs=((abs(cos(pi*m/4*tdvs*x)))^(p2)+(abs(sin(pi*m/4*tdvs*x)))^(p3))^(-1/p1);
ytvs=-Rtdvs*sin(pi*tdvs*x);
dytvs=diff(ytvs,tdvs);
dynew = subs(dytvs,tdvs,cos(s*pi/(2*n+1)));
w=(sin(s*pi/(2*n+1)))^2/(2*n+1);
e=-4/bw^2*symsum(w*a*dynew,s,1,n)

Error using symengine
Division by zero.

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

上のようなエラーが出てしまいました．

どうプログラムを変えれば、上手く動きますか？