function q = q_fsa(P,params)



q0  = params.q0;

P_A = params.E_A/params.tauA;

q = q0./(1+P/P_A);

