import numpy as np
import matplotlib.pyplot as plt

# Veriler
F0 = 2098.091  # N
F1 = 8371  # N
P0 = 476000  # Pa
P1 = 5240000  # Pa
V0 = 3.3264e-3  # m³
V1 = 3.024e-5  # m³

# İzotermal sabit
iso_constant = P0 * V0

# Politropik sabit
poly_constant = P0 * (V0 ** 1.35)

# Piston alanı (m²)
A = 1.6e-3

# Vuruş (stroke) değerleri (m)
strokes = np.linspace(0, 0.163, 100)

# Hacim değerlerini hesaplayalım (m³)
volumes = V0 + strokes * (V1 - V0) / 0.163

# İzotermal basınç değerlerini hesaplayalım (Pa)
iso_pressures = iso_constant / volumes

# Politropik basınç değerlerini hesaplayalım (Pa)
poly_pressures = poly_constant / (volumes ** 1.35)

# Yük (force) değerlerini hesaplayalım (N)
iso_loads = iso_pressures * A
poly_loads = poly_pressures * A


# Grafik oluşturma
plt.plot(strokes, iso_loads, marker='o', label='Isothermal')
plt.plot(strokes, poly_loads, marker='x', label='Polytropic')
plt.title('Load-Stroke Graph')
plt.xlabel('Stroke (m)')
plt.ylabel('Load (N)')
plt.legend()
plt.grid(True)
plt.show()