Вариант 4 / 5 / lab5chm

import math

def fy(x, y):
return (x)/(2*y)

def r(x, y, h, m):
for j in range(m):
k1 = fy(x, y)
k2 = fy(x + h / 2, y + h * k1 / 2)
k3 = fy(x + h / 2, y + h * k2 / 2)
k4 = fy(x + h, y + h * k3)

y = y + (h / 6) * (k1 + 2 * k2 + 2 * k3 + k4)
x = x + h

return y

def diff(x0, y0, b, h0, eps):
n = int((b - x0) / h0)

for i in range(1, n+2):
h = h0
m = 1
y = r(x0, y0, h, m)

while True:
y1 = y
h = h / 2
x = x0
y = y0
m *= 2
y = r(x, y, h, m)

if abs(y - y1) / 15 < eps:
break

x0 += h0
y0 = y

print(f"x0: {round(x0, 2)}")
print(f"y0: {y0}")
print()

x0 = 0
y0 = 1
b = 1.2
h0 = 0.4
eps = 0.0001

diff(x0, y0, b, h0, eps)