def METHOD_NAME(im, sigma, xc, yc):
    "Creating 2D gaussian to be put in a test image"
    e = 1
    angle = 0
    sx = sigma * (1 + e)
    sy = sigma * (1 - e)
    size = int(8 * sigma + 1)
    if size % 2 == 0 :
        size += 1
    x = numpy.arange(0, size, 1, float)
    y = x[:, numpy.newaxis]