xaxis = -1:0.001:2;
yaxis = exp(xaxis);

yaxisapprox = (2-4*sqrt(exp(1))+2*exp(1))*xaxis.^2 + (-3+4*sqrt(exp(1))-exp(1))*xaxis+1;


clf
plot(xaxis,yaxis,'r')
hold on
plot(xaxis,yaxisapprox,'b')
xlim([-1 2])
ylim([0 3])
grid on