def METHOD_NAME(x, epsilon=1e-5):
    if x.dtype is not np.float64:
        x = x.astype(np.float64)
    axes = tuple(range(1,x.ndim))
    mean = np.mean(x, axis=axes, keepdims=True)
    var = np.var(x, ddof=0, axis=axes, keepdims=True)
    return (x - mean) / np.sqrt(var + epsilon)