Two-dimensional equilibrium in coronal magnetostatic flux tubes: an accurate equilibrium solver