multi_simpson

%multi_simpson
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);
ha = fyb-fya
fa = subs(fx,x,a);
fb = subs(fx,x,b);
n = 4;
h = (b-a)/n;
ganjil=0;
genap = 0;
for i=1:(n-1)l
if rem (i,2)== 0
genap=genap+subs(fx,x,a+i*h)
else
ganjil=ganjil+subs(fx,x,a+i*h) 
end
end
hk = h/3*(fxa+fxb+4*ganjil+2*genap)
ea = abs((ha-hk)/hs)