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

```asm
<main>:
  endbr64
  push   %r15
  lea    D86(%rip),%rdi
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  lea    D88(%rip),%rbp
  push   %rbx
  sub    $0x108,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xf8(%rsp)
  xor    %eax,%eax
  call   <printf@plt>
  lea    0x14(%rsp),%rsi
  mov    %rbp,%rdi
  xor    %eax,%eax
  call   <scanf@plt>
  mov    0x14(%rsp),%eax
  test   %eax,%eax
  jle    L27
  lea    0x20(%rsp),%rax
  xor    %ebx,%ebx
  lea    0x50(%rsp),%r14
  xor    %r13d,%r13d
  mov    %rax,0x8(%rsp)
  lea    D44(%rip),%r15
  lea    D83(%rip),%r12
  nopl   0x0(%rax)
L42:
  add    $0x1,%r13d
  mov    %r15,%rdi
  xor    %eax,%eax
  mov    %r13d,%esi
  call   <printf@plt>
  mov    0x8(%rsp),%rax
  lea    (%r14,%rbx,1),%rdx
  mov    %r12,%rdi
  lea    (%rax,%rbx,1),%rsi
  xor    %eax,%eax
  add    $0x4,%rbx
  call   <scanf@plt>
  mov    0x14(%rsp),%eax
  cmp    %eax,%r13d
  jl     L42
L27:
  lea    -0x1(%rax),%r9d
  xor    %r8d,%r8d
  lea    0x58(%rsp),%r10
  xor    %ebx,%ebx
  lea    -0x2(%rax),%esi
  lea    0x84(%rsp),%r11
L158:
  cmp    %ebx,%r9d
  jle    L62
  mov    %r9d,%eax
  sub    %ebx,%eax
  test   %eax,%eax
  jle    L64
  mov    %esi,%edi
  lea    (%r10,%r8,1),%rdx
  lea    0x54(%rsp),%rcx
  sub    %ebx,%edi
  lea    0x80(%rsp),%rax
  lea    (%r11,%rdi,4),%rdi
  jmp    L22
  nop
L63:
  movss  0x4(%rdx),%xmm0
  subss  (%rdx),%xmm0
  movss  %xmm0,(%rax,%r8,1)
L18:
  add    $0x4,%rax
  add    $0x4,%rdx
  add    $0x4,%rcx
  cmp    %rax,%rdi
  je     L64
L22:
  test   %ebx,%ebx
  jne    L63
  movss  (%rcx),%xmm0
  subss  -0x4(%rcx),%xmm0
  movss  %xmm0,(%rax)
  jmp    L18
  nopl   0x0(%rax,%rax,1)
L64:
  add    $0x1,%ebx
  add    $0x28,%r8
  cmp    $0x3,%ebx
  jne    L158
L62:
  lea    D109(%rip),%r13
  lea    0x18(%rsp),%r12
  lea    D38(%rip),%r15
  lea    D94(%rip),%r14
  xchg   %ax,%ax
L169:
  mov    %r13,%rdi
  call   <puts@plt>
  xor    %eax,%eax
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <scanf@plt>
  cmpl   $0x3,0x18(%rsp)
  je     L142
  mov    %r15,%rdi
  xor    %eax,%eax
  call   <printf@plt>
  lea    0x1c(%rsp),%rsi
  mov    %r14,%rdi
  xor    %eax,%eax
  call   <scanf@plt>
  mov    0x18(%rsp),%eax
  movss  0x1c(%rsp),%xmm0
  cmp    $0x1,%eax
  je     L197
  cmp    $0x2,%eax
  jne    L204
  movslq 0x14(%rsp),%rax
  movslq %ebx,%rbx
  movaps %xmm0,%xmm3
  movss  D119(%rip),%xmm5
  movss  0x20(%rsp,%rbx,4),%xmm1
  subss  0x20(%rsp),%xmm1
  xor    %ebx,%ebx
  lea    -0x1(%rax),%edx
  movaps %xmm5,%xmm2
  lea    0x80(%rsp,%rax,4),%rax
  movslq %edx,%rcx
  subss  0x20(%rsp,%rcx,4),%xmm3
  movss  0x50(%rsp,%rcx,4),%xmm6
  movss  %xmm6,0x4(%rsp)
  divss  %xmm1,%xmm3
L31:
  cmp    %ebx,%edx
  jle    L204
  add    $0x1,%ebx
  pxor   %xmm4,%xmm4
  movaps %xmm3,%xmm1
  add    $0x24,%rax
  cvtsi2ss %ebx,%xmm4
  addss  %xmm5,%xmm1
  divss  %xmm4,%xmm1
  mulss  %xmm1,%xmm2
  movss  -0x2c(%rax),%xmm1
  mulss  %xmm2,%xmm1
  addss  0x4(%rsp),%xmm1
  movss  %xmm1,0x4(%rsp)
  cmp    $0x3,%ebx
  jne    L31
L204:
  pxor   %xmm1,%xmm1
  lea    D20(%rip),%rdi
  mov    $0x2,%eax
  cvtss2sd %xmm0,%xmm0
  cvtss2sd 0x4(%rsp),%xmm1
  call   <printf@plt>
  jmp    L169
  nopl   (%rax)
L197:
  movslq %ebx,%rbx
  movss  0x20(%rsp),%xmm3
  movaps %xmm0,%xmm1
  mov    0x14(%rsp),%eax
  movss  0x20(%rsp,%rbx,4),%xmm2
  movss  0x50(%rsp),%xmm4
  subss  %xmm3,%xmm1
  lea    -0x1(%rax),%ebx
  subss  %xmm3,%xmm2
  movss  %xmm4,0x4(%rsp)
  divss  %xmm2,%xmm1
  test   %ebx,%ebx
  jle    L15
  movss  0x80(%rsp),%xmm2
  movaps %xmm1,%xmm3
  xorps  D34(%rip),%xmm3
  mulss  %xmm3,%xmm2
  movaps %xmm2,%xmm7
  addss  %xmm4,%xmm7
  movss  %xmm7,0x4(%rsp)
  cmp    $0x1,%ebx
  je     L204
  movaps %xmm1,%xmm2
  subss  D119(%rip),%xmm2
  divss  D59(%rip),%xmm2
  mulss  %xmm3,%xmm2
  movss  0xa8(%rsp),%xmm3
  mulss  %xmm2,%xmm3
  movaps %xmm3,%xmm5
  addss  %xmm7,%xmm5
  movss  %xmm5,0x4(%rsp)
  cmp    $0x2,%ebx
  je     L204
  subss  D24(%rip),%xmm1
  mov    $0x3,%ebx
  mulss  %xmm2,%xmm1
  mulss  0xd0(%rsp),%xmm1
  addss  %xmm5,%xmm1
  movss  %xmm1,0x4(%rsp)
  jmp    L204
  nopl   0x0(%rax)
L142:
  mov    0xf8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L184
  add    $0x108,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L15:
  xor    %ebx,%ebx
  jmp    L204
L184:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D86
data_type: string
name: unpack_data
arguments: data_label: D88
data_type: string
name: unpack_data
arguments: data_label: D44
data_type: string
name: unpack_data
arguments: data_label: D83
data_type: string
name: unpack_data
arguments: data_label: D109
data_type: string
name: unpack_data
arguments: data_label: D38
data_type: string
name: unpack_data
arguments: data_label: D94
data_type: string
name: unpack_data
arguments: data_label: D119
data_type: byte[4]
name: unpack_data
arguments: data_label: D20
data_type: string
name: unpack_data
arguments: data_label: D34
data_type: byte[16]
name: unpack_data
arguments: data_label: D59
data_type: byte[4]
name: unpack_data
arguments: data_label: D24
data_type: byte[4]

tool: D86: 
No. of pair of points: 

tool: D88: %d

tool: D44: %dth value of x and y wanted: 

tool: D83: %f %f

tool: D109: 
 Give option: Press 1 for Forward
Press 2 for Backward
3 for exit

tool: D38: 
Give value of x for which y to be interpolated: 

tool: D94: %f

tool: {"D119": [0x00,0x00,0x80,0x3f]}

tool: D20: x = %f	y = %f


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

tool: {"D59": [0x00,0x00,0x00,0x00]}

tool: {"D24": [0x00,0x00,0x00,0x40]}

assistant: ```c
int main() {
  int i, n, j, op;
  float x[10], y[10], xx, yy, d[3][10], u, t;
  printf("\nNo. of pair of points: ");
  scanf("%d", &n);
  for (i = 0; i < n; i++) {
    printf("%dth value of x and y wanted: ", i + 1);
    scanf("%f %f", &x[i], &y[i]);
  }
  for (i = 0; i < 3 && i < n - 1; i++) {
    for (j = 0; j < n - i - 1; j++) {
      if (i == 0)
        d[i][j] = y[j + 1] - y[j];
      else
        d[i][j] = d[i - 1][j + 1] - d[i - 1][j];
    }
  }
  while (1) {
    printf("\n Give option: Press 1 for Forward\nPress 2 for Backward\n3 for exit\n");
    scanf("%d", &op);
    if (op == 3)
      break;
    printf("\nGive value of x for which y to be interpolated: ");
    scanf("%f", &xx);
    switch (op) {
    case 1:
      u = (xx - x[0]) / (x[i] - x[0]);
      yy = y[0];
      t = 1;
      for (i = 0; i < 3 && i < n - 1; i++) {
        t *= (u - i) / (i - 1);
        yy = yy + t * d[i][0];
      }
      break;
    case 2:
      u = (xx - x[n - 1]) / (x[i] - x[0]);
      yy = y[n - 1];
      t = 1;
      for (i = 0; i < 3 && i < n - 1; i++) {
        t *= (u + 1) / (i + 1);
        yy = yy + t * d[i][n - i - 2];
      }
    }
    printf("x = %f\ty = %f\n", xx, yy);
  }
  return 0;
}
```