Fit circle to shock to a particular beta value for many inlinations.

Example output is included
plot the inner shock with the circular fit. Taking theta_c as a free parameter

computes other conics.