function dy=second_centered(t,y,cdx)
dy=-(cdx/2)*(circshift(y,-1)-circshift(y,1));