function pFactor = act(a,K,n); % activation pFactor = K.*a.^n./(1+K.*a.^n); % pFactor = (a./K).^n./(1+(a./K).^n); % alternative equation
