[問題]請問 程式跑出來的答案error已回收
前輩們好：
小弟執行了一個程式,找不出問題
原程式如下
syms x
syms k L P A I E x1 x2 real
Bk=k.*pi./L;
Xk=2^(1/2).*sin(Bk.*x);
for m=1:1:7
Bm=(m.*pi)./L;
Xm=2^(1/2).*sin(Bm.*x);
g(1,m)=Xm; ％編列1*m矩陣
end;
f1=Xk.*g.*(Heaviside(x-x1)-Heaviside(x-x2)); ％f1為x的函數矩陣 k為常數
Q1km=int(f1,0,L) ％積分後Q1km為x的函數矩陣
k為常數
symsum(Q1km,k,1,7) ％執行k=1～7的矩陣連加
這邊我是要找k=1~7帶入後疊加出來的答案
這邊的Q1km大概類似(k,2k,3k,4k,5k,6k,7k)
錯誤如下：
??? Error using ==> sym.maple at 87
interval of summation
Error in ==> sym.symsum at 43
r = maple('map','sum',f,[x.s '=' a.s '..' b.s]);
Error in ==> Untitled2 at 13
symsum(Q1km,k,2,3)
===================================================
我有檢查過我的分母帶入K不會有0出現 不過分子會有,
但是我用另外一組Q2km做就不會有問題
f2=diff(Xk).*g.*diff((Heaviside(x-x1)-Heaviside(x-x2)));
Q2km=int(f2,0,L)
symsum(Q2km,k,1,7) ===>這樣答案跑得出來
請問我是那邊弄錯了,這邊搞好久T.T 謝謝！！！！！
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.66.117
※ 編輯: inoran54787 來自: 140.115.66.117 (03/01 21:40)
→
03/02 17:20, , 1F
03/02 17:20, 1F