da.potential = zeros(T,1);
da.activity  = zeros(T,1);
da.self_w    = zeros(1);
da.decay     = 150;
da.noise     = 0;
da.baseline  = 0;
da.th_out    = 0;

da.epsilon   = 1;
da.lambda    = 0;

%da.inp       = zeros(T,1);