d/dx(x^2*exp(−x*(y+c)))