#!/usr/bin/env python3
#
import matplotlib.pyplot as plt
import numpy as np

xpoints = np.array([100, 95, 90, 85, 80, 75, 70, 65, 60, 55, 50, 45, 40, 35, 30, 25, 20, 15, 10, 5])
PHXpoints = np.array([60.84, 56.64, 49.49, 44.92, 39.31, 35.83, 31.81, 27.88, 25.67, 23.15, 21.09, 17.44, 14.36, 11.47, 10.08, 8.32, 5.79, 3.25, 2.05, 0.73])
ATLpoints = np.array([41.02, 36.85, 33.89, 32.02, 29.50, 27.99, 26.82, 25.08, 23.78, 22.50, 20.92, 19.42, 17.83, 16.69, 15.42, 13.16, 10.25, 8.08, 6.34, 4.84])
MDWpoints = np.array([58.83, 54.33, 49.19, 46.06, 43.64, 40.71, 38.50, 36.27, 34.01, 31.66, 29.57, 26.53, 23.61, 21.56, 18.77, 16.89, 13.83, 10.82, 7.54, 3.35])
BWIpoints = np.array([56.99, 52.77, 48.89, 45.93, 42.43, 39.00, 36.77, 34.41, 31.47, 29.28, 26.93, 24.79, 22.58, 20.74, 17.69, 15.48, 12.16, 9.44, 6.86, 4.15])
DALpoints = np.array([56.42, 49.84, 46.90, 43.77, 40.87, 37.95, 34.76, 33.09, 30.31, 27.12, 24.27, 21.80, 17.96, 14.93, 12.24, 9.76, 7.80, 5.73, 3.67, 1.51])


plt.plot(xpoints, PHXpoints)
plt.plot(xpoints, ATLpoints)
plt.plot(xpoints, MDWpoints)
plt.plot(xpoints, BWIpoints)
plt.plot(xpoints, DALpoints)

plt.xlabel("Domicile Sen")
plt.ylabel("System Sen Needed to Hold")

ax = plt.gca()
ax.invert_yaxis()
ax.invert_xaxis()

plt.show()