#!/bin/octave

# Druck über Wasser
function p=pressureoverwater(T)
p = exp(54.842763-6763.22/T-4.210*log(T)+0.000367*T+tanh(0.0415*(T-218.8))*(53.878-1331.22/T-9.44523*log(T)+0.014025*T));
return;
endfunction

# Druck über Eis (Achtung: T muss grösser 110K sein, wird nicht überprüft)
function p=pressureoverice(T)
p = exp(9.550426-5723.265/T+3.53068*log(T)-0.00728332*T);
return;
endfunction

pressureoverwater(273)

pressureoverwater(260)
pressureoverice(260)
