function phi = superbee(a)
    phi = max(0,max(min(1,2*a),min(2,a)));
end