in Matlab and Mathcad for Steam properties such as dH/dP , dH/dT , dS/dP , dS/dT , dV/dP , dV/dT ....
I utilize IAPWS95 model and a process library (Prode Properties),
IAPWS95 is the most accurate but it is complex to code,
you may find more easy IAPWS97 or use an external code as I do.