In the case of the structure function \(F_2\)  the non–singlet anomalous dimensions are \(P^{{\rm NS},+}\)  and the
expansion coefficients of the Wilson coefficient up to \(c_3\)  were calculated in [1]}.
In the case of the structure function \(g_1\)  the anomalous dimensions are \(P^{{\rm NS},-}\) , cf. [2]},
and the Wilson coefficient has been given in [3]}.Lower orders of the non–singlet anomalous
dimensions and
the massless non–singlet Wilson coefficients have been calculated in
[4]}, [5]}, [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, [12]}, [13]}, [14]}
and
[15]}, [16]}, [17]}, [18]}, [12]} and in the polarized
case in [20]}, [21]}, [22]}, [23]}, [24]}, [25]} and
[26]}, [27]}, [28]}, [29]}.
Below we will also use the combination
\(F_2^{\rm h}(N,Q^2) = \left[E_{\rm NS} - \left. E_{\rm NS}\right|_{\rm h=0}\right] F_2(N,Q_0^2).\) 
