
function demo_exp2

x=linspace(0,0.5,1000);
y=myexp(x);
sol=2.^x;
ae=max(sol-y)






function y=myexp(x)

y=1-0.15639e-16+(.69314718055995154416+(.24022650695869054994+...
(0.55504108675285271942e-1+(0.96181289721472527028e-2+...
(0.13333568212100120656e-2+(0.15403075872034576440e-3+...
(0.15265399313676342402e-4+(0.13003468358428470167e-5+...
0.12113766044841794408e-6*x).*x).*x).*x).*x).*x).*x).*x).*x;