function spectrum = dsplaplacian(g,m)
% size of g = 3^(m+1)

M = mlaplacian(g,m) ;

badaddress = 2*ones(m+1,1) ;
badindex = indexsg(badaddress,m) ;

indices = [2:badindex-1 badindex+1:3^(m+1)-1] ;
indices = singleindices(indices,m) ;
M = M(indices,indices) ;

spectrum = eig(M)*(3/2)*5^m;
