n=10
k=3
# count standard deviation
A<-sd(c(73,65,27,71,43,79,78,65,63,58))
B<-sd(c(86,72,34,79,45,86,81,66,72,67))
C<-sd(c(82,74,36,71,50,92,82,62,70,62))


#count mean

a <- mean(c(73,65,27,71,43,79,78,65,63,58))
b <- mean(c(86,72,34,79,45,86,81,66,72,67))
c <- mean(c(82,74,36,71,50,92,82,62,70,62))


#count mean between sample
D <- mean(c(a,b,c))

#count s.d between sample
E<-sd(c(a,b,c))

#count variance between samples
f = n*(E)*(E)

#count variance within samples
G = ((A*A)+(B*B)+(C*C))/k

#calculate test statistic
H = f/G

#numerator&denominator
num = k-1
den = k*(n-1)

#count F crtitical value
qf(.95,df1 = num,df2 = den)