function y = luas_pp(a,b,I)
dx = (b - a)/I; % Step-size delta x
Luas = 0;
for i = 1:I,
    x = a + (i-1)*dx;
    tinggi = absfungsi(x);
    Luas = Luas + (dx*tinggi);
end
y = Luas;