def METHOD_NAME(self):
    x = np.array([1])
    y = np.array([2])
    I0_sersic = 1
    R_sersic = 1
    n_sersic = 1
    center_x = 0
    center_y = 0
    values = self.sersic.function(
        x, y, I0_sersic, R_sersic, n_sersic, center_x, center_y
    )
    npt.assert_almost_equal(values[0], 0.12658651833626802, decimal=6)
    x = np.array([0])
    y = np.array([0])
    values = self.sersic.function(
        x, y, I0_sersic, R_sersic, n_sersic, center_x, center_y
    )
    npt.assert_almost_equal(values[0], 5.1482559148107292, decimal=2)
    x = np.array([2, 3, 4])
    y = np.array([1, 1, 1])
    values = self.sersic.function(
        x, y, I0_sersic, R_sersic, n_sersic, center_x, center_y
    )
    npt.assert_almost_equal(values[0], 0.12658651833626802, decimal=6)
    npt.assert_almost_equal(values[1], 0.026902273598180083, decimal=6)
    npt.assert_almost_equal(values[2], 0.0053957432862338055, decimal=6)
    value = self.sersic.function(
        1000, 0, I0_sersic, R_sersic, n_sersic, center_x, center_y
    )
    npt.assert_almost_equal(value, 0, decimal=8)