global k l m M g A w
k=170; 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;
options=odeset('RelTol',1e-7);
[t1,teta1] = ode45('foncPendule',[0:0.01:50],[x10 x_dot10 teta10 teta_dot10],options);

for i=1:size(t1)
    T(i) = 1/2 * teta1(i,2)^2 * (M+m) + m * l * teta1(i,2) * teta1(i,4) * cos(teta1(i,3)) + 1/2 * m * l^2 * teta1(i,4)^2;
    V(i) = 1/2 * k *(teta1(i,1) - A*sin(w*t1(i)))^2 - m * g * l * cos(teta1(i,3));
end
hold on
plot(t1,V)
plot(t1,T,'Color','r')
plot(t1,V+T,'Color','g', 'LineWidth', 2)

title('Energies potentielle, cinétique et totale du système')
xlabel('temps [s]')
ylabel('Energie [J]')