function dydt=odefun2(t,y)

v=y(2);

yprime=v;
vprime=-v-y(1);

dydt=[yprime; vprime];