X <- matrix(c( 1,1,3,2,3,2,1,1,1,
    1,1,2,2,2,3,3,3,2,1,3,2,
    4,3,5,5,5,5),9,3,byrow=T)
Y <- matrix(c( 4,4,-1.07846,
               3,3,1.214359,
               2,2,0.307180,
               2,3,-0.385641,
               2,1,-0.078461,
               1,1,1.61436,
               1,2,0.814359,
               2,1,-0.0641016,
       1,2,1.535900 ),9,3,byrow=T)