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