Consider the two chemical reactors in series as shown in the picture above ( Q is the flow rate), which are operated at different temperatures to ensure the conversion for the irreversible first order chemical reactions as below (where all species have the same density): ABC There is a recycle stream Q21 from reactor 2 back into the first reactor in order to increase the overall yield of the process. The consumption of a species within the respective reactors depends on the reaction rate constant, which depends on temperature and thus is different in the respective reactors ( kB,1 for the production of B in reactor 1,kB,2 for the production of B in reactor 2,kC,1 for the production of C in reactor 1 and kC,2 for the production of C in reactor 2 ), the respective reactor volumes (V1 and V2 for reactor 1 and 2) and the concentration of the reacting species in the reactor through the general equation: R=kVc For example, for the rate at which A is transformed in reactor 1 , we can write: RAB1=kB,1V1cA,1 During all of our calculations we will assume that the system is at steady state. - Write a Python code for the Gauss-Seidel method to iteratively solve the system of equations to find the concentrations of species in the case of Qin=20L/hr, with an s=109 and initial guesses for all concentrations of 0 moll . Check your solution by comparing it to the solution found under Part 2. (15 points) - How many iterations are required to find the solution using this method? (5 points) - Compare the current method with your previously used method and comment on the differences in both methods. (15 points)