Can you write a function to find Q1, Q2, and Pa in the question 1-5 using matlab preferably or C++, Java, etc

EXAMPLE1S Cossider the parallel fow network of Fig 1-20 with the following specifications | = 0.001 ft Z.-80ft e2 _ 0.0001 ft 0.001 0,00015 Di Find Qi.Q.and P A. -000003 f/s (2 FIGURE 1-20 Parallel flow network for Example 1-5 Desktop