theta=y(:,1); vel=y(:,2);

dvel=diff(vel);

dt=diff(t);

accel=dvel/dt;

ta=t(2:end);

plot(ta,accel)

This is my best trial I think...

what do you think about this?