• Author / Uploaded
• Laura Viviana Martinez

Citation preview

Volumen function y = volumen(longitud,radio,altura) a=longitud*((acos((radio-altura)/radio)*radio^2)); b=longitud*(sqrt((radio^2)-(radio-altura)^2)*(radio-altura)); y=a-b; end

Arrhenius function y= arrhenius(temperatura) y=(2.46E10)*exp((-22000/(1.987*temperatura))); end

ArrheniusLin function y = arrheniusLin(temperatura) dx=0.000001; y1=arrhenius(temperatura+dx); y2=arrhenius(temperatura-dx); dydx=(y1-y2)/(2*dx); y=100+(dydx*(temperatura-573)); end

DensidadLin function y = densidadLin(temperatura,presion,pesomol,CTE) dP=0.000001; dT=0.000001; den1P= densidad(presion+dP,temperatura,pesomol,CTE); den2P= densidad(presion-dP,temperatura,pesomol,CTE); ddenP=((den1P*den2P)/(2*dP)); den1T=densidad(presion,temperatura+dT,pesomol,CTE); den2T=densidad(presion,temperatura-dT,pesomol,CTE); ddenT=((den1T*den2T)/2*dT); y= densidad(101.9,300,pesomol,CTE)+ddenP*(presion-101.9)+ddenT*(temperatura-300); end

densidad function y = densidad(presion,temperatura,MW,R) y=(MW*presion)/(R*temperatura); end