function dy=Rossler(t,y)
a=0.2;  % Rossler 方程參數 a
b=0.2;  % Rossler 方程參數 b
c=5.7;    % Rossler 方程參數 c
dy=zeros(3,1);
dy(1)=(-1)*(y(2)+y(3));
dy(2)=y(1)+a*y(2);
dy(3)=b+y(3)*(y(1)-c);