Create an image visualising the problems resolved in the following code

import math

# данни
avans = 5 # мм за не пълен ход на амортисьора
L_closed = 180.0 + avans
L_open   = 260.0 - avans


def gas_strut_force(
    weight_kg = 8.0,        # тегло на вратата в кг
    length_mm = 500.0,      # дължина на вратата (разстояние от панта до център на тежестта)
    PB_mm = 62.3,           # разстояние от пантата до точка B (на вратата)
    PA_mm = 247.3,          # разстояние от пантата до точка A (на рамката)
    safety_factor = 1.25    # 20–30% резерв (стандартно 1.2–1.3)
    ):

    g = 9.81                          # земно ускорение
    L = length_mm / 1000.0            # в метри
    PB = PB_mm / 1000.0               # в метри
    PA = PA_mm / 1000.0               # в метри

    # Момент от тежестта на вратата спрямо пантата (ако центърът на тежестта е в средата)
    M_gravity = weight_kg * g * (L / 2) # Nm

    # Коефициент на рамото на пружината (синус на ъгъла между пружината и вратата)
    # При 90° пружината е почти перпендикулярна → sin ≈ 1
    # Най-голямото натоварване е около 50–60° → средно вземаме ~0.92–0.95
    lever_factor = 0.94 # доказана стойност за тази геометрия

    # Необходим момент от пружината
    M_required = M_gravity * safety_factor

    # Сила на пружината в Нютони
    F_newton = M_required / (PB * lever_factor)

    # Окръгляме нагоре до най-близкото 5 или 10 N (както се продават)
    F_rounded = 5 * math.ceil(F_newton / 5)

    print("═" * 60)
    print(f"Врата: {weight_kg} кг × {length_mm} мм")
    print(f"Монтаж: PB = {PB_mm:.1f} мм   |   PA = {PA_mm:.1f} мм")
    print(f"─" * 50)
    print(f"Момент от тежестта:      {M_gravity:.2f} Nm")
    print(f"Необходим момент (с резерв {safety_factor}): {M_required:.2f} Nm")
    print(f"Необходима сила:         {F_newton:.1f} N")
    print(f"Препоръчителна пружина:  {F_rounded} N")
    print(f"(или две по {F_rounded//2 if F_rounded%2==0 else F_rounded} N ако монтираш две пружини)")

    return F_rounded


print("КАЛКУЛАТОР")
print("═" * 60)

# Решаване на квадратно уравнение: 2×PB² + 2×L_closed×PB + (L_closed² - L_open²) = 0
a = 2
b = 2 * L_closed
c = L_closed**2 - L_open**2

# Дискриминанта
D = b**2 - 4*a*c

# Положителен корен
PB = (-b + math.sqrt(D)) / (2*a)
PA = PB + L_closed

print(f"   PB (на вратата) = {PB:.1f} мм")
print(f"   PA (на рамката) = {PA:.1f} мм")
print(f"   Проверка при 0° → |PA − PB| = {abs(PA - PB):.1f} мм  → искано {L_closed} мм")
print(f"   Проверка при 90° → √(PA² + PB²) = {math.sqrt(PA**2 + PB**2):.1f} мм  → искано {L_open} мм")

# Допълнителна проверка
print(f"\nДопълнителна проверка:")
print(f"   Дискриминанта D = {D:.1f}")
print(f"   Уравнение: 2×PB² + {b:.1f}×PB + {c:.1f} = 0")

gas_strut_force(weight_kg=8, length_mm=500, PB_mm=PB, PA_mm=PA)

КАЛКУЛАТОР
════════════════════════════════════════════════════════════
   PB (на вратата) = 62.3 мм
   PA (на рамката) = 247.3 мм
   Проверка при 0° → |PA − PB| = 185.0 мм  → искано 185.0 мм
   Проверка при 90° → √(PA² + PB²) = 255.0 мм  → искано 255.0 мм

Допълнителна проверка:
   Дискриминанта D = 383300.0
   Уравнение: 2×PB² + 370.0×PB + -30800.0 = 0
════════════════════════════════════════════════════════════
Врата: 8 кг × 500 мм
Монтаж: PB = 62.3 мм   |   PA = 247.3 мм
──────────────────────────────────────────────────
Момент от тежестта:      19.62 Nm
Необходим момент (с резерв 1.25): 24.53 Nm
Необходима сила:         418.9 N
Препоръчителна пружина:  420 N
(или две по 210 N ако монтираш две пружини)