function [dydt]=odefun1(t,y)

dydt=-t*y;