Simulasi yang berkaitan dengan gugus bola sangat penting dilakukan untuk mengetahui evolusi dan dinamika gugus bola. Dalam penelitian ini, kami mendeskripsikan sebuah model gugus bola dengan densitas Plummer isotropik melalui teknik sampling rejection yang mirip dengan kali code. Kemudian, kami membuat modul program dalam bahasa pemrograman Python. Model yang dihasilkan oleh modul ini sesuai dengan beberapa parameter teoretis Plummer. Walaupun model Plummer ini hanya cocok untuk daerah di sekitar inti dari gugus bola hasil observasi, tetapi model Plummer distribusi bintang-bintang dapat digunakan sebagai kondisi awal gugus bola dalam simulasi gugus bola.