Following Takahashi ("On the concrete construction of hyperbolic structure of 3-manifolds"), I was able to construct the Euclidean cusp cross-section for the 5_2 knot complement (please see http://kias.dyndns.org/topogeoimages/5_2.cusp.png), and determine the vertex invariants

x = 0.12256116687665364 + 0.7448617666197441 i
y = 0.6623589786223728 + 0.5622795120623012 i and
z = 0.7849201454990264 + 1.3071412786820453 i.

I would like to be able to compute the cusp shape from x, y and z, but as a retired physicist and learning this on my own, I have not been able to come up with an algorithm to do this. I have been able to compute the cusp shape following Yokota ("On the cusp shape of hyperbolic knots"), but this is from the knot, not from the cusp cross-section. I feel this should be easy, but the SnapPy documentation (and code, for that matter) is not entirely clear to me, nor is Coulson et. al. ("Computing Arithmetic Invariants of 3-Manifolds").

If I could just see the cusp shape in terms of x1, x2, x3, y1, y2, y3, z1, z2 and z3, I am reasonably sure I can work backwards to an algorithm I can use. Thanks in advance.