謝謝Neisseria大介紹globals函數讓我解決了幫球自動取名的問題
但又碰上新的問題，
我希望在每一個迴圈中，自動把每個球的位置塞入新的list中，
好讓我可以每個迴圈重新計算球與球間的距離。
若不用迴圈我會這樣寫:
b_new_pos_list = []
b_new_pos_list.append(ball_0.pos)
b_new_pos_list.append(ball_1.pos)
b_new_pos_list.append(ball_2.pos)
b_new_pos_list.append(ball_3.pos)
相同的，我想利用for loop自動把每個球的位置放入b_new_pos_list中
因此我嘗試這樣寫。
b_new_pos_list=[]
for N in range(0,4,1):
b_new_pos_list.append(ball_N.pos)
但當然還是不行，因為系統沒辦法自動判斷出ball_N.pos指的就是
ball_0~3.pos
再次感謝。

跑模擬用普通python list跑會很慢喔...你的目的應該是N-body simulation 有很多資源可以找隨便google了一個 https://goo.gl/4ZZmaL裡面有討論基本N個球在重力影響下計算軌跡的例子先參考看看

好的，感謝。不好意思，看完painkiller大的例子，還是不知道如何解決我的問題 0rz..我大概有概念了，謝謝painkiller大

Vpython?