r <- 19.075 #Radio del objeto medido con pie de Rey en cm
z <- 61.1 #Espesor del objeto medido con pie de Rey en cm
m <- 188 #masa del objeto medido en gramos
Deltar <- 0.05 #Resolución del pie de rey en mm

Deltam <- 1 #Resolución de la bascula en gramos

fr <- 2*pi*r*z #Derivada de la función volumen con respecto a r
fz <- pi*r^2 #Derivada de la función volumen con respecto a z

Vol=pi*r^2*z/1000 #Volumen del cilindro en cm³
Error=sqrt((fr*Deltar)^2+(fz*Deltar)^2)/1000 #Error del volumen 
 
fm <- 1/Vol
fv <- -m/Vol^2
Deltav=Error

Den=m/Vol #Densidad del cilindro en g/cm³
Error2=sqrt((fm*Deltam)^2+(fv*Deltav)^2) #Error del volumen 

cat("Volumen:", Vol, "\n")
cat("Error del volumen:", Error, "\n")
cat("Volumen por exceso:", Vol+Error, "\n")
cat("Volumen por defecto:", Vol-Error, "\n")
cat("Error relativo del volumen:", 100*Error/Vol, "\n" )
cat("\n")
cat("Densidad:", Den, "\n")
cat("Error de la densidad:", Error2, "\n")
cat("Densidad por exceso:", Den+Error2, "\n")
cat("Densidad por defecto:", Den-Error2, "\n")
cat("Error relativo de la Densidad:", 100*Error2/Den, "\n" )