I have a college project in which I need to create a program which simulates an Erlang(k) variable. I am not talking about the Erlang language, it's about a standard Gamma(0,1,k), k ϵ N+. I appreciate if someone can help me out or give me advice on how to approach this problem.
I've written the formulas and theorem that comes with this problem: Erlang distribution equations
All the best,