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

```asm
<calculoCoeficienteLinear>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,-0x5c(%rbp)
  mov    %rsp,%rcx
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  movslq %edi,%rax
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rsi
  mov    %rax,%rdx
  and    $0xfffffffffffff000,%rsi
  and    $0xfffffffffffffff0,%rdx
  sub    %rsi,%rcx
  cmp    %rcx,%rsp
  je     L68
L103:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L103
L68:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  jne    L39
L87:
  mov    %rax,%rdx
  mov    %rsp,%rcx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,-0x58(%rbp)
  sub    %rax,%rcx
  and    $0xfffffffffffffff0,%rdx
  cmp    %rcx,%rsp
  je     L86
L23:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L23
L86:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  jne    L66
L28:
  mov    %rsp,%rbx
  lea    D2(%rip),%rsi
  movsd  %xmm2,-0x50(%rbp)
  lea    D24(%rip),%rdi
  movsd  %xmm0,-0x48(%rbp)
  mov    %rbx,-0x68(%rbp)
  call   <fopen@plt>
  movsd  -0x48(%rbp),%xmm0
  movsd  -0x50(%rbp),%xmm2
  mov    %rax,%r15
  mov    -0x5c(%rbp),%eax
  cmp    $0x1,%eax
  jle    L64
  mov    -0x58(%rbp),%rcx
  sub    $0x2,%eax
  pxor   %xmm1,%xmm1
  lea    0x8(%rbx),%r14
  movapd %xmm1,%xmm0
  lea    D29(%rip),%rbx
  lea    0x8(%rcx),%r12
  lea    0x10(%rcx,%rax,8),%r13
  nopl   0x0(%rax,%rax,1)
L94:
  mov    %r12,%rdx
  mov    %rbx,%rsi
  mov    %r15,%rdi
  movsd  %xmm1,-0x50(%rbp)
  movsd  %xmm0,-0x48(%rbp)
  add    $0x8,%r12
  call   <fscanf@plt>
  mov    %r14,%rdx
  mov    %rbx,%rsi
  mov    %r15,%rdi
  call   <fscanf@plt>
  movsd  -0x48(%rbp),%xmm0
  movsd  -0x50(%rbp),%xmm1
  add    $0x8,%r14
  addsd  -0x8(%r12),%xmm0
  addsd  -0x8(%r14),%xmm1
  cmp    %r13,%r12
  jne    L94
  mov    -0x5c(%rbp),%r12d
  pxor   %xmm5,%xmm5
  mov    -0x58(%rbp),%rax
  cvtsi2sd %r12d,%xmm5
  sub    $0x1,%r12d
  movslq %r12d,%r12
  movsd  (%rax,%r12,8),%xmm2
  mov    -0x68(%rbp),%rax
  divsd  %xmm5,%xmm0
  movsd  (%rax,%r12,8),%xmm4
  movapd %xmm4,%xmm3
  divsd  %xmm5,%xmm1
  subsd  %xmm0,%xmm3
  mulsd  %xmm2,%xmm3
  movapd %xmm3,%xmm0
  subsd  %xmm1,%xmm2
  mulsd  %xmm4,%xmm2
L64:
  divsd  %xmm2,%xmm0
  mov    %r15,%rdi
  movsd  %xmm0,-0x48(%rbp)
  call   <fclose@plt>
  movsd  -0x48(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L19
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
  cs nopw 0x0(%rax,%rax,1)
L39:
  orq    $0x0,-0x8(%rsp,%rdx,1)
  jmp    L87
L66:
  orq    $0x0,-0x8(%rsp,%rdx,1)
  jmp    L28
L19:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: label: D2
datatype: string
name: unpack_data
arguments: label: D24
datatype: string
name: unpack_data
arguments: label: D29
datatype: string

tool: D2: rt

tool: D24: dadosMinimosQuadrados.txt

tool: D29: %lf

assistant: ```c
double calculoCoeficienteLinear(int quantidadeVelas) {
  FILE *arquivo;
  double x[quantidadeVelas], y[quantidadeVelas];
  double soma_x = 0, soma_y = 0;
  double numerador, denominador;
  double variacaoLinear;
  int i;
  arquivo = fopen("dadosMinimosQuadrados.txt", "rt");
  for (i = 1; i < quantidadeVelas; i++) {
    fscanf(arquivo, "%lf", &x[i]);
    fscanf(arquivo, "%lf", &y[i]);
    soma_x = soma_x + x[i];
    soma_y = soma_y + y[i];
  }
  for (i = 1; i < quantidadeVelas; i++) {
    numerador = x[i] * (y[i] - soma_x / quantidadeVelas);
    denominador = y[i] * (x[i] - soma_y / quantidadeVelas);
  }
  variacaoLinear = numerador / denominador;
  fclose(arquivo);
  return variacaoLinear;
}
```