j1
j1 = x->input[1];
x1 = j1;
x1 = PLUS(x1,j1);
u32 j0, j1, j2, j3, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15;
j1 = x->input[1];
x1 = j1;
x1 = PLUS(x1,j1);
u32 j0, j1, j2, j3, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15;
return(x*(u/v) + tpi*(j1(x)*log(x)-one/x));
b = j1(x);
w = j1(x);
if(n==1) return(j1(x));
double j1(double);
u32 j0, j1, j2, j3, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15;
j1 = x->input[1];
x1 = j1;
x1 = PLUS(x1,j1);
int i, j, m, i0, i1, j0, j1, nc;
j1 = J[i1 + 1] - 1;
J[i1] = j1;
for (i = i0, j = j0; i <= i1 && j <= j1; i++, j++)
while (j <= j1) {
change(file1, f1, file2, f2, i0, i1, j0, j1, &flags);
j1 = J[i1 - 1] + 1;
J[i1] = j1;
change(file1, f1, file2, f2, i1, i0, j1, j0, &flags);