all (all (abs (sqrt ([4, -4; i, 1-i])
	       - [2, 2i; exp (0.5 * log (i)), exp (0.5 * log (1-i))])
	  < sqrt (eps)))
