New take on Burak's result

Did an independent calculation of PK system in the slice using (1,1,0,0)
and confirmed Burak's results. Made a plot projected in a Gramm-Schmitt
orthogonal basis normal to the group tangent, but it doesn't change
things too much. Uploaded rotate-able Matlab figs.

Burak, please provide numerical root in Cartesian basis and I'll add it
to the figures.

Will comment Matlab codes for human consumption tonight but want to
commit before things get too crazy.