s=[  1     1  1    1  857;
     1     2   2   2  951;
     1     3   3   3  909;
     2     1   2   3  878;
     2     2   3   1  973;
     2     3   1   2  899;
     3     1   3   2  803;
     3     2   1   3  1030;
     3     3   2   1   927]; 
[result,sum0]=zjjc(s,1)