function fac = getFacVec(n)
    fac = prod(1:n);
end