How to simulate Erlang B traffic pattern

Hello everybody. I am a student from HK. May I ask how can I write a program that simulate a cellular environment with Erlang B traffic pattern?

For example, how can I write a Matlab program that can, say, simulate N mobiles, each generate calls and have differnet call hold time(service time) which satisfies Erlang B traffic pattern?

Thank you very much in advance.