2017-02-08 82 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;