def METHOD_NAME(m, u, w, r):
    '''Evaluates the radial profile of the LP modes.
    Parameters
    ----------
    m : int
        The azimuthal order
    u : scalar
        The normalized inner propagation constant.
    w : scalar
        The normalized outer propagation constant.
    r : array_like
        The radial coordinates on which to evaluate the bessel modes.
    Returns
    -------
    array_like
        An array that contains the radial profile.
    '''
    # The scaling factor for the continuity condition
    scaling_factor = jv(m, u) / kn(m, w)
    # Find the grid inside and outside the core radius
    mask = r < 1
    # Evaluate the radial mode profile
    mode_field = np.zeros_like(r)
    mode_field[mask] = jv(m, u * r[mask])
    mode_field[~mask] = scaling_factor * kn(m, w * r[~mask])
    return mode_field