double f(double x_re, double x_im, double y_re, double y_im) {
double r1201661 = x_re;
double r1201662 = r1201661 * r1201661;
double r1201663 = x_im;
double r1201664 = r1201663 * r1201663;
double r1201665 = r1201662 + r1201664;
double r1201666 = sqrt(r1201665);
double r1201667 = log(r1201666);
double r1201668 = y_re;
double r1201669 = r1201667 * r1201668;
double r1201670 = atan2(r1201663, r1201661);
double r1201671 = y_im;
double r1201672 = r1201670 * r1201671;
double r1201673 = r1201669 - r1201672;
double r1201674 = exp(r1201673);
double r1201675 = r1201667 * r1201671;
double r1201676 = r1201670 * r1201668;
double r1201677 = r1201675 + r1201676;
double r1201678 = sin(r1201677);
double r1201679 = r1201674 * r1201678;
return r1201679;
}