A general method using two times Green’s functions together with Linear Response Theory is proposed to calculate the relaxation times in nuclear spin systems. The method is applied to a nuclear spin coupled to a magnetic impurity, which is described by Bloch type equations of motion. the results obtained are in agreement with earlier theories and experimental results.