Div(B) = 0
And update it to say
Div(B) = sigma
Where sigma is a field describing the monopole density. Theres a ready "gap" in Gauss' law for magnetism where you can easily stick monopoles. Of course the divergence would be zero in the absence of monopoles, just as the divergence of the electric field is zero in the absence of electric monopoles, but decidedly non-zero when there's an electron around.
