metodos

# metodo de simpsom
olmer trujillo, jessica tatiana quintero, juan camilo arriola mahecha
# 1. definimos variables
p=c(0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90,95, 100)
z=c(0.9715, 0.9845, 0.9155, 0.8881, 0.8083, 0.7354, 0.5939, 0.4411, 0.2704, 0.3814)
c= coef(lm.fit(outer(p,0:4,'^'),z))
# 2. graficamos con el comando plot (p, z)
plot(p,z)
type="l""main= z vs p"
lines(p,z)
rug(p)
par(bg="white")
col="blue"f
bg=" yellow"
plot(p, z, xlab="Presión (bar)", ylab="Compresibilidad (z)", lines (p,z),pch=19, col="black",bg="blue", bty="l",tcl=0.4,main="Compresibilidad vs Presión", las=1, cex=1.5)
# 3. Resolución de la Integral por el método de Simpson (mayor aproximación al valor real)
Z=6E-11*p^4 -2E-8*p^3 +1E-5*p^2 -6.3E-3*p +1.0137
a=1b=90
h=(b-a)/2
fa=(6E-11*a^4 -2E-8*a^3 +1E-5*a^2 -6.3E-3*a +1.0137 -1)/a
fa #aqui es el resultado de lo anterior

fb=(6E-11*b^4 -2E-8*b^3 +1E-5*b^2 -6.3E-3*b +1.0137-1)/b
fb #aqui es el resultadode lo anterior

fh=(6E-11*h^4 -2E-8*h^3 +1E-5*h^2 -6.3E-3*h +1.0137 -1)/h
fh #aqui es el resultado de lo anterior

I=(b-a)*((fa +4*fh +fb)/6)
I #aqui es el resultado de lo anterior

# 4.Antilogaritmo para poder hallar el coeficiente de fugacidad.
#f=Coeficiente de fugacidad
f=exp(I)

# 5. realizacion de la graficas:
donde:
b=(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100)

a=1b=
h=(b-a)/2
fa=(6E-11*a^4 -2E-8*a^3 +1E-5*a^2 -6.3E-3*a +1.0137 -1)/a
fb=(6E-11*b^4 -2E-8*b^3 +1E-5*b^2 -6.3E-3*b +1.0137-1)/b
fh=(6E-11*h^4 -2E-8*h^3 +1E-5*h^2 -6.3E-3*h +1.0137 -1)/hI=(b-a)*((fa +4*fh +fb)/6)
f=exp(I)

los resultados obtenidos de la fugacidad al reemplazar (b) son:
y=1.004128, 0.9924423, 0.9657716, 0.9473766, 0.8954425, 0.848017, 0.7648644, 0.6740026, 0.5438123,0.453595

#grafica: coeficiente de fugacidad (x) en funcion de la presion (y)
#vectores
x=c(1.004128, 0.9924423, 0.9657716, 0.9473766, 0.8954425, 0.848017, 0.7648644, 0.6740026, 0.5438123,...