2011-05-13 29 views
0

私たちは簡単なODEを解こうとします。テキストbook (warning - russian)を使用します。これは、ブック(18ページ)にサンプルとして提示コードです:Maple ODE解決、グラフィックが完成していません。どうして?

> with(DEtools); with(plots); 
> eq1 := diff(x(t), t) = y(t); eq2 := diff(y(t), t) = -sin(x(t)); 

> init := [[x(0) = 1, y(0) = 0], [x(0) = 0, y(0) = 2]]; 
> DEplot([eq1, eq2], [x, y], t = 0 .. 110, x = -2*Pi .. 2*Pi, init, scene = [x, y], stepsize = 0.5e-1, obsrange = true, color = black, linecolour = black, arrows = 'MEDIUM') 

このコードは描くもの:enter image description here

しかし、楓は私のためにこれを出力しますenter image description here

1)それだけのための機能を解きなぜ四分の一?

2)本のように見えるようにするにはどうすれば解決できますか?

答えて

2

あなたが表示しているプロットに何か間違っているように見えることはありません。あなたが上に貼り付けたものに基づいて、本が与えるプロットはもっと驚いています。あなたの初期条件を考えてください。最初のものは[x(0)= 1、y(0)= 0]で、ポイント(1,0)から始まり、時間の経過とともに110ずつ増加します。これは明らかにインテリアサークルです。

[x(0)= 0、y(0)= 2]の初期条件の2番目のセットは、(0,2)の点のプロットの先頭次にそれはフィールドラインに沿って移動し、小さな内装プロットとほぼ同じ弧長になるように見えます。フィールドラインで判断すると、これは意味があります。 1つの問題は、ポイント(pi、0)の周りのカスプで数値に敏感である可能性があります。マシンのイプシロン内でy値が0までゼロにならない場合は、次に移動するポイントは、プロットが示すように右上がりに戻ります。

少し異なる初期条件と時間の増分で遊んでみましたか?ポイント(-pi、0)から始めて、もっと長い時間(例えばt = 0..300)行って、何が起こるかを見てみましょう。時にはこれらの本の中で、彼らが表示することを選択したプロットは、彼らがあなたに与えるコードから正確に生成されていません。

+1

良い考え。 [x(0)= 0、y(0)= 2-0.001]や[x(0)= 0、y(0)= 2-0.000001]などのICも考えてみましょう。 、x]または[t、y]を読み、そのICのわずかな変化に対して何が起こるかを見てください。 – acer

関連する問題