Hölder continuity of weak solutions is studied for a nondiagonal parabolic system of singular quasilinear differential equations with matrix of coefficients satisfying special structure conditions. A technique based on estimating linear combinations of the unknowns is employed.