import numpy as np
from sklearn.metrics import median_absolute_error as skmedian_absolute_error

def median_absolute_error(y_true, y_pred):
    return np.median(np.abs(y_true - y_pred))

for i in range(10):
    rng = np.random.RandomState(i)
    y_true = rng.rand(10)
    y_pred = rng.rand(10)
    score1 = median_absolute_error(y_true, y_pred)
    score2 = skmedian_absolute_error(y_true, y_pred)
    assert np.isclose(score1, score2)