global k l m M g A w
k=445; l=0.5; m=3; M=8; g=9.81; A=0; w=2;
x10=0.3; x_dot10=0; teta10=pi/4; teta_dot10=0;
x20=0.30001; x_dot20=0; teta20=pi/4; teta_dot20=0;
options=odeset('RelTol',1e-4);
frp = sqrt(l/g); %facteur de la frequence propre du pendule aux petits angles
p = pi*frp;  %periode d'echantillonage
[t1,teta1] = ode45('foncPendule',[0:p:300*p],[x10 x_dot10 teta10 teta_dot10],options);
[t2,teta2] = ode45('foncPendule',[0:p:300*p],[x20 x_dot20 teta20 teta_dot20],options);
box on; hold on
%axis([-3.5 3.5 -5 5]);

for i=1:1:size(t1)
    line(teta1(i,3), teta1(i,4), 'Marker', '.', 'Markersize', 5)
    line(teta2(i,3), teta2(i,4), 'Marker', '.', 'Markersize', 5, 'Color', 'r')
    drawnow
end

title('Section de Poincare stroboscopique  - echantillonnage: \pi*sqrt(l/g)')
xlabel('\theta')
ylabel('d\theta/dt')