/* load["lib/all.k"];; */
Sweyl("x,y,z");
v=[x,y,z];
a=Sminimal(v);
b=a[0];
sm1_pmat(b);
b[1]*b[0]:
b[2]*b[1]:


/*

S^1 -- b[2]--> S^3 -- b[1] --> S^3 -- b[0] --> S^1
S^m is the set of row vectors.


*/