5. (20%) Solve coupled differential equations, subjected to the constraints represent positive flow quantity in the domain is a positive constant. The following boundary conditions are also enforced,