The MATLAB files in this repository contain everything necessary to reproduce the results presented in
"A deterministic two-phase model for an active suspension with non-spherical active particles using Eulerian spatial averaging theory"
by B. Deußen, Y. Wang and M. Oberlack.
All files were created using MATLAB R2017a. Entry point is the file "StartHere.m", which can be run without any changes necessary. 
If you want to change parameter of the simulation (e.g. active forces) you can do this within the first lines of "AVShearPDE.m".
Inital conditions and boundary conditions can be changed in "InitalCond.m" respective "BoundaryCond.m".