f1=inline(['[-8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3);',...
    '-x(1)*x(2)+28*x(2)-x(3)]'],'t','x');

t_final=100; x0=[0;0;1e-10]; [t,x]=ode45(f1,[0,t_final],x0);
plot(t,x), figure;
plot3(x(:,1),x(:,2),x(:,3)); axis([10 42 -20 20 -20 25]);