function [GMatrix,BiggestEigenvector]=GMatrixMaker(AMatrix,q)
n=length(AMatrix);
G1=(q/n).*ones(n,n);
G3=sum(AMatrix').^-1;
G3=repmat(G3,n,1);
G2=(AMatrix.').*((1-q).*G3);
GMatrix=G1+G2;
[EigenVectors,EigenValues]=eig(GMatrix);
BiggestEigenvector=EigenVectors(:,1);
BiggestEigenvector=BiggestEigenvector/sum(BiggestEigenvector);
