clear

t = 50;
t1=t;
dt=0.1;



Y(1) = 0;
Y(2) = 8;
options = odeset('RelTol',1e-7,'OutputFcn','odeplot','OutputSel',[1]);
[t Y] = ode45('oscillatorani1', [0:dt:t], Y, options);


for i=1:t1/dt, xx(1,i)=0; yy(1,i)=0; xx(2,i)=sin(Y(i,1)); yy(2,i)=-cos(Y(i,1)); 
end

plot(xx(:,100),yy(:,100),'ro-');
axis([-1.1 1.1 -1.1 1.1])
% line(xx(:,501),yy(:,501));

n=1;
fig=figure;
set(fig,'DoubleBuffer','on');
set(gca,'NextPlot','replace','Visible','off');
iname=sprintf('mov0.avi', 1)
mov = avifile(iname,'fps',5);
for count=1:400
    h = plot(xx(:,count),yy(:,count),'ro-');
    shading flat;
    axis([-1.1 1.1 -1.1 1.1])
grid off
    set(h,'EraseMode','xor');
    F = getframe(gca);
    mov = addframe(mov,F);
    clear F;
    clear h;
end
mov = close(mov);
 clear mov;