function y = der_fun_2(x)
y = (2*x^2 -1)*exp(-x^2);