computeEntries_F_matrices(FList, TList, numGrp, numVar, index1, index2, const1, const2, const3, const4, const5)

function computeEntries_F_matrices(FList, TList, numGrp, numVar, index1, index2, const1, const2, const3, const4, const5)

m=index1;
j=index2;

c=const1;
r=const2;
s=const3;
t1=const4;
t2=const5;

for i=1:numGrp
	FList{i}(m,m)=c*(TList{i}(1,1)+2*t1*TList{i}(1,2)+t2*TList{i}(2,2));
	FList{i}(j,j)=c*(t2*TList{i}(1,1)-2*t1*TList{i}(1,2)+TList{i}(2,2));  
	FList{i}(m,j)=c*(t1*(TList{i}(2,2)-TList{i}(1,1))+(1-t2)*TList{i}(1,2));
	FList{i}(j,m)=FList{i}(m,j);
    
    for l=1:numVar
	    if (l ~= m )&(l ~=j) 
		    b1=FList{i}(l,m);
		    b2=FList{i}(l,j);
            
		    FList{i}(l,m) = b1 +  s * (b2 - r * b1);
			FList{i}(l,j) = b2 -  s * (b1 + r * b2);
     		FList{i}(m,l) = FList{i}(l,m);
            FList{i}(j,l) = FList{i}(l,j);
        end;
    end;
end