double f(double x_re, double x_im, double y_re, double y_im) {
double r1308882 = x_re;
double r1308883 = r1308882 * r1308882;
double r1308884 = x_im;
double r1308885 = r1308884 * r1308884;
double r1308886 = r1308883 + r1308885;
double r1308887 = sqrt(r1308886);
double r1308888 = log(r1308887);
double r1308889 = y_re;
double r1308890 = r1308888 * r1308889;
double r1308891 = atan2(r1308884, r1308882);
double r1308892 = y_im;
double r1308893 = r1308891 * r1308892;
double r1308894 = r1308890 - r1308893;
double r1308895 = exp(r1308894);
double r1308896 = r1308888 * r1308892;
double r1308897 = r1308891 * r1308889;
double r1308898 = r1308896 + r1308897;
double r1308899 = sin(r1308898);
double r1308900 = r1308895 * r1308899;
return r1308900;
}