B=[23,15,26,13,21,25,20,21,16,18,21,17,16,24,27,14,17,19,20,24;
    28,22,25,19,26,30,26,26,20,28,19,24,19,25,29,17,21,18,26,23;
    18,10,12,22,13,15,21,22,14,12,23,25,19,13,22,16,12,23,22,19];
anova2(B',5);

C=[];
for i=1:3
    for j=1:4
        C(i,j)=mean(B(i,[1:5]+(j-1)*5));
    end, end
C=[C; mean(C)];  C=[C mean(C')']