---

# ============================================================
#  Сравнение Системы 1 и Системы 2 при h = 0.025
# ============================================================
print("=" * 72)
print("Сравнение Системы 1 и Системы 2 при h = 0.025")
print("=" * 72)
print(f"{'x':>8} {'Система 1':>16} {'Система 2':>16} {'|разность|':>14}")
print("-" * 72)

# Сетка у обеих систем одинаковая (x1f == x2f), берём x1f
for xi, ya, yb in zip(x1f, y1f, y2f):
    print(f"{xi:>8.4f} {ya:>16.10f} {yb:>16.10f} {abs(ya - yb):>14.3e}")
print()

# Максимальная разница между системами при h = 0.025
diff_coarse = max(abs(ya - yb) for ya, yb in zip(y1c, y2c))
diff_fine   = max(abs(ya - yb) for ya, yb in zip(y1f, y2f))

print("=" * 72)
print("Как меняется разница между Системой 1 и Системой 2")
print("при уменьшении шага")
print("=" * 72)
print(f"max|С1 - С2| при h = 0.05 : {diff_coarse:.6e}")
print(f"max|С1 - С2| при h = 0.025: {diff_fine:.6e}")
print(f"Отношение (должно быть ≈ 4 для O(h²)): {diff_coarse / diff_fine:.2f}")
print()