Consider this pseudo-code (partly from python, you'll get the idea, I think):

Code :
A44 = matrix([\
     [  0.96126169,   0.,           0.27563736,          0],\
     [  0.07597595,   0.96126169,  -0.26495963,          0],\
     [ -0.26495963,   0.27563736,   0.92402405,          0],\
     [  0,            0,            0,                   1])
 
glMultMatrixf(Anew)
// Drawing some stuff works ok
model = glGetDoublev(GL_MODELVIEW_MATRIX)
print "model = \n", model
 
// This prints the following out to screen:
//model = 
//[[ 1.  0.  0.  0.]
// [ 0.  1.  0.  0.]
// [ 0.  0.  1.  0.]
// [ 0.  0.  0.  1.]]
// But why is it still identity matrix, like it's unaffected by multiplication???

And then I have another problem:
Code :
A44 = matrix([\
     [  0.96126169,   0.,           0.27563736,          0],\
     [  0.07597595,   0.96126169,  -0.26495963,          0],\
     [ -0.26495963,   0.27563736,   0.92402405,          0],\
     [  0,            0,            0,                   1])
glLoadMatrixd(Anew)
// Trying to draw some new things - don't work - nothing happens!
 
// And once again, when I try to print out current modelview matrix, I get 4x4 identity - why???

In the last case, I suspect that the A44 matrix I've given isn't precise enough and then I suspect that the matrix isn't orthogonal - hence nothing is drawn. Is that guess correct (that I need more digits to A44) ???