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