To simulate the characteristics of magnetic exchange coupling at the interface of MR/TbCo thin films, the directions of magnetizations were calculated by minimizing energy in the films. Newton method and Gauss-Seidel method were used. The width of M-H curve increased with TbCo anisotropy constant, and with the thickness of the transition region of TbCo layer. Hysteresis loop width became extremely narrow (less than 10 Oe of coercivity), when the TbCo transition region length was $400;AA$. Also the hysteresis loop of films with low interfacial exchange coupling constant was similiar to that of short transition region length. When interfacial exchange coupling constant was 1/100 of perfect coupling, hysteresis loop showed a coercivity of less than 10 Oe. Comparing the measured hysteresis loop of a fabricated sample with that of simulated one, exchange coupling con¬stant could be estimated.