% ファンデルポール方程式の位相図を描く
interval = 0.00:0.01:10.00;
qinit = [2.00;0.00];
[time,q] = ode45(@van_der_Pol,interval,qinit);

% line style   solid -   broken -.   chain --   dotted :

set(0,'defaultAxesFontSize',16);
set(0,'defaultTextFontSize',16);

% t-x グラフをファイルに保存
plot(time,q(:,1),'-');
saveas(gcf, 'van_der_Pol_position.png', 'png');
fprintf("t-x グラフ\n");
%fprintf("一時停止：何かのキーを押してください\n");
%pause;

% t-v グラフをファイルに保存
plot(time,q(:,2),'-');
saveas(gcf, 'van_der_Pol_velocity.png', 'png');
fprintf("t-v グラフ\n");
%fprintf("一時停止：何かのキーを押してください\n");
%pause;

% x-v グラフ（位相図）をファイルに保存
plot(q(:,1),q(:,2),'-');
saveas(gcf, 'van_der_Pol_phase.png', 'png');
fprintf("x-v グラフ\n");
%fprintf("一時停止：何かのキーを押してください\n");
%pause;
