program trapesium

%program metode trapesium
disp('PROGRAM METODE TRAPESIUM')
disp('========================')
f=input('Masukkan fungsi f (x) simbolik: ');
a=input('Masukkan batas bawah           : ');
b=input('Masukkan batas atas            : ');
n=input('Masukkan cacah interval        : ');
fa=subs(f,'x',a);
fb=subs(f,'x',b);
Trap=fa+fb;h=(b-a)/n;

if n > 1
    for i=1:n-1
        nilf=subs(f,'x',a+i*h);
        Trap=Trap+2*nilf;
    end
end
Trap=(h/2)*Trap;
error=(-1/12)*(b-a)*h^2*subs(diff(f,2),'x',(b-a)/2);
hasil=sprintf('Nilai Integrasi Pendekatan   = %5.6g ',Trap);
nilerror=sprintf('Estimasi Error            = %5.6g ',error);
disp('    ')
disp(hasil)
disp(nilerror)