Mathematica
In week 2 it was decided that it would be useful to create a function for NDSolve, so that different data sets can easily be computed and verified.
This page will hold relevant data on the issue.
Mathematica on defining functions.
The idea of the work is to model the ODE systems in a more efficient way. The theory is that the blue region on the graph can be modeled using an exponential, and the red region using a different linear system. This would make the calculations many times faster (and potentially more accurate than a floating point general ODE solver).
The yellow highlighted region would be modeled again in a different way due to its highly dynamic behaviour. When plotted backwards it is possible to determine the threshold of voltage at which a spike is inevitable. This looks to be largely independent of other variable (around -42mV)