1.5

def find_eps(t):
    M = 100
    delta_M = 1
    r_f = 0.09
    while (True):
        if (np.abs(NewtonHeat(0,85,22,r_f, t, M)[1][-1] - HeatEq(0,85,22,r_f, t, M)[1][-1]) <= 0.01):
            return M,t/M
        M += delta_M

    return False
res1 = find_eps(1)
res2 = find_eps(5)
print("Ответ")
print(f"при t = 1 M = {res1[0]}, Δt = {res1[1]}")
print(f"при t = 5 M = {res2[0]}, Δt = {res2[1]}")