P0 on E                                                      identity: True  (denominator nonzero: True)
P1 on E                                                      identity: True  (denominator nonzero: True)
side P0P1 tangent to C                                       identity: True  (denominator nonzero: True)
side P1P2 tangent to C                                       identity: True  (denominator nonzero: True)
reflection law at P1 (squared cosines)                       identity: True  (denominator nonzero: True)
inversion of P0 in C: midpoint = classical (x)               identity: True  (denominator nonzero: True)
inversion of P0 in C: midpoint = classical (y)               identity: True  (denominator nonzero: True)
inversion of P'0 in E: midpoint = classical (x)              identity: True  (denominator nonzero: True)
inversion of P'0 in E: midpoint = classical (y)              identity: True  (denominator nonzero: True)
inversion of P1 in C: midpoint = classical (x)               identity: True  (denominator nonzero: True)
inversion of P1 in C: midpoint = classical (y)               identity: True  (denominator nonzero: True)
inversion of P'1 in E: midpoint = classical (x)              identity: True  (denominator nonzero: True)
inversion of P'1 in E: midpoint = classical (y)              identity: True  (denominator nonzero: True)
inversion of P2 in C: midpoint = classical (x)               identity: True  (denominator nonzero: True)
inversion of P2 in C: midpoint = classical (y)               identity: True  (denominator nonzero: True)
inversion of P'2 in E: midpoint = classical (x)              identity: True  (denominator nonzero: True)
inversion of P'2 in E: midpoint = classical (y)              identity: True  (denominator nonzero: True)
inversion of P3 in C: midpoint = classical (x)               identity: True  (denominator nonzero: True)
inversion of P3 in C: midpoint = classical (y)               identity: True  (denominator nonzero: True)
inversion of P'3 in E: midpoint = classical (x)              identity: True  (denominator nonzero: True)
inversion of P'3 in E: midpoint = classical (y)              identity: True  (denominator nonzero: True)
k110 = A*A2 = 8a^4b^4/(a^2+b^2)^2                            identity for every member: True  (denominator nonzero: True)
k807 = A*Aot = 4a^4b^4/(a^2+b^2)^2                           identity for every member: True  (denominator nonzero: True)
k809 = A1*Aom = 4a^2b^2                                      identity for every member: True  (denominator nonzero: True)
k106 = A*A1 = 8a^2b^2                                        identity for every member: True  (denominator nonzero: True)
k113 = A1/A2 = (a^2+b^2)^2/(a^2b^2)                          identity for every member: True  (denominator nonzero: True)
Aot = A2/2 (Varignon)                                        identity for every member: True  (denominator nonzero: True)
Aom = A/2 (Varignon)                                         identity for every member: True  (denominator nonzero: True)
WRONG: k110 = 2a^4b^4/(a^2+b^2)^2 (GR-SI printed value)      identity for every member: False  (should be False)
WRONG: k807 = 2a^4b^4/(a^2+b^2)^2                            identity for every member: False  (should be False)
WRONG: k111 = A1*A2 constant (vs. symmetric member value 16 a^4 b^4/(a^2+b^2)^2... test vs 0 shift) identity for every member: False  (should be False)
rhombus member:   A'A'' = 16*a**4*b**4/(a**4 + 2*a**2*b**2 + b**4) ;  A A'' = 8*a**4*b**4/(a**4 + 2*a**2*b**2 + b**4)
rectangle member: A'A'' = 4*a**2*b**2 ;  A A'' = 8*a**4*b**4/(a**2 + b**2)**2
k111 differs between the two members: True  difference = -4*a**2*b**2*(a - b)**2*(a + b)**2/(a**2 + b**2)**2
