Generic Givens Rotation Matrix in EQN

[Damian McGuckin]

Would anybody have a version of this in EQN. I do not mean the simple 2x2 
matrix but the generic nxn square matrix

        G(i, j, theta)

which is, when i > j,

        G sub kk = 1  ... k != i, j
        G sub ii = G sub jj = cos(theta)
        G sub ji = - G sub ij = - sin(theta)

See:
        https://en.wikipedia.org/wiki/Givens_rotation

Thanks heaps - Damian

Pacific Engineering Systems International, 277-279 Broadway, Glebe NSW 2037
Ph:+61-2-8571-0847 .. Fx:+61-2-9692-9623 | unsolicited email not wanted here
Views & opinions here are mine and not those of any past or present employer