In reply to Wei Chen.

Wei, I agree with you that SAS/IML is a great product -- sort of a SAS/STAT on steroids for matrix operations. My colleague Rick Wicklin's blog The DO Loop offers hundreds of examples of things you can do in SAS/IML.