function act = act(in_conc, k, m)
% in_conc = input concentration of state a which activates state b
% k = threshold for regulatory influence
% m = steepness parameter
act = ((in_conc*k)^m)/(1+(in_conc*k)^m);