Lignes de champs magnétiques
Author: Daniel Mentrard

d = 1.3 (slider from -5 to 5, step 0.1)
e = 0.5 (slider from -5 to 5, step 0.1)
m = 7.8 (slider from -5 to 50, step 0.1)
B = (6.674854093165036, 0.12910534437467686)
A = (-5.589598733699832, 0.2259289979289246)
b = Segment(B, A)
a = LineBisector(B, A)
C = Intersect(a, b)
f(x, y) = (x - x(A))^(2) + (y - y(A))^(2)
g(x, y) = (x - x(B))^(2) + (y - y(B))^(2)
liste3 = Sequence(f(x, y) * g(x, y) = i * f(x, y)^(2), i, 0.1, m, e)
liste4 = Sequence(f(x, y) * g(x, y) = i * g(x, y)^(2), i, 0.1, m, e)
liste1 = Sequence(Translate[C, Vector[i^(2) / d * (0, 1)]], i, -m, m, 2 * e)
C' = Translate(C, (0, 1))
liste2 = Sequence(CircumcircleArc[A, Element[liste1, i], B], i, 1, Length[liste1])
liste5 = Sequence(Translate[C, Vector[i^(2) / d * (0, -1)]], i, -m, m, 2 * e)
liste6 = Sequence(CircumcircleArc[A, Element[liste5, i], B], i, 1, Length[liste5])