% Y. X. Li, G. R. Chen, K. S. Wallace, IEEE Trans. on Cir. Sys.��52 (2005) 204
% Author: Thomas Lee
% E-mail: lixf1979@126.com
% Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
function xdot=lorrenz(t,x)
for t=0:0.1:100;
if(cos(5.3*t)>=0)
p=1;
else
p=-1;
end
end
xdot=[-8/3*x(1)+p*x(2)*x(3);
(-x(2)+x(3))*(25-10*cos(5.3*t));
-p*x(1)*x(2)-(17.5*cos(5.3*t)+10.5)*x(2)+(13.3-14*cos(5.3*t))*x(3)];
tic;
x0=[0;0;1e-10];tspan=[0,100];
[t,x]=ode45(@lorrenz,tspan,x0);
plot(t,x(:,1),'k');
r=length(t)-5000:length(t);
figure; plot(x(r,2),x(r,3),'k')
toc;
