Write program to output exact results for Kitaev model