u = np.linspace(start, stop, num)  # Discretization levels