In this work, equations of state of hard, anisotropic particles are studied employing the complementary methods of virial theory and Monte Carlo simulations. To calculate higher-order virial coefficients, an algorithm is developed and optimised for hard particles. Using this algorithm, virial coefficients of hard lenses and hard, oblate ellipsoids of revolution are calculated. The method is extended to the 4D space. The development of a cluster Monte Carlo algorithm gives access to precise equation-of-state data, which is compared to the virial series.