function psat=psat(t,const)
  for i=1:rows(const)
    psat(i)=exp(const(i,1)-const(i,2)/(t+const(i,3)));
  end
endfunction