function y = der2_fun_2(x)
y = (6*x-4*x^3)*exp(-x^2);