def METHOD_NAME(str_method, qin_arr, axes):
    eps = 1.0e-8
    arr = np.random.random((3, 4, 5))
    qin_arr = [0.1, 37.0, 42.0, 67.0, 83.0, 100.0, 49.0, 50.0, 0.0]
    axes = 0
    q_out = num.percentile(arr, qin_arr, method=str_method, axis=axes)
    np_q_out = np.percentile(arr, qin_arr, method=str_method, axis=axes)
    assert q_out.shape == np_q_out.shape
    assert q_out.dtype == np_q_out.dtype
    assert allclose(np_q_out, q_out, atol=eps)