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

```asm
<plgndr_a>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %edi,-0x54(%rbp)
  mov    %esi,-0x58(%rbp)
  movsd  %xmm0,-0x60(%rbp)
  movsd  D5(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  D9(%rip),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x60(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  cmpl   $0x0,-0x58(%rbp)
  jle    L5
  movsd  -0x18(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x10(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x8(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movl   $0x1,-0x50(%rbp)
  jmp    L80
L45:
  movsd  -0x38(%rbp),%xmm0
  movq   D25(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movsd  -0x8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x48(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x50(%rbp)
L80:
  mov    -0x50(%rbp),%eax
  cmp    -0x58(%rbp),%eax
  jle    L45
L5:
  mov    -0x54(%rbp),%eax
  cmp    -0x58(%rbp),%eax
  jne    L82
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  jmp    L7
L82:
  mov    -0x58(%rbp),%eax
  add    %eax,%eax
  add    $0x1,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x58(%rbp),%eax
  add    $0x1,%eax
  cmp    %eax,-0x54(%rbp)
  jne    L25
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  jmp    L7
L25:
  mov    -0x58(%rbp),%eax
  add    $0x2,%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L23
L73:
  mov    -0x4c(%rbp),%eax
  add    %eax,%eax
  sub    $0x1,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  mov    -0x4c(%rbp),%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  sub    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x4c(%rbp),%eax
  sub    -0x58(%rbp),%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  addl   $0x1,-0x4c(%rbp)
L23:
  mov    -0x4c(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jle    L73
  movsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
L7:
  movsd  -0x28(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: label: D5
datatype: f64
name: retrieve_data
arguments: label: D9
datatype: f64
name: retrieve_data
arguments: label: D25
datatype: byte[8]

tool: D5: 2

tool: D9: 1

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

assistant: ```c
double plgndr_a(l, m, t)
int l, m;
double t;
{
  int i, ll;
  double x, fact, pll, pmm, pmmp1, somx2, plgndr;
  const double two = 2.0;
  const double one = 1.0;
  x = cos(t);
  pmm = one;
  if (m > 0) {
    somx2 = sqrt((one - x) * (one + x));
    fact = one;
    for (i = 1; i <= m; i++) {
      pmm = -pmm * fact * somx2;
      fact = fact + two;
    }
  }
  if (l == m)
    plgndr = pmm;
  else {
    pmmp1 = x * (2 * m + 1) * pmm;
    if (l == m + 1)
      plgndr = pmmp1;
    else {
      for (ll = m + 2; ll <= l; ll++) {
        pll = (x * (2 * ll - 1) * pmmp1 - (ll + m - 1) * pmm) / (ll - m);
        pmm = pmmp1;
        pmmp1 = pll;
      }
      plgndr = pll;
    }
  }
  return plgndr;
}
```