An approach is presented to simulate energy filtered inelastic images of crystalline materials with atomic resolution. To this purpose the multislice method (used to calculate multiple elastic scattering events) was adopted in such a way that it includes inelastic scattering events. The approach is described with a density matrix formalism. The implementation into the software program is explained and simulation results for SrTiO3 and Si are presented. The program shows a simple way to parallelize the calculations and proves that it is possible to run simulations on home computers in affordable time.