Now can simulate multiple neurons with gap junctions in vclamp/iclamp mode and solved using RK/CVode classes