user: Please show me the C implementation that matches the assembly code listed below:

```asm
<bcuint>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,0x18(%rsp)
  mov    %rsi,0x20(%rsp)
  mov    %rdx,%r14
  mov    %rcx,%r15
  movss  %xmm0,(%rsp)
  movss  %xmm1,0x10(%rsp)
  movss  %xmm2,0x4(%rsp)
  movss  %xmm3,0x14(%rsp)
  movss  %xmm4,0x28(%rsp)
  movss  %xmm5,0x2c(%rsp)
  mov    %r8,%r13
  mov    %r9,%rbp
  mov    0x70(%rsp),%r12
  mov    $0x4,%ecx
  mov    $0x1,%edx
  mov    $0x4,%esi
  mov    $0x1,%edi
  call   <matrix@plt>
  mov    %rax,%rbx
  movss  0x10(%rsp),%xmm7
  subss  (%rsp),%xmm7
  movaps %xmm7,%xmm0
  movss  %xmm7,0x8(%rsp)
  movss  0x14(%rsp),%xmm3
  subss  0x4(%rsp),%xmm3
  movaps %xmm3,%xmm1
  movss  %xmm3,0xc(%rsp)
  mov    %rax,%r8
  mov    %r15,%rcx
  mov    %r14,%rdx
  mov    0x20(%rsp),%rsi
  mov    0x18(%rsp),%rdi
  call   <bcucof@plt>
  movss  0x10(%rsp),%xmm7
  ucomiss (%rsp),%xmm7
  setnp  %al
  mov    $0x0,%edx
  cmovne %edx,%eax
  test   %al,%al
  jne    L1
  movss  0x14(%rsp),%xmm3
  ucomiss 0x4(%rsp),%xmm3
  setnp  %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L1
L34:
  movss  0x28(%rsp),%xmm2
  subss  (%rsp),%xmm2
  divss  0x8(%rsp),%xmm2
  movss  0x2c(%rsp),%xmm1
  subss  0x4(%rsp),%xmm1
  divss  0xc(%rsp),%xmm1
  movl   $0x0,0x0(%rbp)
  movl   $0x0,(%r12)
  movl   $0x0,0x0(%r13)
  mov    $0x10,%eax
  pxor   %xmm5,%xmm5
  cvtss2sd %xmm1,%xmm5
  movsd  D11(%rip),%xmm4
  pxor   %xmm3,%xmm3
  cvtss2sd %xmm2,%xmm3
L0:
  mov    (%rbx,%rax,2),%rdx
  movaps %xmm1,%xmm0
  mulss  0x10(%rdx),%xmm0
  addss  0xc(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  addss  0x8(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movaps %xmm2,%xmm6
  mulss  0x0(%r13),%xmm6
  addss  %xmm6,%xmm0
  addss  0x4(%rdx),%xmm0
  movss  %xmm0,0x0(%r13)
  mov    (%rbx,%rax,2),%rdx
  pxor   %xmm0,%xmm0
  cvtss2sd 0x10(%rdx),%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm0
  pxor   %xmm6,%xmm6
  cvtss2sd 0xc(%rdx),%xmm6
  addsd  %xmm6,%xmm6
  addsd  %xmm6,%xmm0
  mulsd  %xmm5,%xmm0
  movaps %xmm2,%xmm6
  mulss  (%r12),%xmm6
  cvtss2sd %xmm6,%xmm6
  addsd  %xmm6,%xmm0
  pxor   %xmm6,%xmm6
  cvtss2sd 0x8(%rdx),%xmm6
  addsd  %xmm6,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r12)
  mov    0x20(%rbx),%rdx
  pxor   %xmm0,%xmm0
  cvtss2sd (%rdx,%rax,1),%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  %xmm3,%xmm0
  mov    0x18(%rbx),%rdx
  pxor   %xmm6,%xmm6
  cvtss2sd (%rdx,%rax,1),%xmm6
  addsd  %xmm6,%xmm6
  addsd  %xmm6,%xmm0
  mulsd  %xmm3,%xmm0
  movaps %xmm1,%xmm6
  mulss  0x0(%rbp),%xmm6
  cvtss2sd %xmm6,%xmm6
  addsd  %xmm6,%xmm0
  mov    0x10(%rbx),%rdx
  pxor   %xmm6,%xmm6
  cvtss2sd (%rdx,%rax,1),%xmm6
  addsd  %xmm6,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp)
  sub    $0x4,%rax
  jne    L0
  divss  0x8(%rsp),%xmm0
  movss  %xmm0,0x0(%rbp)
  movss  (%r12),%xmm0
  divss  0xc(%rsp),%xmm0
  movss  %xmm0,(%r12)
  mov    $0x4,%r8d
  mov    $0x1,%ecx
  mov    $0x4,%edx
  mov    $0x1,%esi
  mov    %rbx,%rdi
  call   <free_matrix@plt>
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L1:
  lea    D2(%rip),%rdi
  call   <nrerror@plt>
  jmp    L34
```

tool_calls: name: parse_data
arguments: data_label: D11
type: f64
name: parse_data
arguments: data_label: D2
type: string

tool: D11: 3

tool: D2: Bad input in routine bcuint

assistant: ```c
void bcuint(float y[], float y1[], float y2[], float y12[], float x1l, float x1u, float x2l, float x2u, float x1, float x2, float *ansy, float *ansy1, float *ansy2) {
  void bcucof(float y[], float y1[], float y2[], float y12[], float d1, float d2, float **c);
  int i;
  float t, u, d1, d2, **c;
  c = matrix(1, 4, 1, 4);
  d1 = x1u - x1l;
  d2 = x2u - x2l;
  bcucof(y, y1, y2, y12, d1, d2, c);
  if (x1u == x1l || x2u == x2l)
    nrerror("Bad input in routine bcuint");
  t = (x1 - x1l) / d1;
  u = (x2 - x2l) / d2;
  *ansy = (*ansy2) = (*ansy1) = 0.0;
  for (i = 4; i >= 1; i--) {
    *ansy = t * (*ansy) + ((c[i][4] * u + c[i][3]) * u + c[i][2]) * u + c[i][1];
    *ansy2 = t * (*ansy2) + (3.0 * c[i][4] * u + 2.0 * c[i][3]) * u + c[i][2];
    *ansy1 = u * (*ansy1) + (3.0 * c[4][i] * t + 2.0 * c[3][i]) * t + c[2][i];
  }
  *ansy1 /= d1;
  *ansy2 /= d2;
  free_matrix(c, 1, 4, 1, 4);
}
```