user: What is the c source code of the assembly code below?

```asm
<polint>:
  endbr64
  push   %r15
  mov    %rdi,%r15
  push   %r14
  push   %r13
  xor    %r13d,%r13d
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  push   %rbx
  mov    $0x1,%ebx
  sub    $0x78,%rsp
  mov    %rdx,0x20(%rsp)
  lea    0x30(%rsp),%r14
  lea    0x50(%rsp),%rbp
  mov    %rcx,0x28(%rsp)
  mov    %fs:0x28,%rdx
  mov    %rdx,0x68(%rsp)
  xor    %edx,%edx
  movss  %xmm0,0x8(%rsp)
  subss  (%rdi),%xmm0
  call   <fabs@plt>
  movss  0x8(%rsp),%xmm4
  movaps %xmm0,%xmm1
L69:
  movaps %xmm4,%xmm0
  subss  -0x4(%r15,%rbx,4),%xmm0
  movss  %xmm1,0xc(%rsp)
  movss  %xmm4,0x8(%rsp)
  call   <fabs@plt>
  movss  0xc(%rsp),%xmm1
  movss  0x8(%rsp),%xmm4
  comiss %xmm0,%xmm1
  minss  %xmm1,%xmm0
  movaps %xmm0,%xmm1
  movss  -0x4(%r12,%rbx,4),%xmm0
  cmova  %ebx,%r13d
  movss  %xmm0,-0x4(%r14,%rbx,4)
  movss  %xmm0,-0x4(%rbp,%rbx,4)
  add    $0x1,%rbx
  cmp    $0x6,%rbx
  jne    L69
  lea    -0x1(%r13),%edx
  movslq %r13d,%r13
  mov    0x20(%rsp),%rax
  pxor   %xmm5,%xmm5
  movss  (%r12,%r13,4),%xmm0
  lea    0x4(%r15),%r13
  mov    $0x4,%r12d
  movss  %xmm0,(%rax)
L13:
  mov    $0x1,%ebx
L2:
  movss  -0x4(%r15,%rbx,4),%xmm0
  movss  -0x4(%r13,%rbx,4),%xmm2
  movss  (%r14,%rbx,4),%xmm1
  subss  -0x4(%rbp,%rbx,4),%xmm1
  subss  %xmm4,%xmm0
  subss  %xmm4,%xmm2
  movaps %xmm0,%xmm3
  subss  %xmm2,%xmm3
  ucomiss %xmm5,%xmm3
  jp     L43
  je     L61
L43:
  divss  %xmm3,%xmm1
  mulss  %xmm1,%xmm2
  mulss  %xmm1,%xmm0
  movss  %xmm2,-0x4(%rbp,%rbx,4)
  movss  %xmm0,-0x4(%r14,%rbx,4)
  add    $0x1,%rbx
  cmp    %ebx,%r12d
  jge    L2
  lea    (%rdx,%rdx,1),%eax
  cmp    %r12d,%eax
  jge    L10
  movslq %edx,%rax
  movss  0x30(%rsp,%rax,4),%xmm0
L27:
  mov    0x28(%rsp),%rax
  add    $0x4,%r13
  movss  %xmm0,(%rax)
  mov    0x20(%rsp),%rax
  addss  (%rax),%xmm0
  movss  %xmm0,(%rax)
  sub    $0x1,%r12d
  jne    L13
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L31
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L61:
  lea    D5(%rip),%rdi
  mov    %edx,0x8(%rsp)
  movss  %xmm4,0x1c(%rsp)
  movss  %xmm3,0x18(%rsp)
  movss  %xmm1,0x14(%rsp)
  movss  %xmm2,0x10(%rsp)
  movss  %xmm0,0xc(%rsp)
  call   <printf@plt>
  movss  0x1c(%rsp),%xmm4
  mov    0x8(%rsp),%edx
  pxor   %xmm5,%xmm5
  movss  0x18(%rsp),%xmm3
  movss  0x14(%rsp),%xmm1
  movss  0x10(%rsp),%xmm2
  movss  0xc(%rsp),%xmm0
  jmp    L43
  nopl   0x0(%rax)
L10:
  sub    $0x1,%edx
  movslq %edx,%rax
  movss  0x50(%rsp,%rax,4),%xmm0
  jmp    L27
L31:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: label: D5
data_type: string

tool: D5: problem in polint


assistant: ```c
void polint(float *xa, float *ya, float x, float *y, float *dy) {
  int i, m, ns = 0;
  float den, dif, dift, ho, hp, w;
  float c[5], d[5];
  dif = fabs(x - xa[0]);
  for (i = 1; i <= 5; i++) {
    if ((dift = fabs(x - xa[i - 1])) < dif) {
      ns = i;
      dif = dift;
    }
    c[i - 1] = ya[i - 1];
    d[i - 1] = ya[i - 1];
  }
  *y = ya[ns--];
  for (m = 1; m < 5; m++) {
    for (i = 1; i <= (5 - m); i++) {
      ho = xa[i - 1] - x;
      hp = xa[i + m - 1] - x;
      w = c[i] - d[i - 1];
      if ((den = ho - hp) == 0.0)
        printf("problem in polint\n");
      den = w / den;
      d[i - 1] = hp * den;
      c[i - 1] = ho * den;
    }
    *y += (*dy = (2 * ns < (5 - m) ? c[ns] : d[(ns--) - 1]));
  }
}
```