It's one thing to estimate the amount of phosphate versus hydrogen phosphate from the Henderson-Hasselbalch equation, but how do I set up an expression derived from multiple Ksps, if I want to have an expression for the maximum pH-dependent concentration of any phosphate species ( tribasic phosphate and dibasic phosphate?