function G=carbon
    G=graph(bucky);
    %G=rmnode(G,[31:60]);
end