function y = f1(t,x)

y = -x.*t.**2;

end