function [firingMtx]=pfcSim(inputSig,attnSig,weightMtx)

exec('sigmoid.sci');
sigAlpha=1.5;
sigBeta=3;

if length(inputSig)~=length(attnSig) then
  error('first two inputs not the same length!');
end;

firingMtx=zeros(size(weightMtx,1),length(inputSig));

for k=2:size(firingMtx,2),
 firingMtx(:,k)=sigmoid(weightMtx*firingMtx(:,k-1)+inputSig(k)+attnSig(k),sigAlpha,sigBeta);
end;