name doubling-oriented XXYZZR coordinates
parameter c
assume a^2+c^2=1
variable X
variable XX
variable Y
variable Z
variable ZZ
variable R
satisfying x = X/Z
satisfying y = Y/ZZ
satisfying XX = X^2
satisfying ZZ = Z^2
satisfying R = 2 X Z