clear;
x=[]; y=[]; z=[];
x(1)=3;
y(1)=15;
z(1)=1;
for i=1:10000
    x(i+1)=10*0.01*y(i)-x(i)*(10*0.01-1);
    y(i+1)=0.01*x(i)*(28-z(i))-y(i)*(0.01-1);
    z(i+1)=0.01*x(i)*y(i)-z(i)*(0.01*8/3-1);
end
figure(1)
%subplot(222); % Time»~®t­È0.01
plot3(x,y,z,'gx');
hold on 
%xlabel('-X-');
%ylabel('-Y-');
%zlabel('-Z-');
view(62,5)
set(gca,'color','k')
set(gca,'xcolor','b','ycolor','b','zcolor','b')
%grid on;


%@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
x(1)=3;
y(1)=15;
z(1)=1.000001;
for i=1:10000
    x(i+1)=10*0.01*y(i)-x(i)*(10*0.01-1);
    y(i+1)=0.01*x(i)*(28-z(i))-y(i)*(0.01-1);
    z(i+1)=0.01*x(i)*y(i)-z(i)*(0.01*8/3-1);
end
%figure(2)
%subplot(222); % Time»~®t­È0.01%
plot3(x,y,z,'wx');
xlabel('-X-');
ylabel('-Y-');
zlabel('-Z-');
grid on;