Question: Wigner Rotation in Maple

Has anyone gotten Wigner Rotation to work in Maple. My implementation based on a formula derived in Edmonds, which is equivalent to whatever is on wikipedia, just isn't working. And I don't know where the problem lies. I know for sure that my z axis rotaition works. But I guess my Y-axis rotation is failing for some reason. I'm sure that it matches the formula. I do a Z-Y-Z rotation, so the resulting rotation matrix is real. 

Ideas? Thanks fo any help you can provide!

Please Wait...