cis(x_) := cos(x)+i*sin(x);