conjugate(x_) := re(x)-i*im(x);