oleksandr-danylchenko

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint

# Визначаємо функцію рівняння із умови
def func(y, x):
    return (np.sin(2 * x) / (2 - np.sin(x) ** 2) - 2 * y) / 2

# Реалізація методу Ейлера-Коші де крок ітерації визначено як:
# x[i+1] = x[i] + h