import numpy as np from pcw import dls n = 5 alpha = np.logspace(-1, -5, n) x = 0.5/alpha uw = np.logspace(-9, 1) t = 0.25/uw for i in range(len(x)): s, _