s=0:.01:1; X=6*s.^2-5*s.^3; Y=6*s-12*s.^2+6*s.^3; plot(X,Y) hold all
set(gca,'XLim',[-0.2 1.5],'YLim',[-0.2 1.5], 'Drawmode', 'fast', 'visible', 'on'); axis square ball=line('color','r','Marker','o','MarkerSize',10,'LineWidth',2,'erase','xor','xdata',[],'ydata',[]); ball2=line('color','b','Marker','o','MarkerSize',10,'LineWidth',2,'erase','xor','xdata',[],'ydata',[]);
t=.00:.05:20;
t2 = zeros(20,1);
for i=0:20 t2(i+1)=newtonhjalp(i*.05,.001); end
a=6*t.^2-5*t.^3; b=6*t-12*t.^2+6*t.^3;
c=6*t2.^2-5*t2.^3; d=6*t2-12*t2.^2+6*t2.^3;
for i = 1:21 set(ball,'xdata',a(i),'ydata',b(i)); set(ball2,'xdata',c(i),'ydata',d(i)); drawnow; pause(0.01); frame = getframe(1);
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
end