def METHOD_NAME(n, x, ell):
    supportX1 = numpy.ndarray(shape=(n, ))
    supportX2 = numpy.ndarray(shape=(n, ))
    support_consts = numpy.ndarray(shape=(n, ))
    xScaled = float(2*x-n)/float(n)
    xCos = xScaled * ell.cosAlpha
    xSin = -xScaled * ell.sinAlpha
    supportX1.fill(xCos)
    supportX2.fill(xSin)
    supportY1 = numpy.arange(n)
    support_consts.fill(2.)
    supportY1 = numpy.multiply(support_consts, supportY1)
    support_consts.fill(n)
    supportY1 = numpy.subtract(supportY1, support_consts)
    support_consts.fill(n)
    supportY1 = numpy.divide(supportY1, support_consts)
    supportY2 = numpy.array(supportY1)
    support_consts.fill(ell.sinAlpha)
    supportY1 = numpy.add(supportX1,
                          numpy.multiply(supportY1, support_consts))
    support_consts.fill(ell.cosAlpha)
    supportY2 = numpy.add(supportX2,
                          numpy.multiply(supportY2, support_consts))
    support_consts.fill(ell.x0)
    supportY1 = numpy.subtract(supportY1, support_consts)
    support_consts.fill(ell.y0)
    supportY2 = numpy.subtract(supportY2, support_consts)
    support_consts.fill(ell.a)
    supportY1 = numpy.power((numpy.divide(supportY1, support_consts)),
                            2)
    support_consts.fill(ell.b)
    supportY2 = numpy.power(numpy.divide(supportY2, support_consts),
                            2)
    return numpy.add(supportY1, supportY2)