Скрипт моделирует упругое взаимодействие частиц разной массы.

Разрешение дисплея смартфона-240x320 pixels.

import appuifw,e32
from key_codes import*
from random import choice
from math import hypot

appuifw.app.screen='full'
c=appuifw.Canvas()
appuifw.app.body=c
c.clear(0x000000)

class M:
def __init__(self,x,y,u,v,r):
self.x,self.y=x,y
self.u,self.v=u,v
self.c,self.r=c,r
def dryak(self):
self.c.ellipse((self.x-self.r,\
self.y-self.r,self.x+self.r,self.y+\
self.r),fill=0xff0000)

X,Y,U,V,Un,Ut,R=[],[],[],[],[],[],[]
for i in range(21):
X.append(choice(range(15,215)))
Y.append(choice(range(15,305)))
U.append(choice(range(1,9)))
V.append(choice(range(1,9)))
R.append(choice(range(7,27)))

while True:
for i in range(21):
for j in range(21):
if j>i:
a,b=X[j]-X[i],Y[j]-Y[i]
d=hypot(a,b)
if d<=R[i]+R[j]:
Uni=(U[i]*a+V[i]*b)/d
Uti=(V[i]*a-U[i]*b)/d
Unj=(U[j]*a+V[j]*b)/d
Utj=(V[j]*a-U[j]*b)/d
mi,mj=R[i]**2,R[j]**2
u0=(2*mj*Unj+mi*Uni-mj*\
Uni)/(mi+mj)
u1=(2*mi*Uni+mj*Unj-mi*\
Unj)/(mi+mj)
Uni,Unj=u0,u1
U[i]=(a*Uni-b*Uti)/d
V[i]=(a*Uti+b*Uni)/d
U[j]=(a*Unj-b*Utj)/d
V[j]=(a*Utj+b*Unj)/d
X[i]+=U[i]
Y[i]+=V[i]
X[j]+=U[j]
Y[j]+=V[j]

for i in range(21):
if X[i]-R[i]<=0 or X[i]+R[i]>=\
240:
U[i]=-U[i]
if Y[i]-R[i]<=0 or Y[i]+R[i]>=\
320:
V[i]=-V[i]
X[i]+=U[i]
Y[i]+=V[i]
for i in range(21):
M(X[i],Y[i],U[i],V[i],R[i]).dryak()
e32.ao_sleep(0.0001)
c.clear(0x000000)