
function dydt = fpendulo1(t,y)
dydt = [y(2); -sin(y(1))];