  We develop a systematic expansion method of physical quantities for the SK
model and the finite-dimensional $\pm J$ model of spin glasses in
non-equilibrium states. The dynamical probability distribution function is
derived from the master equation using a high temperature expansion. We
calculate the expectation values of physical quantities from the dynamical
probability distribution function. The theoretical curves show satisfactory
agreement with Monte Carlo simulation results in the appropriate temperature
and time regions. A comparison is made with the results of a dynamics theory by
Coolen, Laughton and Sherrington.
