Смешано V1 литров воды температуры t1 с V2 литрами воды температуры t2. Найти объём и температуру образовавшейся смеси

v1, t1, v2, t2 = map(int, input().split())
v = v1 + v2
print(f"{v}\n{(v1 * t1 + v2 * t2) / v}")

import numpy as np
 
spec_volume_poly=np.poly1d([-1.39698E-11,  5.98196E-09, 
                            -2.29081E-08, 9.99970E-04] )
 
spec_Cp_poly=np.poly1d([-3.525641E-11, 1.199009E-08, -1.554341E-06, 
                        1.034630E-04, -3.277896E-03, 4.216955E+00])
V1=1E-3*float(input('first volume, L: '))
V2=1E-3*float(input('second volume, L: '))
t1=float(input('first temperature, °C: '))
t2=float(input('second temperature, °C: '))
mass1=V1/spec_volume_poly(t1)
mass2=V2/spec_volume_poly(t2)
mass=mass1+mass2 # summary mass water
Q1=mass1*(np.polyval(np.polyint(spec_Cp_poly), t1)) #Energy water
Q2=mass2*(np.polyval(np.polyint(spec_Cp_poly), t2))
Q=Q1+Q2 
poly=np.polyint(spec_Cp_poly)
poly.c[-1]=-Q/mass
t=poly.r #roots
t=t[(t<max(t1,t2)) & (min(t1,t2)<t)].real[0] # delete complex roots
print('Temp=', t)

import numpy as np
 
spec_volume_poly=np.poly1d([-1.39698E-11,  5.98196E-09, 
                            -2.29081E-08, 9.99970E-04] )
 
spec_Cp_poly=np.poly1d([-3.525641E-11, 1.199009E-08, -1.554341E-06, 
                        1.034630E-04, -3.277896E-03, 4.216955E+00])
V1=1E-3*float(input('first volume, L: '))
V2=1E-3*float(input('second volume, L: '))
t1=float(input('first temperature, °C: '))
t2=float(input('second temperature, °C: '))
mass1=V1/spec_volume_poly(t1)
mass2=V2/spec_volume_poly(t2)
mass=mass1+mass2 # summary mass water
Q1=mass1*(np.polyval(np.polyint(spec_Cp_poly), t1)) #Energy water
Q2=mass2*(np.polyval(np.polyint(spec_Cp_poly), t2))
Q=Q1+Q2 
poly=np.polyint(spec_Cp_poly)
poly.c[-1]=-Q/mass
t=poly.r #roots
t=t[(t<max(t1,t2)) & (min(t1,t2)<t)].real[0] # delete complex roots
print(f'Temp={t:.4f}°C')
V=1000*spec_volume_poly(t)*mass
print(f'Volume={V:.4f} liters')