function I = Mtrap(a,b,n,fx)

    sum = fx(1,1);
    for i = 1:n-1
        sum = sum+2*fx(i,1);
    end
    sum = sum+fx(n,1);
    I = ((b-a)/2)*sum;
end