The partition function of two dimensional massless staggered fermions interacting with U(N) gauge fields is rewritten in terms of loop variables in the strong coupling limit. We use this representation of the theory to devise a non-local Metropolis algorithm to calculate the chiral susceptibility. For small lattices our algorithm reproduces exact results quite accurately. Applying this algorithm to large volumes yields rather surprising results. In particular we find m π ≠ 0 for all N and it increases with N. Since the talk was presented we have found reasons to believe that our algorithm breaks down for large volumes questioning the validity of our results.