Use nid, eid, sid and rid as node, element, support and reaction counter and ID names respectively. Also more simplifications, improved aesthetics of tables, code comments.
Use speye functions instead of sparse(eye()) combination.