def forward_euler(initial_x, initial_y, num_steps, h):
    x = np.zeros([num_steps + 1, 2])
y = np.zeros([num_steps + 1, 2])
a = 1.0
x[0, 0] = initial_x[0]
y[0, 0] = initial_y[0]
x[0, 1] = initial_x[1]
y[0, 1] = initial_y[1]
for step in range(int(num_steps)):
    x[step + 1] = x[step] + h * y[step]
y[step + 1] = y[step] + h * (-x[step] - a * y[step] * (x[step] ** 2 + y[
    step] ** 2 - 1))
return x, y
Векторизовать прямой метод Эйлера для системы дифференциальных уравнений