user: 这段汇编代码用C语言表示是什么样的?

```asm
<betacf>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  movsd  %xmm0,-0x58(%rbp)
  movsd  %xmm1,-0x60(%rbp)
  movsd  %xmm2,-0x68(%rbp)
  movsd  -0x58(%rbp),%xmm0
  addsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x58(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x58(%rbp),%xmm0
  movsd  D246(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  D246(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x20(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L42
  movsd  D97(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
L42:
  movsd  D246(%rip),%xmm0
  divsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movl   $0x1,-0x48(%rbp)
  jmp    L33
L52:
  mov    -0x48(%rbp),%eax
  add    %eax,%eax
  mov    %eax,-0x44(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x48(%rbp),%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2sdl -0x48(%rbp),%xmm2
  movsd  -0x60(%rbp),%xmm0
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x68(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x44(%rbp),%xmm1
  movapd %xmm1,%xmm2
  addsd  -0x18(%rbp),%xmm2
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x44(%rbp),%xmm1
  addsd  -0x58(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L27
  movsd  D97(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
L27:
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x40(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L85
  movsd  D97(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
L85:
  movsd  D246(%rip),%xmm0
  divsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movsd  -0x30(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x48(%rbp),%xmm0
  addsd  -0x58(%rbp),%xmm0
  movq   D250(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x48(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x68(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x44(%rbp),%xmm1
  movapd %xmm1,%xmm2
  addsd  -0x58(%rbp),%xmm2
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x44(%rbp),%xmm1
  addsd  -0x20(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L83
  movsd  D97(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
L83:
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x40(%rbp),%xmm1
  movsd  D246(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L18
  movsd  D97(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
L18:
  movsd  D246(%rip),%xmm0
  divsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x30(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movsd  D246(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movq   D156(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D260(%rip),%xmm0
  comisd %xmm1,%xmm0
  ja     L98
  addl   $0x1,-0x48(%rbp)
L33:
  cmpl   $0x3e8,-0x48(%rbp)
  jle    L52
  jmp    L8
L98:
  nop
L8:
  cmpl   $0x3e8,-0x48(%rbp)
  jle    L3
  lea    D27(%rip),%rax
  mov    %rax,%rsi
  mov    $0x5d,%edi
  call   <error2@plt>
  mov    $0x3,%edi
  call   <exit@plt>
L3:
  movsd  -0x30(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D246
type_guessed: f64
name: extract_symbolic_data
arguments: label: D156
type_guessed: qword
name: extract_symbolic_data
arguments: label: D97
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D250
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D260
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D27
type_guessed: string

tool: D246: 1

tool: D156: 9223372036854775807

tool: {"D97": [0xa0,0xc2,0xeb,0xfe,0x4b,0x48,0xb4,0x39]}

tool: {"D250": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: {"D260": [0x76,0x83,0x0d,0xf4,0xf5,0x21,0x94,0x3e]}

tool: D27: a or b too big, or MAXIT too small in betacf



assistant: ```c
double betacf(double a, double b, double x) {
  int m, m2;
  double aa, c, d, del, h, qab, qam, qap;
  qab = a + b;
  qap = a + 1.0;
  qam = a - 1.0;
  c = 1.0;
  d = 1.0 - qab * x / qap;
  if (fabs(d) < 1.0e-30)
    d = 1.0e-30;
  d = 1.0 / d;
  h = d;
  for (m = 1; m <= 1000; m++) {
    m2 = 2 * m;
    aa = m * (b - m) * x / ((qam + m2) * (a + m2));
    d = 1.0 + aa * d;
    if (fabs(d) < 1.0e-30)
      d = 1.0e-30;
    c = 1.0 + aa / c;
    if (fabs(c) < 1.0e-30)
      c = 1.0e-30;
    d = 1.0 / d;
    h *= d * c;
    aa = -(a + m) * (qab + m) * x / ((a + m2) * (qap + m2));
    d = 1.0 + aa * d;
    if (fabs(d) < 1.0e-30)
      d = 1.0e-30;
    c = 1.0 + aa / c;
    if (fabs(c) < 1.0e-30)
      c = 1.0e-30;
    d = 1.0 / d;
    del = d * c;
    h *= del;
    if (fabs(del - 1.0) < 3.0e-7)
      break;
  }
  if (m > 1000) {
    error2(93, "a or b too big, or MAXIT too small in betacf\n\n");
    exit(3);
  }
  return h;
}
```