exec('boussole.sci');
a1=0.5;
a2=0.2;
tmax=500;
initial=[0.1,0;0.5,0;1,0;2,0];
                
plotA=scf();
poincare(a1,a2,initial,tmax,[-3,-3,3,3])
plotB=scf();
spectre(a1,a2,0.5,0,tmax,0.05,-100,100)
