function main()
    xmin = input('xmin = ');
    xmax = input('xmax = ');

    fun = @(x) 5 * (exp(1) ^ -((x - 3)^2));
    area = riemann(fun, xmin, xmax);

    area
end

function y = riemann(f, xmin, xmax)
    delta = 0.0001;

    y = 0;

    for t=xmin:delta:xmax
        y = y + delta * f(t);
    end
end