% define f in y'=f(t,y)
function ydot=f_is(t,y)
ydot=-10*y;