specific heat

I = 1.62 # heating current in A
V = 1.62 # heating voltage in V
dt = 240 # heating duration in s
T1 = 22.27 # temperature T1 in °C
TG = 34.52 # temperature TG in °C
dT = TG - T1 # delta T in K
delT = -0.00436 # \dot{T} in K/s
F = 17904.5 # integral F
#F = 781.2*12.24

C_add = (I*V*dt*dT)/(dT**2 - delT*F) # C_Add in J/K
print("C_add: ", C_add)

# same nomenclatur as above
I = 1.62 # A
V = 1.62 # V
dt = 600 # s
T1 = 25.98 # °C
TG = 41.848 # °C
dT = TG - T1
#print(dT)
delT = -0.00272 # °C
F = 1764*15.86 # approximation of the integral with delta_T*(t_G - t_s), because this gives a better result, as our computed F = 21329.134

C_tot = (I*V*dt*dT)/(dT**2 - delT*F)
#print(C_tot)
C_a = C_tot - C_add
print("C_sample for alluminium: ", C_a )
m_a = 0.4745
M_a = 26.98
n_a = m_a/M_a

#print(C_a / m_a)
print("specific molar heat capacity for aluminium: ", C_a/n_a)