A=[-8/3 0 0;0 -10 10;0 28 -1];
y=[10 15 -11]';
h=.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:5000
   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))
   pause(.1)
   drawnow
   i=i+1;
end
hold off