On a stochastic approach to model the double phosphorylation/dephosphorylation cycle