def casoratian(seqs, n, zero=True):
    from .dense import Matrix

    seqs = list(map(sympify, seqs))

    if not zero:
        f = lambda i, j: seqs[j].subs(n, n + i)
    else:
        f = lambda i, j: seqs[j].subs(n, i)

    k = len(seqs)

    return Matrix(k, k, f).det()
Given linear difference operator L of order 'k' and homogeneous
dataset/ETHPy150Open sympy/sympy/sympy/matrices/dense.py/casoratian