function f=te1(t,x)
f=3.*x-x.^2+sin(t);