def METHOD_NAME(const, dl, ds, St):
    dsCM = ds / const.CM
    dlCM = dl / const.CM
    Dss2 = (
        0.254 * (dsCM ** (0.413)) * np.exp(3.53 * dsCM ** (2.51) * (dlCM - dsCM))
    )  # (4.27)
    bstar = 14.2 * np.exp(-17.2 * dsCM)
    Ps20 = 0.23 * dsCM ** (-3.93) * dlCM ** (bstar)  # (4.29)
    sigmas2 = 10 * Dss2  # as in (4.22)
    mus2 = np.log(Dss2) + sigmas2**2  # (4.32)
    Hs2 = Ps20 * Dss2 / np.exp(-0.5 * sigmas2**2)  # (4.28)
    Fs = 5 * math.erf((St - 2.52e-6) / (1.85e-6)) + 6  # (3.7)
    for _ in range(10):
        sigmas2 = (
            np.sqrt(2 / np.pi)
            * (Fs - 1)
            / Hs2
            / (1 - math.erf((np.log(0.01) - mus2) / np.sqrt(2) / sigmas2))
        )
        mus2 = np.log(Dss2) + sigmas2**2  # (4.32)
        Hs2 = Ps20 * Dss2 / np.exp(-0.5 * sigmas2**2)  # (4.28)
    return (Hs2, mus2, sigmas2)