clear
clf
%?membrane?constants
tau_m?=?0.030;
Rm?=?3e7;
thresh=0.05;
t_spikes=[];
No_spikes=0;
in_abs_ref=0;
abs_ref_count=0;
%?time
delta_t?=?0.001;?%?integration?time?step
duration?=?0.5;%?duration?of?the?simulation
abs_T=0.003;%absolute?refractory?period絕對不應期
abs_ref_steps=round（abs_T./delta_t）;%四舍五入
No_steps?=?round（duration./delta_t）;
times?=?linspace（0?duration?No_steps?+?1）;?%?the?time?points?for?update
%?membrane?potential
V_0=?0.03;
V?=?zeros（1?No_steps?+1）;
V（1）=?V_0;

%%?injection?current
I_0?=2e-9;
I_duration?=?0.25;
I_start=?0.05;
I_start_index=?round（I_start./?delta_t）;
I_finish_index?=?round（（I_start+I_duration）./delta_t）;
Iinj=zeros（1No_steps+1）;
Iinj（I_start_index:?I_finish_index）=?I_0;
%%?noise?current
std_dev_noise?=?1.5e-10;
I_noise