Code using Python

HW6_3 For a system of nine interconnected tanks (shown at right), the massbalance equations are: r01c01+r41c4r02c02+r12c1+r52c5r03c03+r23c2r04c04+r74c7r15c1+r45c4r06c06+r26c2+r36c3r07c07+r87c8r08c08+r58c5+r48c4r59c5+r69c6+r89c8=(r12+r15)c1=(r23+r26)c2=r36c3=(r41+r48+r45)c4=(r52+r58+r59)c5=r69c6=r74c7=(r87+r89)c8=r9 The concentration of chemical inflow into tank Ti is designated ci and the initial concentration in tank Ti is designated c0i. The flow rate from tank Ti to Tj is designated rij. Furthermore, the total flow into each tank must equal the flow out, r (from tank 9 ). Given initial concentrations: c01=1,c02=2,c03=1,c04=2,c06=1,c07=1 and c08=1, in g/L given flow rate between tanks: r12=6,r15=2,r23=5,r26=5,r36=8,r41=5,r45=1,r48=1,r52=1,r58=2,r59=1,r69=15,r74=4,r87=2,r89=3, in L/min and given external flow rates (i.e., inflow into the tanks from outside): r01=3,r02=3,r03=3,r04=3,r06=2,r07=2,r08=3, in L/min Note: r=r01+r02+r03+r04+r06+r07+r08 With these given parameters we can solve for the concentrations in the tanks. Write a PYTHON code to solve for these concentrations using Gauss-Seidel, to an accuracy of 5 sig figs. Use print statement to print the final answer for the concentrations