Tsai, I forgot to mention that Sb should be multiplied by the number of pcitures per class (in this case 2), you are right on that too. Refer to section 2.4 of the following paper to note the error: http://www.face-rec.org/algorithms/LDA/belhumeur96eigenfaces.pdf
If you guys would like to look at another sample Matlab code for fisher faces, look at: