Virtual Element approximation of 2D magnetostatic problems