Euclidean iterative method to find solutions of ax + by = GCD(a,b)