Question: How to compute flow rate, impedance, and shear stress?

I am solving a hybrid nanofluid flow problem in a bifurcated artery using Maple. The governing equations for velocity and temperature are solved using dsolve(..., numeric, method=bvp[midrich]).

My Maple code successfully produces velocity profiles for different Hartmann numbers Ha for both the artery  parentartery_and_daughter_artery_error.mw.

The velocity profiles are obtained correctly using odeplot.

However, I want to compute additional physical quantities and generate plots similar to the velocity profiles.

Specifically I want to plot:

1. Flow rate Q versus axial distance z

2. Impedance (flow resistance) λ versus z

3. Wall shear stress τ versus z

for different values of Hartmann number Ha.

The formulas I am using are

Flow rate:

Q=2π(R∫01ηw(η) dη+R2∫01w(η) dη)Q = 2\pi \left( R \int_0^1 \eta w(\eta)\,d\eta + R_2 \int_0^1 w(\eta)\,d\eta \right)Q=2π(R∫01​ηw(η)dη+R2​∫01​w(η)dη)

Wall shear stress:

τ=μ∣dwdr∣\tau = \mu \left|\frac{dw}{dr}\right|τ=μ​drdw​​

Impedance:

λ=∣dp/dz∣Q\lambda = \frac{|dp/dz|}{Q}λ=Q∣dp/dz∣​
Please help me to solve this question.