with(LinearAlgebra):
with(VectorCalculus):
BasisFormat(false):
f:=x1+x2*exp(-x3*t)-y;
Gradient(f,[x1,x2,x3]);
Hessian(f,[x1,x2,x3]);