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