run PRC_script_curr.hoc first and then run Fig2B.m