invalid syntax

import matplotlib.pyplot as plt
import numpy as np
import sys
 
x0=float(input("Введите начальное значение x "))
y0=float(input("Введите начальное значение y "))
a=int(input("Введите начало отрезка "))
b=int(input("Введите конец отрезка "))
h=float(input("Введите шаг "))
n = round((b - a) / h) + 1
x=[]
y=[]
y.insert(0,y0)
x.insert(0,x0)
 
e=sys.float_info.epsilon
 
 
def f (x0,y0):
    return pow(e, x0)*np.cos(2*pow(y0,2))
 
def Euler(x, y, h, n):
    for i in range (0,n-1):
        y.insert(i + 1, y[i] + h * f(x[i], y[i])
        x.insert(i + 1, x[i] + h)
    return(y)
 
def ModifiedEulerMethod(x, y, h, n):
    for i in range (0,n-1):
        y1 = y[i] + h / 2 * f(x[i], y[i])
        y.insert(i + 1, y[i] + h * f(x[i] + h / 2, y1))
        x.insert(i + 1, x[i] + h)
    return(y)
 
def EulerMethodWithRecalculation(x, y, h, n):
    for i in range (0,n-1):
        y1 = y[i] + h * f(x[i], y[i])
        y.insert(i + 1, y[i] + h / 2 * (f(x[i], y[i]) + f(x[i] + h, y1)))
        x.insert(i + 1, x[i] + h)
     return(y)
 
x = np.linspace(-4, 4, 81)
fig, (ax1, ax2, ax3) = plt.subplots(1, 3, figsize=(20, 5))
ax1.set_xlim([-1, 1])
ax1.set_ylim([-1, 1])
ax2.set_xlim([-1, 1])
ax2.set_ylim([-1, 1])
ax3.set_xlim([-1, 1])
ax3.set_ylim([-1, 1])
 
ax1.set_xlabel('Euler')
ax1.plot(x, f(x0,y0), color='red', linewidth=0.5)
ax1.plot(x, Euler(x0, y0, a, b, h), color='black', linewidth=0.5)
                    
ax2.set_xlabel('ModifiedEulerMethod')
ax2.plot(x, f(x0,y0), color='red', linewidth=0.5)
ax2.plot(x, ModifiedEulerMethod(x0, y0, a, b, h), color='black', linewidth=0.5)
 
ax3.set_xlabel('EulerMethodWithRecalculation')
ax3.plot(x, f(x0,y0), color='red', linewidth=0.5)
ax3.plot(x, EulerMethodWithRecalculation(x0, y0, a, b, h), color='black', linewidth=0.5)
plt.show()

x0=float(input("Введите начальное значение x "))
y0=float(input("Введите начальное значение y "))
x=[]
y=[]
y.insert(0,x0)
x.insert(0,y0)
print(x)
print(y)

y.insert(i + 1, y[i] + h * f(x[i], y[i]))