single_trapezoid_rule_integration

%single_trapezoid_rule_integration
syms x
fx =x^2+x+2;
a = 0;
b = 2;
fy = int(fx,x);
fya = subs(fy,x,a);
fyb = subs(fy,x,b);
hs = fyb-fya
fxa = subs(fx,x,a);
fxb = subs(fx,x,b);
hk = (b-a)*(fxa+fxb)/2
ea= abs((hs-hk)/hs)