Cooling Tower Problem Coding In Python Part 1

Solve the differential equations derived in B1 for a 20 m tall cooling tower where the air-water heat transfer coefficient is 10 W(m2 K), specific area of 5 (m2/m3). Hot water enters from the top at 70 celsius with a flow rate of 0.5 kg/(m2 s). Air enters from the bottom with a velocity of 5 m/s, at a temperature and relative humidity of 25 celsius and 50% respectively. Make plots of the water temperature, air temperature, and water partial pressure as a function of position in the cooling tower.