function vp = sat_pr_antoine(A, B, C, T) % Vap_Pressure = sat_pr_antoine(A, B, C, T) % Vapor Pressure for a given Temperature % Input: Antoine Constants (A,B,C), T % Output: Vapor Pressure % ln P = A - B/(T+C) % A, B, C constants of Antoine equation. % T temperature % Vapor pressure from Antoine Equation % Done on 09-Sep-2001; M.Subramanian % vp = Vapor pressure % Vectorized (the dot next to B) vp = exp(A - B./(T+C)); % end