xvec = np.arange(250.) / 100.