function phi = minmod(a,b)
    phi = (a.*b>0).*( (abs(a)<abs(b)).*a + (abs(b)<=abs(a)).*b ) ;
end