This Letter presents a new method for modeling of complex optical setups illuminated by quasi monochromatic spatially incoherent light. The algorithm provides better performance and quality than other modeling methods both for isoplanatic and nonisoplanatic systems. The algorithm maintains energy relations, image orientation, and magnification of the system. Computer modeling and experimental results are presented.