Análisis experimental del principio de Arquímedes

import numpy as np

# Datos
ms = 92.7  # masa de los objeto: plomo 227.3, hierro 433.9, aluminio 149.3 y laton 92.7
me = 10.9  # masa de empuje: plomo 20.1, hierro 55.7, aluminio 55 y laton 10.9
pa = 0.997  # densidad del agua en g/cm³

# Sensibilidades de la balanza
Deltams = 0.1
deltame = 0.1

# Cálculo de la densidad del objeto problema
ps = (ms / me) * pa

# Derivadas parciales para el cálculo del error
fms = pa / me
fme = -(ms / me**2) * pa

# Error de la densidad
deltaps = np.sqrt((fms * Deltams)**2 + (fme * deltame)**2)

# Impresión de resultados
print(f"Densidad de la sustancia en (g/cm³): {ps:.4f}")
print(f"Error de la densidad: {deltaps:.4f}")
print(f"Error relativo: {100 * deltaps / ps:.2f} %")
print(f"Densidad por exceso: {ps + deltaps:.4f}")
print(f"Densidad por defecto: {ps - deltaps:.4f}")