%dde_example.m
%解延遲微分方程的例子
lags = [1 3];                                     %延遲向量
sol = dde23(@ddefun,lags,@ddefun_history,[0,1]);  %解方程
hold on;
plot(sol.x,sol.y(1,:),'b-');                        %畫出結果
plot(sol.x,sol.y(2,:),'r-.');
title('延遲微分方程的解');                            
xlabel('t');
ylabel('y');
legend('y_1','y_2',2);