Вот переменные и сама формула (по сути)
m1=6**20
m2=6**20
x=400
y=250
xx=200
yy=200
G=6,67430*(10**-11) #Гравитационная постоянная
r
(x-xx)**2+(y-yy)**2)**0.5 #Подсчет расстояния по координатамr=round(r,0) #Округление что-бы было проще
F=G*((m1*m2)/r**2) #Собственно на этом этапе ошибка