Numerical solution of NLO Q^2 evolution equations for spin-dependent   structure functions
  Numerical solution of DGLAP $Q^2$ evolution equations is studied for polarized parton distributions by using a ``brute-force" method. NLO contributions to splitting functions are recently calculated,and they are included in our analysis. Numerical results in polarized parton distributions and in the structure function $g_1$ are shown. In particular, we discuss how numerical accuracy depends on number of steps in the variable $x$ and in $Q^2$. 
M. Hirai, S. Kumano, and M. Miyama