function [t_syn, g_syn] = alpha_syn(t_peak,g_max)
%
%[t_syn, g_syn] = alpha_syn(t_peak,g_max)
%
%alpha_syn generates an alpha function with a sampling step of 
%0.01 msec. The peak time is in msec and the conductance in nS. 
%

%Compute only up to 8 times the peak value since by 8 times 
%t_peak the value of the alpha function is already less than
%one percent of its peak value
t_syn = 0:0.01:ceil(8*t_peak);
g_syn = g_max*(exp(1)/t_peak) * t_syn.*exp(-t_syn/t_peak);