[0,1)を生成する乱数源rand()があったとして、
確率Kの事象は 式rand() < Kで生成できる。
あとはKを加速度的に増やしていけばOK.
力学的には加速度はa=d²p/dt²なので、
K=a*t²+b*t+cで導出するのが筋かな？