function z = flux(u)
     z  = .5*u.^2;