clear;
clc;

f = @(x) (exp(x));
C = 1.6;
W = 3.2;

n = 10^4;
h = (W - C) / n;

result = f(W) + f(C);
for i = (C + h):h:(W - h)
    result = result + 2*f(i);
end
result = h/2 * result;

disp("Result: " + result);