function?dy=?ey1（ty）
%UNtitleD2?此處顯示有關此函數的摘要
%???此處顯示詳細說明
b=1;
w=0.3;
r=0.35;
a=0.2;
dy=zeros（41）;
dy（1）=-b*y（1）*y（3）+a*y（4）;
dy（2）=b*y（1）*y（3）-w*y（2）;
dy（3）=w*y（2）-r*y（3）;
dy（4）=r*y（3）-a*y（4）;


end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????239??2019-03-10?14:49??seirs\ey1.m

?????文件??????30048??2019-03-10?15:45??seirs\result.jpg

?????文件????????214??2019-03-10?14:50??seirs\seirs.m

?????目錄??????????0??2019-03-10?16:13??seirs

-----------?---------??----------?-----??----

????????????????30501????????????????????4