def METHOD_NAME(l, alpha, beta, *, normalization: str = "integral"):
    r"""Spherical harmonics of :math:`\vec r = R_y(\alpha) R_x(\beta) e_y`
    .. math:: Y^l(\alpha, \beta) = S^l(\alpha) P^l(\cos(\beta))
    where :math:`P^l` are the `Legendre` polynomials
    Parameters
    ----------
    l : int or list of int
        degree of the spherical harmonics.
    alpha : `torch.Tensor`
        tensor of shape ``(...)``.
    beta : `torch.Tensor`
        tensor of shape ``(...)``.
    Returns
    -------
    `torch.Tensor`
        a tensor of shape ``(..., 2l+1)``
    """
    sh = SphericalHarmonicsAlphaBeta(l, normalization=normalization)
    return sh(alpha, beta)