次のシリーズの n=1 から n=infinity までの合計を見つけようとしています:
A(n)=T/(n+1)*{a1*T^2*A(n-3)/((n-1)*n) + a2*T*A(n-2)/n + a3*A(n-1)}
と
A(-1)=A(0)=0
A(1)=T*a1/2
ムパッドを使用。出来ますか?
次のコードを試しました:
A := proc(n)
begin
if(n<=0) then return(0) end_if:
if(n=1) then return(a_1*T/2) end_if:
return((T/(n+1))*(a_1*T^2*A(n-3)/((n-1)*n)
+a_2*T*A(n-2)/n
+a_3*A(n-1)))
end_proc:
これは、シリーズ A の単一要素に対しては正常に機能しますが、n=1..2 の部分和を見つけようとすると、次のようになります。
sum(A(n),n=1..2)
mupad は次のエラーを返します。
Error: Can't evaluate to boolean [_leequal];
during evaluation of 'A'
この問題に対するより簡単なアプローチはありますか?