Nonlinear spinor field in Bianchi type-I Universe filled with viscous fluid: numerical solutions

We consider a system of nonlinear spinor and a Bianchi type I gravitational fields in presence of viscous fluid. The nonlinear term in the spinor field Lagrangian is chosen to be λ F, with λ being a self-coupling constant and F being a function of the invariants I an J constructed from bilinear spinor forms S and P. Self-consistent solutions to the spinor and BI gravitational field equations are obtained in terms of τ, where τ is the volume scale of BI universe. System of equations for τ and ε, where ε is the energy of the viscous fluid, is deduced. This system is solved numerically for some special cases.