  We present a method for direct hybrid Monte Carlo simulation of graphene on
the hexagonal lattice. We compare the results of the simulation with exact
results for a unit hexagonal cell system, where the Hamiltonian can be solved
analytically.
