An integrated numerical system is established to model time-dependent behavior of the neoclassical tearing mode (NTM) in a tokamak which solves the modified Rutherford equation (MRE) by coupling with plasma transport, equilibrium, heating and current drive self-consistently. The MRE is formulated by considering the electron cyclotron heating (ECH) effect as well as the electron cyclotron current drive (ECCD) effect to stabilize the NTM which is well-suited for time-dependent simulations including a predictive purpose. The integrated numerical system is applied to experiments in which NTMs are stabilized by ECCD and by ECH at ASDEX Upgrade and TCV, respectively. It is found that the simulation results catch the main experimental trend.