(* Simulation: To Estimate the Probability Of Hitting M Before Zero <<Knox *)



AbsorptionPctAtM[p_, n_, M0_, numreps_] :=

 Module[{absatM, empprob, poslist, lastpos, i},

  		absatM = 0;

  		Do[poslist = SimulatedRandomWalk[p, n, M0];

   			lastpos = Last[poslist];

   If[lastpos == M0, absatM = absatM + 1, Null], {i, 1, numreps}];

  empprob = N[absatM/numreps]]







SeedRandom[342197]

AbsorptionPctAtM[.4, 3, 5, 100]

AbsorptionPctAtM[.4, 3, 5, 100]