import numpy as np

def null(a, rtol=1e-05):
    u, s, v = np.linalg.svd(a)
rank = (s > rtol * s[0]).sum()
return rank, v[rank:].T.copy()
простой, похожий на Matlab способ поиска нулевого пространства небольшой матрицы в numpy (и форматирование чисел)