いつになったらみんなが電気自動車を運転するようになるのでしょうか？

import matplotlib.pyplot as plt
import numpy as np
import math

x = np.linspace(2012, 2019, 8)
y1 = np.array([60936407, 63429200, 65708230, 66314155, 69464432, 70694834, 68690468,  64341693]) # кол-во произведенных легковых машин
y2 = np.array([52605, 97507, 320713, 550297, 777495, 1227117, 2018247,  1940147]) # кол-во произведенных легковых электромобилей
y = y2/y1 #доля электромобилей в общем производстве автомобилей

ymax=1 #первоначальное максимальное отклонение статистических данных от значений функции
Gmax=2025 #год для начало поиска "нулевого года"
rmax=0.35 #начальный коэффициент
k=1 #принят "1" из предпосылки, что электромобили полностью заменят легковые автомобили с ДВС
p0=0.5 # процент рынка в "нулевой год"
for j in range(10): # цикл перебора "нулевых годов"
    x0=2025+j
    r=0.35
    
    for i in range(10): # цикл перебора коэффициента в каждом "нулевом году"
            r=0.25+0.02*i
            y4=k*p0*math.e**(r*(x-x0))/(k+p0*(math.e**(r*(x-x0))-1))-y 
           # print(str(x0).ljust(20), str(r).ljust(20), max(abs(y4))) 
            if max(abs(y4))<=ymax: # поиск минимального из максимальных отклонений внутри каждого года при каждом коэффициенте r
                ymax=max(abs(y4))
                Gmax=x0
                rmax=r
print(str(Gmax).ljust(20), str(rmax).ljust(20), ymax) # вывод "нулевого года", коэффициента r и максимального из отклонений от функции

import matplotlib.pyplot as plt
import numpy as np
import math

x = np.linspace(2012, 2019, 8)
y1 = np.array([60936407, 63429200, 65708230, 66314155, 69464432, 70694834, 68690468,  64341693])
y2 = np.array([52605, 97507, 320713, 550297, 777495, 1227117, 2018247,  1940147])
y = y2/y1

k=1
p0=0.5

x0=2028   
r=0.37 
y1=k*p0*math.e**(r*(x-x0))/(k+p0*(math.e**(r*(x-x0))-1))
#Строим график функции на отрезке между 2012 и 2019 годами
fig, ax = plt.subplots(figsize=(30, 20), facecolor="#f5f5f5")
plt.grid()
ax.plot(x, y, 'o', color='tab:brown') 
ax.plot(x, y1)
#Строим график функции на отрезке между 2010 и 2050 годами
x = np.linspace(2010, 2050)
y2 = [k*p0*math.e**(r*(i-x0))/(k+p0*(math.e**(r*(i-x0))-1)) for i in x]
y3 = 0.99+0*x
fig, ax = plt.subplots(figsize=(30, 20), facecolor="#f5f5f5") 
ax.set_xlim([2010, 2050])
ax.set_ylim([0, 1])
plt.grid()             
plt.plot(x, y2, x, y3)