function dydt = example(t,y)
dydt = 4*exp(0.8*t)-0.5*y;