clear
k=input('amortissement : ');
Nb = input('Nombre de simulations : ');

Simu_Phase_4
figure(2)
clf
axis([-10 10 -5 5])
hold on

for i =1 : Nb
[x0 ,y0]=ginput(1)
sim('Simu_Phase_4',[0 20])
figure(2)
t= plot(x,y,'r',x0,y0,'*g')

end

pause
figure(2)
clf
axis([-10 10 -5 5])
hold on
 x0 =-5;
 y0=3;
Nb = input('Nombre de simulations en k variant: ');

for i =1 : Nb
   k = 0+2*i/Nb;
   home
   disp(k)
   sim('Simu_Phase_4',[0 20])
	figure(2)
	t= plot(x,y,'r',x0,y0,'*g')
end