A=[1 5 6 8
1/5 1 3 5
1/6 1/3 1 4
1/8 1/5 1/4 1];
[a,b]=eig(A);
maxeignvalue=max(max(b))
index=find(b==max(max(b)));
eigenvector=a(:,index)
maxeignvalue =
4.2627
eigenvector =
0.9349
0.3088
0.1616
0.0671