Simu_Phase_1   % ouverture du fichier simulink
figure(1)
clf
axis([-.01 .01 -.01 .01]) % fixe  le jeu d'axe autour de l'origine

for i = 1: 10
    [x0,y0]=ginput(1)   % saisie de la condition inititale sur le fenetre graphique
    sim('Simu_Phase_1',[0 10]) % simulation de système d'équations 
    plot(x,y) % tracé de la trajectoire de phase pour la CI choisie
    hold on   % mise en "surimpression" des différents tracés successifs
    axis([-.01 .01 -.01 .01]) % ramène le jeu d'axe autour de l'origine
    drawnow
end

axis([-2 2 -5 5])  % fixe  le jeu d'axe sur un plan de phase élargi
   
for i = 1: 50
    [x0,y0]=ginput(1)
    sim('Simu_Phase_1',[0 10])
    plot(x,y)
    hold on
    drawnow
end