PlotData(vstack([di_low * ones(nu), di_up * ones(nu)]).T, u,