Abstract: In this paper we describe new algorithm of calculating total density of electromagnetic states in one-dimensional photonic dielectric crystals. Our algoritm do not use unitary transformations and have low computational complexity. We are implemented it on Fortran in sequential and parallel variants.