2014-09-23 96 views
0

假設Mathematica中此代碼:斜面情節衍生物

w[w1_] := NDSolve[{y''[x] + y[x] == 2, y'[0] == 0, y[0] == w1},y[x], {x, 0, 30}]; 
Manipulate[Plot[Evaluate[{y[x], y'[x]} /. w[w1]], {x, 0, 30}], {w1, -2, 3}] 

的代碼以這種方式工作:一個滑桿(在操縱)控制一個初始條件微分方程的保存([0] y的值)在變量w中,當它改變時,繪製解決方案。這工作正常,但我有一個解決方案的一階導數(y'[x])的問題。相反的函數圖形的存在只是簡單的空白,並沒有被繪製,即使我檢查了它在數學的幫助,這裏的示例代碼

s = NDSolve[{y''[x] + Sin[y[x]] y[x] == 0, y[0] == 1, y'[0] == 0}, y, {x, 0, 30}] 
Plot[Evaluate[{y[x], y'[x], y''[x]} /. s], {x, 0, 30},PlotStyle -> Automatic] 

作品。有人知道問題在哪裏嗎?

回答

1

我發現這個問題,它只是在NDSolve它是需要寫

w[w1_] := NDSolve[{y''[x] + y[x] == 2, y'[0] == 0, y[0] == w1},y, {x, 0, 30}] 

得到Y而不是Y [X]和一切正常。