def METHOD_NAME(x, y, theta_E, gamma, e1, e2, center_x=0.0, center_y=0.0):
    """
    :param x: x-coordinate (angle)
    :param y: y-coordinate (angle)
    :param theta_E: Einstein radius (angle), pay attention to specific definition!
    :param gamma: logarithmic slope of the power-law profile. gamma=2 corresponds to isothermal
    :param e1: eccentricity component
    :param e2: eccentricity component
    :param center_x: x-position of lens center
    :param center_y: y-position of lens center
    :return: lensing potential
    """
    z, b, t, q, ang = param_transform(
        x, y, theta_E, gamma, e1, e2, center_x, center_y
    )
    alph = alpha(z.real, z.imag, b, q, t)
    return 1 / (2 - t) * (z.real * alph.real + z.imag * alph.imag)