function z = g(x)
z = x*exp(-x);