factorial mat

function r=fact(n)
    r = 1;
    if n == 0
        return
    end
    for i=2:n
        r = r * i;
    end
end