def METHOD_NAME(v: np.ndarray) -> np.ndarray:
    """Find a unitvector in the orthogonal plane to v.
    Parameters
    ----------
    v : vector with shape (3,)
         axis ([x, y ,z])
    Returns
    -------
    vector which is either crossproduct with [0,1,0] or [1,0,0].
    Examples
    --------
    two vectors each of shape (3,)
    >>> find_random_orthogonal(np.array([1, 0, 0]))
    array([0, 0, 1])
    """
    if is_almost_parallel_or_antiparallel(v, np.array([1.0, 0, 0])):
        result = np.cross(v, [0, 1, 0])
    else:
        result = np.cross(v, [1, 0, 0])
    return normalize(result)