f_u=function(ub,vb,u,v){ # # Function f_u computes the rate of binding/unbinding for u # # Return rate function return(c(-u*(ub-u)+v*(vb-v))); }