|p||One of the real diagonals corresponding to the prescribed offdiagonal.|
|q||The offdiagonal in question.|
|r||The other real diagonal corresponding to the prescribed offdiagonal.|
|return value||Rotation matrix zeroing out q.|
For a matrix M with real diagonal, there is a rotation S such that the symmetic SMS' has a prescribed zero offdiagonal.
M = [1,2,3;4,5,6;7,8,9] M = M+M'; col = 2 row = 1 J = eye(numRows(M)); J[col,row;col,row] = jacobi(M[row;row], M[row;col],M[col;col]); J*M*J'