import numpy as np
def find_density(T, P):
coeficientes = np.array([
569.502E-6 * 2.1224E-6 * T - 513.2105E-9, # r^3
738.22E-6 - 569.502E-6 * 2.507E-3 * T, # r^2
-569.502E-6 * T, # r^1
P # r^0
])
soluciones = np.roots(coeficientes)
for s in soluciones:
if s.imag == 0.0:
return s.real
return np.nan()