Approximation of Hopkins equation in Hilbert space and its application in polarized illumination modelling