Finished code for figure4. Switched convention so that cathodic voltage values are negative, producing positive MDF values at center node. Added error in computeVoltages for infinite voltage values (due to electrode and node at same location).