We present a numerical method to solve the Vlasov–Maxwell equations for spin-1/2 particles, in semiclassical approximation where orbital motion is treated classically while spin variable fully quantum. Unlike spinless case, phase-space distribution function $2\times 2$ matrix, which can also be represented, Pauli basis, as one scalar $f_0$ and three-component vector $\boldsymbol f$ . The relati...