Maple ODE 해결 중, 그래픽이 완전하지 않습니다. 왜?

Question

그래서 우리는 간단한 ODE를 풀려고합니다. 우리는 텍스트 book (warning - russian)을 사용합니다. 책에서 샘플로 제시 이것은 코드 (18 페이지) :

을하지만 나를 위해 메이플 출력이 :

> 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') 

이 코드 작성하여야

1) 왜에만 기능을 해결 사분의 일?

2) 책 에서처럼 보이게 만드는 방법?

Answer 1

내가 보여주는 플롯에 문제가있는 것처럼 보이지 않습니다. 당신이 위에 붙여 넣은 것에 따라 저에게 책이주는 줄거리가 더 놀랍습니다. 초기 조건을 생각해보십시오. 첫 번째 것은 [x (0) = 1, y (0) = 0]이며, 점 (1,0)에서 시작하여 110 번씩 앞으로 이동합니다 (제공된 코드에 따라). 이것은 분명히 내부 동그라미입니다.

[x (0) = 0, y (0) = 2]의 초기 조건의 두 번째 세트는 올바른 위치 (지점 (0,2)에서 플롯의 상단) 그 다음에 그것은 작은 줄거리와 거의 동일한 원호 길이로 보이기 위해 줄을 따라 이동합니다. 필드 라인으로 판단 할 때 이것은 더 합리적입니다. 하나의 문제는 점 주위의 첨두 (pi, 0)에서 수치 적으로 민감 할 수 있습니다. 기계 엡실론 내에서 y 값을 0으로 줄이지 않으면 다음으로 이동할 지점이 플롯이 표시하는대로 오른쪽 상단으로 등반합니다.

약간 다른 초기 조건 및 시간 단위로 재생 해 보았습니까? 포인트 (-pi, 0)부터 시작하여 훨씬 오래, 말하자면 t = 0..300으로 진행하고 어떤 일이 일어나는지보십시오. 때로는 이와 같은 책에서 그들이 보여 주기로 한 음모가 자신이 제공 한 코드에서 정확하게 생성되지는 못했지만이 점에 대해서는 언급하지 않았습니다.