The code is developed for the calculation of the magneto-crystalline anisotropy (MAE) in thin films using a classical Heisenberg hamiltonian with a correction developed by Van Vleck. A Metropolis style Monte Carlo algorithm was used with adequate corrections to accelerate the calculation. The MAE was calculated for the case of...