# Cálculo de la entalpía de vaporización a partir de datos de Temperatura y presión.
P=c(17.54,31.82,55.32,92.51,149.38,233.7) # Introduce valores de Presión (Hgmm)
T= c(20,30,40,50,60,70) # Introduce valores de Temperatura (°C)
Tk=T+273.15 # Transforma grados Celcius a Kelvin
lnP=log(P) # Calcula el logaritmo natural de la Presión
iT=1/Tk # Calcula el inverso de temperatura (K)
reg=lm(lnP~iT) # Calcula la regresión lineal
plot(iT,lnP, abline(reg, col = "green"), col="red", xlab = "1/T (°C)",
ylab="Ln(P-mmHg)") # Grafica los valores transformados 1/T vs lnP y la pendiente
summary(reg) # Visualiza los valores del modelo linear
c=coef(reg) # Extrae los coeficientes de la regresión linear
b=signif(c[1]) # Asigna valor a "b"(Intercept) a partir de "c"
m=signif(c[2]) # Asigna valor a "m"(pendiente)
R= 8.314 # J*K^-1*mol^-1 , introduce el valor de R
dH= m*R*-1 # calcula la entalpía de vaporización
dH # visualiza la entalpía de vaporizacion en Julios
dHkJ=dH/1000 # convierte los Julios en kiloJulios
dHkJ # visualiza el valor entalpía de vap en kilojulios
No hay comentarios.:
Publicar un comentario