We study the flow of supercurrents between two superconducting contacts connected by a 2d layer of graphene. We use a Markov chain Monte Carlo method to find Andreev bound states for circular electron trajectories. Using sample trajectories we estimate the currents as function of the superconducting phase difference between the contacts and the magnetic field.