function y = discharge_rate(x);
% estimate discharge rate of spike time train

x = diff(x);
x = x.*(x>0);
x = nonzeros(x)';

T = sum(x);
y = size(x, 2) / T;