Solving a system of coupled PDEs in fenicsx

I would suggest you have a look at: Default absolute tolerance and relative tolerance - #4 by nate