h=15;
z=0:0.01:(2+h);
g=@(x) -0.00002*x.*exp(x).^(2*x)+0.25*x.^2+0.5*x+3/8;
plot(z,g(z)), hold on