clear all
clc
M=[7 7 15 11 9;12 17 12 18 18; 14 18 18 19 19;19 25 22 19 23;7 10 11 15 11]
yb=mean(M,2)
ybb=mean(yb)
SStotal=sum(sum((M-ybb).^2))
MSStotal=SStotal/24
SStreat=sum((yb-ybb).^2)*5
MSStreat=SStreat/4
for j=1:5
    Myb(:,j)=mean(M,2);
end
SSerro=sum(sum(((M-Myb).^2)))
MSSerro=SSerro/(25-5)
F=MSStreat/MSSerro