B = (3, 1.24)
A = (0, 0)
AB = Vector(A, B)
a = 3 (slider from -5 to 5, step 0.1)
b = 3 (slider from -5 to 5, step 0.1)
k = 20 (slider from -20 to 20, step 0.1)
ABe2: AB * AB
FBAi: k * a * b / ABe2
ABu = UnitVector(AB)
FBA: FBAi * ABu
C = (6, 0)
c = 3 (slider from -5 to 5, step 0.1)
CB = Vector(C, B)
CBu = UnitVector(CB)
FBCi: k * b * c / (CB * CB)
FBC: FBCi * CBu
D = (3, 5.16)
d = 3 (slider from -5 to 5, step 0.1)
DB = Vector(D, B)
DBe2: DB * DB
DBu = UnitVector(DB)
FBDi: k * b * d / DBe2
FBD: FBDi * DBu
FB: FBA + FBC + FBD
poly3, e, f, g, F = Polygon(A, C, 3.0)
E = (4.46, 0.86)
G = (2.98, 3.44)
poly1, h, i, j, H = Polygon(E, G, 3.0)