  function [L] = mylaplacian(A)
% function [L] = laplacian(A) 
% converts A to a Laplacian L
    L = spones(A);
    L = full(L);
    L = L - diag(diag(L));    
    L = spdiags(sum(L,2),0,-L); 
%

