Graphics 800,600,16,2
SetBuffer BackBuffer()
SeedRnd MilliSecs()
grav#=0.1
Type ball
 Field id
 Field x#,y#
 Field sx#,sy#
 Field rad#
End Type

For i=0 To 9
b.ball=New ball
b\id=i
b\x#=Rand(0,500)
b\y#=Rand(0,500)
b\rad#=10
Next

t=CreateTimer(60)
While Not KeyHit(1)
Cls

For b.ball=Each ball
 Plot b\x#+b\rad#,b\y#+b\rad#
 Oval b\x#,b\y#, b\rad#*2,b\rad#*2,0
 b\x#=b\x#+b\sx#
 b\y#=b\y#+b\sy#
 b\sx#=b\sx#
 b\sy#=b\sy#+grav#

 If b\y>=GraphicsHeight()-b\rad#*2 Then b\y=GraphicsHeight()-b\rad#*2:b\sy=0

 For b1.ball=Each ball
d=Sqr(((b1\x#+b1\rad#)-(b\x#+b\rad#))*((b1\x#+b1\rad#)-(b\x#+b\rad#))+((b1\y#+b1\rad#)-(b\y#+b\rad#))*((b1\y#+b1\rad#)-(b\y#+b\rad#)))
   

  If d<20 And b\id<>b1\id

d=Sqr(((b1\x#+b1\rad#)-(b\x#+b\rad#))*((b1\x#+b1\rad#)-(b\x#+b\rad#))+((b1\y#+b1\rad#)-(b\y#+b\rad#))*((b1\y#+b1\rad#)-(b\y#+b\rad#)))
   
    alpha = ASin(Abs(b\x-b1\x)/d)

   
  EndIf


 Next
Next




WaitTimer(t)
Flip
Wend

Local Treff_winkel = Winkel_2D#(k\X, k\Y, k2\X, k2\Y)            
k\w = Treff_winkel
k2\w = -Treff_winkel

[...]

Function Winkel_2D#(x1#, y1#, x2#, y2#)
   Return (ATan2(y1 - y2, x1 - x2) + 270) Mod 360
End Function