import pylab as pl
pl.plot(all_l, qtr)
pl.grid()
pl.show()
sys.exit(0)
if debug:
assert(l >= 0)
After Change
"to one.")
if debug:
assert alpha >= 0, alpha
y[:, m] = (alpha * c) / (1. + alpha * q) // x*
// These lines can be put out of this loop. Used only to
// compute criterion for debugging.