def polynomial\_interpolation(x\_data, y\_data, x\_new):
pol = np.poly1d(np.polyfit(x\_data, y\_data, 3))
return pol(x\_new)
x\_data = np.array([date2num(d) for d in ["01/01/2022", "02/01/2022", "03/01/2022",
"04/01/2022", "05/01/2022"]])
y\_data = np.array([10000, 12000, 14000, 16000, 18000])
x\_new = date2num("06/01/2022")
pol = np.polyfit(x\_data, y\_data, 3)
new\_value = polynomial\_interpolation(x\_data, y\_data, x\_new)
new\_date = num2date(x\_new)