from pygasflow.interactive.diagrams import ShockPolarDiagram
d = ShockPolarDiagram(mach_number=4, gamma=1.4, theta=30)
d.show_figure()