This code implements the evolution equations derived in Refs. [1] & [2]. It computes the iterative time evolution for N real scalar fields with quartic self-interaction coupled to two flavours of fermions by a Yukawa coupling. [1] @article{Berges:2001fi, author = "Berges, Juergen", title = "{Controlled nonperturbative dynamics of quantum fields out-of-equilibrium}", eprint = "hep-ph/0105311", archivePrefix = "arXiv", reportNumber = "HD-THEP-01-26", doi = "10.1016/S0375-9474(01)01295-7", journal = "Nucl. Phys. A", volume = "699", pages = "847--886", year = "2002" } [2] @article{Aarts:2002dj, author = "Aarts, Gert and Ahrensmeier, Daria and Baier, Rudolf and Berges, Juergen and Serreau, Julien", title = "{Far from equilibrium dynamics with broken symmetries from the 2PI - 1/N expansion}", eprint = "hep-ph/0201308", archivePrefix = "arXiv", reportNumber = "HD-THEP-02-5", doi = "10.1103/PhysRevD.66.045008", journal = "Phys. Rev. D", volume = "66", pages = "045008", year = "2002" }