a = input('Masukkan nilai a sembarang =')
b = input('Masukkan nilai b sembarang =')
TOL = input('Masukkan nilai toleransi (dalam %) =')
N = input('Masukkan maksimum iterasi =')
Luas_Analitik = input('Masukkan Luas Analitik = ')
% Awal iterasi
for i = 1:N
    % Merekam iterasi ke i:
        k(i)=i-1; % Ini iterasi ke k
        N_interval(i) = 5*i; % jumlah interval 5*iterasi
        delta_x(i)= (b-a)/N_interval(i); 
    % Perhitungan luas empat persegi panjang
        Luas_4pp(i) = 0;
        for j = 1:N_interval(i)
            Luas_4pp(i) = Luas_4pp(i) + abs((delta_x(i)*fungsi(a + (j-1)*delta_x(i))));
        end
    % Perhitungan ERROR
        ERROR = (abs((Luas_Analitik - Luas_4pp(i))/Luas_Analitik))*100; % Error dalam persen
    if ERROR <= TOL
        'Luas Numerik diperoleh pada:' 
        ITERASI_KE = k(i)
        break
    end
end
plot(N_interval,Luas_4pp), grid on, xlabel('Jumlah Interval'), ylabel ('Luas_4pp')
ERROR_TERAKHIR =ERROR