def METHOD_NAME(
    displacement: Union[float, np.ndarray],
    c1: float = 0.335,
    c2: float = -0.048,
    c3: float = 1.565,
) -> Union[float, np.ndarray]:
    """
    Computes the probability of ground failure using a Weibull
    model based on the predicted Newmark displacements.
    Constants from Jibson et al., 2000 for the Northridge earthquake.
    :param displacement:
        Predicted Newmark displacement at a site, in meters (will be converted
        to cm in the function).  Can be scalar or array.
    :returns:
        Scalar or array of ground failure probability.
    """
    Dn = displacement * 100.0
    return c1 * (1 - np.exp(c2 * Dn ** c3))