% PeceRungeKuttaMain.m
%
% PECE法およびルンゲクッタ法のスクリプトファイル
% 呼び出される関数
% pece.m あるいは　rungekutta.m
% 
% Copyright (C) 2007 清水 信行, 曽我部 潔, 小金澤 鋼一 All Rights Reserved.

% =================================================================================
% プログラムの使用条件
% 1. 本プログラムおよびデータを商用で，頒布してはならない．
% 2. 本プログラムおよびデータを著作者の許可なく営利目的で利用してはならない．
% 3. 本プログラムおよびデータの改変は，営利目的でない限り自由である．
% 4. 本プログラムおよびデータを使用することにより生じた損害などについては，著者，
%    プロクラム著作者，(社)日本機械学会，および(株)コロナ社は一切の責任を負わない．
% 5. 本プログラムの著作権は，本書編著者ならびにプログラム著作者に帰す．
% =================================================================================

clear all;

% 微分方程式の定義
    eq={'xdot' '10*(1-x^2)*xdot-x'};
%　時間積分
    tic; sol=pece(eq, [1,0],20,0.01); toc;
% 結果の表示
    plot(sol(1,:), sol(2,:));
    xlabel('x');
    ylabel('xdot');
    grid;