function [A,indices,rooturls]=loadSBMatrix;

truesize=683446;

load stanford-berkeley-bool-sorted.dat;
A=spconvert(stanford_berkeley_bool_sorted);
% make the matrix square
n=max(size(A));
A(n,n)=0;
A=A(1:truesize,1:truesize);
% normalize the rows to sum to 1
A=normalize(A,1);
% take the transpose A^T (which is what we multiply by)
A=A';
load stanford-berkeley-sorted-roots.dat;
indices=stanford_berkeley_sorted_roots;
indices=indices(find(indices<truesize));
rooturls=textread('rooturls.txt','%s');
rooturls=rooturls(1:max(size(indices)));