function y = hw3der(x)
% 
% 4th derivative of function for HW 31

a = (4*x.^6 - 28*x.^4 + 39*x.^2-6);
y = 4*exp(-x.^2).*a + 81/2*sin(3*x);