function dy=first_biased(t,y,cdx)
dy=-cdx*(y-circshift(y,1));