How do you simulate an Erlang(k) variable using a programming language?

Question

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 in either C++ or Java.

I've written the formulas and theorem that comes with this problem: link

Answer 1

Try this: http://www.ee.ucl.ac.uk/~mflanaga/java/Stat.html