lowerPascal[size_] := MatrixExp[ SparseArray[{Band[{2, 1}] -> Range[size - 1]}, {size, size}]]]