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

```asm
<main>:
  endbr64
  push   %r15
  mov    $0x4,%ecx
  mov    $0x1,%edx
  xor    %eax,%eax
  push   %r14
  mov    $0x4,%esi
  mov    $0x1,%edi
  push   %r13
  push   %r12
  push   %rbp
  lea    D128(%rip),%rbp
  push   %rbx
  sub    $0x28,%rsp
  call   <matrix@plt>
  mov    $0x4,%esi
  mov    $0x1,%edi
  cltq
  mov    %rax,0x8(%rsp)
  xor    %eax,%eax
  call   <vector@plt>
  mov    $0x4,%esi
  mov    $0x1,%edi
  movslq %eax,%r14
  xor    %eax,%eax
  call   <vector@plt>
  movl   $0xffffffff,0x14(%rsp)
  movslq %eax,%r13
  mov    %r13,%rbx
  sub    %r14,%rbx
L46:
  pxor   %xmm5,%xmm5
  movl   $0x1,0x10(%rsp)
  cvtsi2sdl 0x14(%rsp),%xmm5
  movsd  %xmm5,0x18(%rsp)
L38:
  mov    0x10(%rsp),%r15d
  lea    D103(%rip),%rdi
  xor    %eax,%eax
  lea    0x4(%r13),%r12
  mov    %r15d,%esi
  call   <printf@plt>
  pxor   %xmm0,%xmm0
  pxor   %xmm4,%xmm4
  cvtsi2sd %r15d,%xmm0
  mulsd  D97(%rip),%xmm0
  mov    $0x1,%r15d
  mulsd  0x18(%rsp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  cvtss2sd %xmm0,%xmm4
  movsd  %xmm4,(%rsp)
L45:
  pxor   %xmm0,%xmm0
  mov    %r15d,%edx
  mov    $0x1,%eax
  add    $0x4,%r12
  cvtsi2sd %r15d,%xmm0
  mulsd  D104(%rip),%xmm0
  lea    D56(%rip),%rcx
  add    $0x1,%r15d
  addsd  (%rsp),%xmm0
  lea    D55(%rip),%rsi
  lea    D112(%rip),%rdi
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,-0x4(%r12)
  cvtss2sd %xmm0,%xmm0
  call   <printf@plt>
  cmp    $0x5,%r15d
  jne    L45
  mov    $0xa,%edi
  call   <putchar@plt>
  movl   $0x5,(%rsp)
L43:
  mov    D29(%rip),%rax
  mov    $0x4,%edx
  mov    %r13,%rsi
  lea    0x4(%r14),%r12
  mov    $0x1,%edi
  mov    $0x1,%r15d
  movq   %rax,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <mnewt@plt>
  mov    0x8(%rsp),%rcx
  mov    %r14,%rdx
  mov    %r13,%rdi
  mov    $0x4,%esi
  xor    %eax,%eax
  call   <usrfun@plt>
  lea    D12(%rip),%rcx
  xor    %eax,%eax
  lea    D107(%rip),%rdx
  lea    D60(%rip),%rsi
  lea    D35(%rip),%rdi
  call   <printf@plt>
L1:
  mov    %r15d,%esi
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  mov    %rbp,%rdi
  mov    $0x2,%eax
  add    $0x1,%r15d
  cvtss2sd (%rbx,%r12,1),%xmm0
  add    $0x4,%r12
  cvtss2sd -0x4(%r12),%xmm1
  call   <printf@plt>
  cmp    $0x5,%r15d
  jne    L1
  lea    D15(%rip),%rdi
  call   <puts@plt>
  xor    %eax,%eax
  call   <getchar@plt>
  subl   $0x1,(%rsp)
  jne    L43
  addl   $0x1,0x10(%rsp)
  mov    0x10(%rsp),%eax
  cmp    $0x4,%eax
  jne    L38
  addl   $0x2,0x14(%rsp)
  mov    0x14(%rsp),%eax
  cmp    $0x3,%eax
  jne    L46
  mov    %r13,%rdi
  mov    $0x4,%edx
  mov    $0x1,%esi
  xor    %eax,%eax
  call   <free_vector@plt>
  mov    %r14,%rdi
  mov    $0x4,%edx
  xor    %eax,%eax
  mov    $0x1,%esi
  call   <free_vector@plt>
  mov    0x8(%rsp),%rdi
  mov    $0x4,%r8d
  mov    $0x1,%ecx
  mov    $0x4,%edx
  mov    $0x1,%esi
  xor    %eax,%eax
  call   <free_matrix@plt>
  add    $0x28,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: label: D128
type_infered: string
name: parse_data
arguments: label: D103
type_infered: string
name: parse_data
arguments: label: D97
type_infered: f64
name: parse_data
arguments: label: D104
type_infered: f64
name: parse_data
arguments: label: D56
type_infered: string
name: parse_data
arguments: label: D55
type_infered: string
name: parse_data
arguments: label: D112
type_infered: string
name: parse_data
arguments: label: D29
type_infered: byte[8]
name: parse_data
arguments: label: D12
type_infered: string
name: parse_data
arguments: label: D107
type_infered: string
name: parse_data
arguments: label: D60
type_infered: string
name: parse_data
arguments: label: D35
type_infered: string
name: parse_data
arguments: label: D15
type_infered: string

tool: D128: %5d %14.6f %15.6f


tool: D103: Starting vector number %2d


tool: D97: 0.2001

tool: D104: 0.2

tool: D56: ] = 

tool: D55: x[

tool: D112: %7s%1d%s %5.2f


tool: {"D29": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xb0,0x3e]}

tool: D12: f

tool: D107: x[i]

tool: D60: i

tool: D35: %5s %13s %13s


tool: D15: 
press RETURN to continue...

assistant: ```c
int main(void) {
  int i, j, k, kk;
  float xx, *x, *fvec, **fjac;
  fjac = matrix(1, 4, 1, 4);
  fvec = vector(1, 4);
  x = vector(1, 4);
  for (kk = 1; kk <= 2; kk++) {
    for (k = 1; k <= 3; k++) {
      xx = 0.2001 * k * (2 * kk - 3);
      printf("Starting vector number %2d\n", k);
      for (i = 1; i <= 4; i++) {
        x[i] = xx + 0.2 * i;
        printf("%7s%1d%s %5.2f\n", "x[", i, "] = ", x[i]);
      }
      printf("\n");
      for (j = 1; j <= 5; j++) {
        mnewt(1, x, 4, 1.0e-6, 1.0e-6);
        usrfun(x, 4, fvec, fjac);
        printf("%5s %13s %13s\n", "i", "x[i]", "f");
        for (i = 1; i <= 4; i++)
          printf("%5d %14.6f %15.6f\n", i, x[i], fvec[i]);
        printf("\npress RETURN to continue...\n");
        (void)getchar();
      }
    }
  }
  free_vector(x, 1, 4);
  free_vector(fvec, 1, 4);
  free_matrix(fjac, 1, 4, 1, 4);
  return 0;
}
```