program multifactorial;
    loop for d in [1..5] do
        print(d, ":", [multifac(n, d) : n in [1..10]]);
    end loop;

    proc multifac(n, d);
        return */{n, (n-d)..1};
    end proc;
end program;