def METHOD_NAME(self, spcmin, spcmax, spclen=1000):
    # Dimensions of the spectrum
    delta = ( spcmax - spcmin ) / ( spclen - 1 )
    spc = numpy.zeros((spclen, 2), dtype=float)
    for ix in range(spclen):
        spc[ix, 0] = spcmin + delta*ix
    return spc