编写2D相图

Question

我正在尝试编写带有签名f(pressure, temperature)的函数，该函数返回一个相（即，固体，液体或蒸气）。

问题在于此二维函数非常复杂。例如，铁具有以下相图：

表面上我可以将其写为许多if语句，但是有更好的方法吗？这是使用偏微分方程的合适应用程序吗？如果可以，我将如何用（任何一种）编程语言来表示它们？