Discretization errors at free boundaries of the Grad-Schlüter-Shafranov equation