Re: [問題] deqn 問題再哪裡呢??

[chungyuandye (養花種魚數月亮賞星星)]

※ 引述《ntust661 (TOEFL_5!)》之銘言：
: http://ppt.cc/ipoL
: 我都刪除變數了
: 他一直說我的方程式有錯@@
: why??
: 各位大大救救小弟阿QQ

{Pa,Qa,Maa}={1,1,1};
sys={\[Theta]'[s]==4Pi^2*(Maa-Pa-Qa*X[s]),
X'[s]==Cos[\[Theta][s]],
y'[s]==Sin[\[Theta][s]],
M'[s]==Fx[s]*Sin[\[Theta][s]]-Fy[s]*Cos[\[Theta][s]],
Fx'[s]==0,
Fy'[s]==0};
BC={\[Theta][0]==0,X[0]==0,y[0]==0,Fx[0]==-Pa,
Fy[0]==Qa,M[0]==Maa};
ans=NDSolve[
Flatten@{sys,BC},{\[Theta][s],X[s],y[s],M[s],Fx[s],
Fy[s]},{s,0,1}][[1]]

GraphicsGrid@Partition[Plot[ans[[#,2]],{s,0,1}]&/@Range[6],3]

"=" 是 Set 不是 Equal

--
養花種魚數月亮賞星星

http://chungyuandye.twbbs.org


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 218.173.129.126
推 ntust661:老師好._./ 我不知道我錯在哪裡耶QQ...
→ chungyuandye:一階導數那些要用==，BC那邊最後應該是M[0]==Ma
推 ntust661:我使用你的打法 他最後出現
→ ntust661:{0==4*pi^2*x[s],True,True, ...}
→ ntust661:除了第一個不是True,其他都True!
→ ntust661:好多True,可是沒半個interpolation function QQ
推 ntust661:重開Kernel true 就沒出來了
→ ntust661:然後發現一個重要的東東!!!
推 ntust661:IC = {x'[t],y'[t],z'[t],...} 直接放進去NDSolve
→ ntust661:他竟然不會讀到 X'[t] == ...