function y=fmif(x)
y=(3*x+x.^2).*exp(x);