function f=rtf1(R,T)

for j=1:6
fr(j)=sum(R(j,1,1:R(j,1,13)))/R(j,1,13)-sum(R(j,2,1:R(j,2,13)))/R(j,2,13)-sum(R(j,3,1:R(j,3,13)))/R(j,3,13)+sum(R(j,4,1:R(j,4,13)))/R(j,4,13);
ft(j)=sum(T(j,1,1:R(j,1,13)))/R(j,1,13)-sum(T(j,2,1:R(j,2,13)))/R(j,2,13)-sum(T(j,3,1:R(j,3,13)))/R(j,3,13)+sum(T(j,4,1:R(j,4,13)))/R(j,4,13);

end

f=mean(fr)/mean(ft);