Can you write a function to find Q1, Q2, and Pb in the question 1-5 using python using this solution as a roadmap? I can't figure it out.

Consider the parallel flow network of Fig. 1-20 with the following specifications: L1=3,000ftD1=1ft1=0.001ftD11=0.001PA=80psiaZA=100ftZB=80ftL2=3,000ftD2=8in.2=0.0001ftD22=0.00015 Find Q1,Q2, and PB. QA=5.3ft3/s ,UJRE 1-20 Parallel flow network for Example 1-5. Solution. This is a type 2 problem. Step 1: Assume Q1=3ft3/s. Apply the energy equation along line 1 from A to B to obtain Wsggc=zbza+28gD14Q12[f1D1L1]=zbza+hf1ggc assuming that VA=VB. Finding hf1 with Q1 specified is a category I problem, and we have V1ReD1=A1Q1=3.82ft/s=1V1D1=1.273105 So f1=0.022, and hf1=28gD14Q12[f1D1L1]=14.97lbmftlbf Step 2: The loss for pipe 2 must then become hf2=hf1=14.97lbmftlbf Finding Q2 for hf2=14.97ftlbf/lbm is a category II problem. We shall use V2 as the iteration variable. The results are given in Table 1-4. Step 3: From step 2, Q2=V2A2=1.141ft3/s and i=12Qi=1.141+3=4.141ft3/s. We find the corrected values by using Q1=4.14135.3=3.84ft3/sQ2=4.1411.1415.3=1.46ft3/s Step 4: Using Q1=3.84ft3/s and Q2=1.46ft3/s, we compute V1,V2,ReD1,ReD2,f1, and f2, to find hf1=23.85lbmftlbfhf2=23.98lbmftlbf The head losses in the pipes then agree to 0.54 percent, which is sufficient. With the flow rates in each line and the head loss for the parallel segments known, Ws is computed from the first equation in step 1. WsggcWs=zbza+hf1ggc=80ft100ft+23.91lbmftlbf32.174s2ft32.174lbfs2ftlbm=20ft+23.91ft=3.91ft=3.91lbmftlbf Thus, a pump with an increase in head of 3.91ftlbf/lbm would be required for the system to pass 5.3ft3/s while maintaining PA=PB