The time-evolution of the damage in the Ising model with the heat bath and Glauber dynamics has been investigated with the damage spreading method. The damage density D( t) and the average distance r( t) of damaged sites from the origin have been measured numerically. From the simulational data, the scaling forms have been established as D(t) ∼ t x, r(t) ∼ t 1 z and D( t) ∼ r( t) d f and the corresponding exponents x, z and d f have been estimated.