Weight the "greater than" constraints very heavily, so rather than PO4>0, use PO4*10^13>0. However, if you add a constraint HPO4*10^13>0 it can't find a solution, so I think you have an error somewhere (perhaps a sign error?)
Пришлите, плз, сам файл с расчётом. Я его хочу в Прайме с единицами измерения прикрутить. Там, кстати, в последнем уравнении, как мне кажется, с.0 нужно убрать.
Да, Валерий Федорович, наше активное сотрудничество стало возможным именно благодаря этим логарифмам!
Прикрепляю расчетные файлы. Я совсем забыл, что химическую символику нельзя использовать в решающих блоках Mathcad 14, 15. Файл работать будет, но при попытке его сохранения появится ошибка:
Если проигнорировать это сообщение, то в дальнейшем такой файл уже открыть не удастся. Поэтому я отдельно прикрепляю файл без химической нотации.
По поводу C.0 в последнем уравнении. У меня ошибки нет. Я специально сверился с подобной задачей из Дж. Батлер. Ионные равновесия. - Л.:Химия, 1973.