clear; % Menghapus data sebelomnya
display('CONTOH PROGRAM MENACRI LUAS BIDANG')
toleransi = input('masukkan nilai toleransi (dalam %): ')
N = input('masukkan jumlah INTERVAL maximum N = ')
a = input('masukkan batas bawah = ')
b = input('masukkan batas atas = ')
if a > b % Kalo' batas bawah lebih besar dari batas atas
    c = b; % Tukarkan
    b = a; % nilai a dan
    a = c; % nilai b
    for I = 1:N,
        INT(I)=I;
        Lpp(I)= luas_pp(a,b,I); % Metode Persegi Panjang
        Ltr(I)= luas_tr(a,b,I); % Metode Trapesium
        Error(I) = (abs(Lpp(I)-Ltr(I))/(Ltr(I)))*100; % Error dalam %
        if Error(I) <= toleransi
            break
        end
    end
    if length(Error) == N
        display('Batas toleransi tidak dicapai, Error masih besar')
    end
    subplot(311), plot(INT,Lpp), grid on, xlabel('Jumlah Interval'), ylabel('Luas Persegi Panjang')
    subplot(312), plot(INT,Ltr), grid on, xlabel('Jumlah Interval'), ylabel('Luas Trapesium')
    subplot(313), plot(INT,Error), grid on, xlabel('Jumlah Interval'), ylabel('Error (%)')
    format short
    display('Hasil Komputasi Luas')
    display('Jumlah Interval    Luas Persegi Panjang    Luas Trapesium  Error')
    [INT'   Lpp'    Ltr'    Error']
elseif a == b
    display('Kok batas atas sama dengan batas bawah ???')
    display('Luasnya pasti NOL dong')
else
    for I = 1:N,
        INT(I)=I;
        Lpp(I)= luas_pp(a,b,I); % Metode Persegi Panjang
        Ltr(I)= luas_tr(a,b,I); % Metode Trapesium
        Error(I) = (abs(Lpp(I)-Ltr(I))/(Ltr(I)))*100; % Error dalam %
        if Error(I) <= toleransi
            break
        end
    end
if length(Error) == N
    display('Batas toleransi tidak dicapai, Error masih besar')
end
subplot(311), plot(INT,Lpp), grid on, xlabel('Jumlah Interval'), ylabel('Luas Persegi Panjang')
subplot(312), plot(INT,Ltr), grid on, xlabel('Jumlah Interval'), ylabel('Luas Trapesium')
subplot(313), plot(INT,Error), grid on, xlabel('Jumlah Interval'), ylabel('Error (%)')
format short
display('Hasil Komputasi Luas')
display('Jumlah Interval    Luas Persegi Panjang    Luas Trapesium  Error')
[INT'   Lpp'    Ltr'    Error']
end