# Tarea#11 
# Fisica computacional (2:00-3:00 pm)
# Programa elaborado por Nalleli Iridian Avila Garcia (1941603)

import numpy as np
from scipy.integrate import quad

print('Este programa calcula los coeficientes de un polinomio de segundo grado que se ajusta a los datos de capacidad calorífica molar en función de la temperatura. Luego utiliza este polinomio para calcular el calor total absorbido por 23 moles de nitrógeno en un proceso isobárico cuasiestático. Finalmente calcula una cota del error absoluto cometido en el inciso anterior.')

# Aproximación de la función Cp(T) mediante una expresión polinomial de 2do grado
T = np.array([300,500,2100])
Cp = np.array([6.81495,8.25375,8.60055])
coeff = np.polyfit(T, Cp, 2)

# Calor total absorbido por 23 moles de nitrógeno en un proceso isobárico cuasiestático
Ti = 350
Tf = 800
Q, _ = quad(np.poly1d(coeff), Ti, Tf)
Q *= 23

# Cota del error absoluto cometido en el inciso anterior
error = abs((Tf - Ti) * 10**-13)

print(f'a) Coeficientes del polinomio: {coeff}')
print(f'b) Calor total absorbido: {Q}')
print(f'c) Cota del error absoluto: {error}')