function dydt=zh21de(t,y)
dydt=-3*y;