clear
Simu_Phase_5
figure(1)
clf
axis([-0.5 4 -0.5 4])
plot([1 1],[-0.5 4],':g');
hold on
plot([-0.5 4],[1 1],':g');
plot([-0.5 4],[0 0],':r');
plot([0 0],[-0.5 4],':r');

axis([-0.5 4 -0.5 4])

Nb = 10
for i =1 : Nb
[x0 ,y0]=ginput(1)
sim('Simu_Phase_5',[0 10])
t= plot(x,y,'b',x0,y0,'*g')

end