user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<cdflib_beta_frac>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  mov    %r8,%r13
  push   %r12
  mov    %r9,%r12
  push   %rbp
  mov    %rcx,%rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x38,%rsp
  call   <cdflib_beta_cdflib_rcomp@plt>
  pxor   %xmm7,%xmm7
  ucomisd %xmm7,%xmm0
  movsd  %xmm0,0x28(%rsp)
  jp     L42
  jne    L42
L20:
  movsd  0x28(%rsp),%xmm0
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L42:
  mov    (%r15),%rdx
  movsd  D5(%rip),%xmm7
  movsd  (%r14),%xmm10
  movsd  0x0(%r13),%xmm4
  movq   %rdx,%xmm6
  movapd %xmm7,%xmm13
  movsd  (%r12),%xmm5
  movsd  (%rbx),%xmm11
  divsd  %xmm10,%xmm6
  addsd  %xmm7,%xmm4
  movapd %xmm7,%xmm14
  movsd  %xmm5,0x8(%rsp)
  movapd %xmm7,%xmm8
  movapd %xmm7,%xmm5
  movapd %xmm4,%xmm9
  movsd  %xmm4,0x10(%rsp)
  divsd  %xmm10,%xmm13
  movsd  %xmm6,0x18(%rsp)
  movsd  0x0(%rbp),%xmm6
  addsd  %xmm7,%xmm6
  movsd  %xmm6,0x20(%rsp)
  movapd %xmm10,%xmm6
  addsd  %xmm7,%xmm6
  addsd  %xmm7,%xmm13
  divsd  %xmm13,%xmm9
  movapd %xmm13,%xmm12
  divsd  %xmm4,%xmm12
  pxor   %xmm4,%xmm4
  movapd %xmm4,%xmm1
  jmp    L49
  nopl   0x0(%rax)
L30:
  movapd %xmm8,%xmm1
  movapd %xmm9,%xmm14
  movapd %xmm0,%xmm8
  divsd  %xmm3,%xmm1
  movapd %xmm0,%xmm12
  movapd %xmm7,%xmm9
  divsd  %xmm3,%xmm14
L49:
  movapd %xmm10,%xmm0
  addsd  %xmm7,%xmm4
  movsd  0x18(%rsp),%xmm2
  movq   %rdx,%xmm3
  addsd  %xmm5,%xmm2
  movapd %xmm4,%xmm15
  subsd  %xmm4,%xmm3
  mulsd  %xmm5,%xmm2
  mulsd  %xmm4,%xmm3
  mulsd  %xmm11,%xmm3
  divsd  %xmm6,%xmm0
  divsd  %xmm10,%xmm15
  mulsd  %xmm0,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm11,%xmm2
  mulsd  %xmm3,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm13,%xmm2
  divsd  %xmm6,%xmm3
  addsd  D7(%rip),%xmm6
  addsd  %xmm15,%xmm2
  movapd %xmm15,%xmm5
  addsd  %xmm7,%xmm5
  addsd  %xmm15,%xmm2
  movapd %xmm5,%xmm15
  divsd  %xmm2,%xmm15
  movsd  0x20(%rsp),%xmm2
  mulsd  %xmm4,%xmm2
  addsd  0x10(%rsp),%xmm2
  mulsd  %xmm15,%xmm2
  movapd %xmm3,%xmm15
  movapd %xmm14,%xmm3
  addsd  %xmm4,%xmm15
  mulsd  %xmm0,%xmm3
  mulsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm15
  movapd %xmm9,%xmm2
  mulsd  %xmm15,%xmm2
  mulsd  %xmm8,%xmm15
  addsd  %xmm2,%xmm3
  movsd  0x8(%rsp),%xmm2
  addsd  %xmm15,%xmm0
  divsd  %xmm3,%xmm0
  mulsd  %xmm0,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm12,%xmm1
  andpd  D20(%rip),%xmm1
  comisd %xmm1,%xmm2
  jb     L30
  mulsd  0x28(%rsp),%xmm0
  movsd  %xmm0,0x28(%rsp)
  jmp    L20
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D5
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D7
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D20
data_type: qword[2]

tool: D5: 1

tool: D7: 2

tool: {"D20": [0x7fffffffffffffff,0x0000000000000000]}

assistant: ```c
double cdflib_beta_frac(const double *a, const double *b, const double *x, const double *y, const double *lambda, const double *eps) {
  double bfrac, alpha, an, anp1, beta, bn, bnp1, c, c0, c1, e, n, p, r, r0, s, t, w, yp1;
  bfrac = cdflib_beta_cdflib_rcomp(a, b, x, y);
  if (bfrac == 0.0e0) {
    return bfrac;
  }
  c = 1.0e0 + *lambda;
  c0 = *b / *a;
  c1 = 1.0e0 + 1.0e0 / *a;
  yp1 = *y + 1.0e0;
  n = 0.0e0;
  p = 1.0e0;
  s = *a + 1.0e0;
  an = 0.0e0;
  bn = anp1 = 1.0e0;
  bnp1 = c / c1;
  r = c1 / c;
S10:
  n = n + 1.0e0;
  t = n / *a;
  w = n * (*b - n) * *x;
  e = *a / s;
  alpha = p * (p + c0) * e * e * (w * *x);
  e = (1.0e0 + t) / (c1 + t + t);
  beta = n + w / s + e * (c + n * yp1);
  p = 1.0e0 + t;
  s += 2.0e0;
  t = alpha * an + beta * anp1;
  an = anp1;
  anp1 = t;
  t = alpha * bn + beta * bnp1;
  bn = bnp1;
  bnp1 = t;
  r0 = r;
  r = anp1 / bnp1;
  if (fabs(r - r0) <= (*eps) * r) {
    goto S20;
  }
  an /= bnp1;
  bn /= bnp1;
  anp1 = r;
  bnp1 = 1.0e0;
  goto S10;
S20:
  bfrac = bfrac * r;
  return bfrac;
}
```