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$.