begin # RCohen(), file a block, 2005 March 08, Jim Worthey function RCohen(A) begin SmallMatrix = A'*A SM1 = inv(SmallMatrix) R = A*SM1*A' return R end # fn RCohen() end # file a block