y = [ 551, 595, 639, 417, 563; 457, 580, 615, 449, 631; 
   450, 508, 511, 517, 522; 731, 583, 573, 438, 613;
   499, 633, 648, 415, 656; 632, 517, 677, 555, 679]

anova1(y)