SHG and MSHG of bcc Fe, hep Co and fee Ni magnetic films with different surface, direction of magnetization and input-output polarization combinations are derived. The image of nonlinear magneto-optical response of a modeling domain structure is obtained, which can explain some recent experimental facts.