// Question 3

Sol = eulerE(f,[0;1],5,0.001)

// Question 4

scf(0)
clf()

plot2d([0:0.001:5],Sol(1,:))

scf(1)
clf()

plot2d([0:0.001:5],Sol(2,:))

// Question 5

X=5*sin(Sol(1,:))
Y=-5*cos(Sol(1,:))

scf(2)
clf()

plot2d(X,Y,rect=[-5,-6,5,0])
// suivant les versions de Scilab la commande isoview peut ne pas fonctionner de la même manière (cf documentation)
isoview

// Question 6

Sol2 = eulerE(f,[0;1.5],5,0.001)

X2=5*sin(Sol2(1,:))
Y2=-5*cos(Sol2(1,:))

scf(3)
clf()

plot2d(X2,Y2,rect=[-5,-6,5,0])
// suivant les versions de Scilab la commande isoview peut ne pas fonctionner de la même manière (cf documentation)
isoview

// Lorsque la vitesse de départ est plus élevée, on observe que le pendule a une trajectoire plus "ample" (il monte plus haut). Ceci est assez naturel d'un point de vue physique.
