formula (10) applies to ``sin(p* pi / q)`` for ``p=1,..,q-1``