A=[-8/3, 0, 0;0, -10, 10;0, 28, -1];
y=[35 -10 -7]';
h=0.01;
p=plot3(y(1),y(2),y(3),'.','EraseMode','none','MarkerSize',5);
axis([0 50 -25 25 -25 25]);
hold on;
for i=1:4000
   A(1,3)=y(2);
   A(3,1)=-y(2);
   ydot=A*y;
   y=y+h*ydot;
   set(p,'XData',y(1),'YData',y(2),'ZData',y(3));
   drawnow;
end