Сравнить значения реальной численности населения в определенные годы со значениями, вычисленными по формуле в те же годы

import math
def arcctg(x):
    return math.pi / 2 - math.atan(x)
def compute_lambda(t):
    return 172 / 45 * arcctg((2000 - t) / 45)
 
sears = [1000, 1750, 1800, 1850, 1900, 1950, 1955, 
         1960, 1965, 1970, 1975, 1980, 1985, 1990, 1995, 
         2000, 2005, 2010, 2011, 2012, 2013, 2014, 2015,
         2016, 2017, 2018, 2019]
 
opulation =[0.400, 0.791, 1.000, 1.262, 1.650, 2.519,
             2.756, 3.021, 3.335, 3.692, 4.068, 4.435, 4.831,
             5.264, 5.674, 6.071, 6.344, 6.933, 7.015, 7.100,
             7.162, 7.271, 7.358, 7.444, 7.530, 7.669, 7.763]
 
per, vtor = int(input()), int(input())
t_list  = sears[per:vtor]
lambda_exp_list = opulation[per:vtor]
 
lambda_list =[compute_lambda(t) for t in t_list]
error_list = [abs((lambda_exp_list[i] - lambda_list[i]) / lambda_exp_list[i] )
              for i in range(len(t_list))]
 
max_error = max(error_list)
index_max_error = error_list.index(max_error)
min_error = min(error_list)
index_min_error = error_list.index(min_error)
avg_error = (sum(error_list) / len(t_list)) *100
print("Погрешность - минимальная, год: %4d, максимальная, год: %4d, средняя, процент: %6.3f" 
      % (t_list[index_min_error], t_list[index_max_error], avg_error))

t_list  = sears[per:vtor]
lambda_exp_list = opulation[per:vtor]

t_list  = sears[per:vtor+1]
lambda_exp_list = opulation[per:vtor+1]

import math
 
 
def arcctg(x):
    return math.pi / 2 - math.atan(x)
 
def calculate(t):
    return 172 / 45 * arcctg((2000 - t) / 45)
 
 
years = [1000, 1750, 1800, 1850, 1900, 1950, 1955,
         1960, 1965, 1970, 1975, 1980, 1985, 1990, 1995,
         2000, 2005, 2010, 2011, 2012, 2013, 2014, 2015,
         2016, 2017, 2018, 2019]
 
population = [0.400, 0.791, 1.000, 1.262, 1.650, 2.519,
              2.756, 3.021, 3.335, 3.692, 4.068, 4.435, 4.831,
              5.264, 5.674, 6.071, 6.344, 6.933, 7.015, 7.100,
              7.162, 7.271, 7.358, 7.444, 7.530, 7.669, 7.763]
 
start, end = int(input()), int(input()) + 1
 
error = [abs((p - calculate(y)) / p) for y, p in zip(years[start: end],
                                                     population[start: end])]
 
minimum = years[start + error.index(min(error))]
maximum = years[start + error.index(max(error))]
avg = sum(error) / len(error)
print("Погрешность - "
      f"минимальная, год: {minimum:4d}, "
      f"максимальная, год: {maximum:4d}, "
      f"средняя, процент: {avg:6.3f}")

import math
def c_l(t):
    return 172 / 45 * (math.pi / 2 - math.atan((2000 - t) / 45))
 
years, population  = [1000, 1750, 1800, 1850, 1900, 1950, 1955, 
         1960, 1965, 1970, 1975, 1980, 1985, 1990, 1995, 
         2000, 2005, 2010, 2011, 2012, 2013, 2014, 2015,
         2016, 2017, 2018, 2019], [0.400, 0.791, 1.000, 1.262, 1.650, 2.519,
             2.756, 3.021, 3.335, 3.692, 4.068, 4.435, 4.831,
             5.264, 5.674, 6.071, 6.344, 6.933, 7.015, 7.100,
             7.162, 7.271, 7.358, 7.444, 7.530, 7.669, 7.763]
 
per, vtor = int(input()), int(input()) +1
lil =[c_l(t) for t in years[per:vtor]]
err_list = [abs((population[per:vtor][i] - lil[i]) / population[per:vtor][i]) for i in range(len(years[per:vtor]))]
ime, imie, errs = err_list.index(max(err_list)), err_list.index(min(err_list)), (sum(err_list) / len(years[per:vtor])) *100
 
print("Погрешность - минимальная, год: %4d, максимальная, год: %4d, средняя, процент: %6.3f" 
      % (years[per:vtor][imie], years[per:vtor][ime], errs))