The configurational-bias Monte Carlo technique is applied to simulate the adsorption of long chain alkanes in zeolites. This simulation technique is several orders of magnitude more efficient than conventional methods that can be used to simulate the adsorption of long chain alkanes. The calculated heats of adsorption are found to be in excellent agreement with experimental data. The results show a surprising chain length dependence of the heats of adsorption. This dependence has a simple molecular explanation in terms of preferential siting of the long chain alkanes.