first run current_analysis.hoc to generate the data including current trajectory for Matlab code.
then run current_analysis.m