Using the recently proposed concept of B(PI), the individual photon's magnetostatic flux density operator, it is argued in semi-classical approximation that magnetization by a circularly polarised laser (the inverse Faraday effect) takes place through terms proportional to all positive integral powers of B(PI), the classical equivalent of the quantum field operator B(PI). The original theory of the IFE considers only the term in the square of B(PI).