A new numerical method for div-curl systems with low regularity assumptions. (15th May 2022)