2017-02-08 28 views
0

我有一个递归方程的量w1 [k,n],w2 [k,n]其中k是正整数或负整数,n是正整数。公式如下:Maxima CAS - 如何设置东西(为零)

eq1_1 : w1[ K, N ] + R*sum(sum(w1[ K-k, N-n ]*'diff(w2[ k, n-1 ], R) + w2[ K-k, N-n ]*'diff(w2[ k, n-1 ], z) + %i*k*(w1[ K-k, N-n ]*w2[ k, n ]*'diff(S, R) + w2[ K-k, N-n ]*w2[ k, n ]*'diff(S, z)) , k, -n-1, n+1), n, 0, N) = 0; 

我写的规则

rule0: [N=0, K=0] 

写的第一项。我只想写它,而不是解决它。我代替这个规则,设置

simpsum:true 

,我得到的输出方程(很抱歉,我不能在这里复制)与像W1条件[-1,-1],它应该等于零,因为n是正数。我试图将它设置为零,subst,ratsubsttellimp但它们都没有工作(也许我使用了错误的语法)。有没有办法将这些条件设置为零,或甚至更好地告诉最大值在执行任何计算之前不要写入它们?

+0

什么是'W1 [K,N]''和W2 [K,N]'的定义? –

回答

0
w2[-1,0] : 0$ /* repeat as needed for other specific elements */ 
subst(rule0,eq1_1); 
ev(%); 

或,等价地,

w2[-1,0] : 0$ 
subst(rule0,eq1_1),eval;