x1= -2 ln(u1) cos(2 u2) +