function y = hw31fun(x)
%
% Function for HW 31

y = x.^2.*exp(-x.^2) + 0.5*sin(3*x);