An accurate scalar potential finite element method for linear, two‐dimensional magnetostatics problems