1.1

import matplotlib.pyplot as plt
from math import log, exp
import pandas as pd
import numpy as np

T_s = 22.0 #температура окружающей среды
# экспериментальные данные
T = np.array([83, 77.7,75.1,73,71.1,69.4,67.8,66.4,64.7,63.4,62.1,61,59.9,58.7,57.8,56.6])
t = np.arange(0,16,1)

def NewtonHeat(t0,T0,Ts,r,t,M):
    t = np.linspace(t0,t0+t,M)
    return t, Ts - (Ts-T0)*np.exp(-r*t)

# Вычислим коэффициенты r по формуле
r = np.divide(np.log((T[0]-T_s)/(T[1:]-T_s)),t[1:])
# Проверим
T_r = T_s + (T[0] - T_s)*np.exp(-r*t[1:])

plt.plot(t[1:], r)
plt.ylabel('коэффициент r')
plt.xlabel('время t')
plt.show()