2

次のシリーズの 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'

この問題に対するより簡単なアプローチはありますか?

4

0 に答える 0