def METHOD_NAME(self, x, y, kappa_0, theta_c, center_x=0, center_y=0):
    """
    :param x: angular position (normally in units of arc seconds)
    :param y: angular position (normally in units of arc seconds)
    :param kappa_0: central convergence of profile
    :param theta_c: core radius (in arcsec)
    :param center_x: center of halo (in angular units)
    :param center_y: center of halo (in angular units)
    :return: lensing potential (in arcsec^2)
    """
    x_ = x - center_x
    y_ = y - center_y
    r = np.sqrt(x_**2 + y_**2)
    r = np.maximum(r, self._s)
    Integral_factor = 0.5 * exp1((r / theta_c) ** 2) + np.log((r / theta_c))
    METHOD_NAME = kappa_0 * theta_c**2 * Integral_factor
    return METHOD_NAME