your version of NumPy does not support the rtol keyword in the np.linalg.matrix_rank so you can upgrade it or use this code
import numpy as np
A = np.random.randn(4, 4)
tol = 1e-5
u, s, vh = np.linalg.svd(A)
rank = np.sum(s > tol * np.max(s))
print("Rank:", rank)