function dydt=linkde(t,y)
dydt=-2*y;