function f = f18_4_1(t,x)
f = -exp(10*t)*x*x;