docs: def _build_matrix(p, q, deriv):
      """Constructs the equation system matrix for the finite difference coefficients"""
      A = [([1 for _ in range(-p, q+1)])]
      for i in range(1, p + q + 1):
          A.append([j**i for j in range(-p, q+1)])
      return np.array(A)