a=0; b=2; k=16; 
I = trapzc(a,b,k,'f1')

e= I - (exp(b)-1),