function dY = oscillatorani1(t,Y)

c = 0.1;
g = 10;
L = 1;




dY(1) = Y(2);
dY(2) = -(g/L)*sin(Y(1)) - c*Y(2);

dY = dY';