A Monte Carlo simulation of a nematic Lebwohl-Lasher liquid crystal model with a certain percentage of randomly oriented frozen spins is presented. The site-site distance dependence of the nematic correlation functions obtained are generally well fitted by stretched exponentials, which evolve into simple exponential at low temperature. The dependence of the correlation length on the disorder density is presented and compared to corresponding experimental results.