Solving coupled equations
 

Hi all,
I am developing a 2dv hydrodynamic solver and I have a problem on coupling EARSM and EAHFM in turbulence modelling. Can anyone help me on simultaneously solving of 5 coupled equations for 1-stress anisotropy tensor, 2-P/eps, 3-scaled heat flux, 4-Ptheta/epstheta and 5-G/eps in iterative manner?
thank you