user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<main>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x90,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x88(%rsp)
  xor    %eax,%eax
  lea    D136(%rip),%rsi
  lea    D65(%rip),%rdi
  call   <fopen@plt>
  mov    %rax,%rbx
  test   %rax,%rax
  je     L58
L64:
  lea    0x30(%rsp),%rdi
  mov    %rbx,%rdx
  mov    $0x50,%esi
  call   <fgets@plt>
  lea    D70(%rip),%r12
  lea    0x30(%rsp),%rbp
  lea    D48(%rip),%r13
L37:
  mov    $0x1a,%edx
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <strncmp@plt>
  test   %rax,%rax
  je     L34
  mov    %rbx,%rdx
  mov    $0x50,%esi
  mov    %rbp,%rdi
  call   <fgets@plt>
  mov    %rbx,%rdi
  call   <feof@plt>
  test   %rax,%rax
  je     L37
  mov    %r13,%rdi
  call   <nrerror@plt>
  jmp    L37
L58:
  lea    D66(%rip),%rdi
  call   <nrerror@plt>
  jmp    L64
L34:
  lea    0x14(%rsp),%rdx
  lea    D44(%rip),%rsi
  mov    %rbx,%rdi
  call   <fscanf@plt>
  lea    0x30(%rsp),%rsi
  lea    D111(%rip),%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  sub    $0x8,%rsp
  lea    D46(%rip),%rax
  push   %rax
  lea    D2(%rip),%r9
  lea    D162(%rip),%r8
  lea    D67(%rip),%rcx
  lea    D139(%rip),%rdx
  lea    D61(%rip),%rsi
  lea    D106(%rip),%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  add    $0x10,%rsp
  cmpl   $0x0,0x14(%rsp)
  jle    L66
  mov    $0x1,%ebp
  lea    0x20(%rsp),%r14
  lea    D118(%rip),%r13
  lea    D88(%rip),%r12
L45:
  lea    0x1c(%rsp),%rcx
  lea    0x18(%rsp),%rdx
  mov    %r14,%r8
  mov    %r13,%rsi
  mov    %rbx,%rdi
  mov    $0x0,%eax
  call   <fscanf@plt>
  movss  D141(%rip),%xmm6
  subss  0x18(%rsp),%xmm6
  lea    0x2c(%rsp),%rdx
  lea    0x28(%rsp),%rsi
  lea    0x24(%rsp),%rdi
  movss  %xmm6,0xc(%rsp)
  movaps %xmm6,%xmm1
  movss  0x1c(%rsp),%xmm0
  call   <sncndn@plt>
  movss  0x24(%rsp),%xmm3
  movss  0x2c(%rsp),%xmm1
  movss  0x28(%rsp),%xmm0
  movaps %xmm3,%xmm5
  mulss  0x18(%rsp),%xmm5
  mulss  %xmm3,%xmm5
  mulss  %xmm1,%xmm1
  addss  %xmm1,%xmm5
  movaps %xmm3,%xmm4
  mulss  %xmm3,%xmm4
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm4
  pxor   %xmm0,%xmm0
  cvtss2sd 0xc(%rsp),%xmm0
  cvtss2sd %xmm5,%xmm5
  cvtss2sd %xmm4,%xmm4
  cvtss2sd %xmm3,%xmm3
  pxor   %xmm2,%xmm2
  cvtss2sd 0x20(%rsp),%xmm2
  pxor   %xmm1,%xmm1
  cvtss2sd 0x1c(%rsp),%xmm1
  mov    %r12,%rdi
  mov    $0x6,%eax
  call   <printf@plt>
  add    $0x1,%ebp
  cmp    %ebp,0x14(%rsp)
  jge    L45
L66:
  mov    %rbx,%rdi
  call   <fclose@plt>
  mov    0x88(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L72
  mov    $0x0,%eax
  add    $0x90,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L72:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: tag: D136
data_type: string
name: extract_symbolic_data
arguments: tag: D65
data_type: string
name: extract_symbolic_data
arguments: tag: D70
data_type: string
name: extract_symbolic_data
arguments: tag: D48
data_type: string
name: extract_symbolic_data
arguments: tag: D66
data_type: string
name: extract_symbolic_data
arguments: tag: D44
data_type: string
name: extract_symbolic_data
arguments: tag: D111
data_type: string
name: extract_symbolic_data
arguments: tag: D46
data_type: string
name: extract_symbolic_data
arguments: tag: D2
data_type: string
name: extract_symbolic_data
arguments: tag: D162
data_type: string
name: extract_symbolic_data
arguments: tag: D67
data_type: string
name: extract_symbolic_data
arguments: tag: D139
data_type: string
name: extract_symbolic_data
arguments: tag: D61
data_type: string
name: extract_symbolic_data
arguments: tag: D106
data_type: string
name: extract_symbolic_data
arguments: tag: D118
data_type: string
name: extract_symbolic_data
arguments: tag: D88
data_type: string
name: extract_symbolic_data
arguments: tag: D141
data_type: f32

tool: D136: r

tool: D65: fncval.dat

tool: D70: Jacobian Elliptic Function

tool: D48: Data not found in fncval.dat


tool: D66: Data file fncval.dat not found


tool: D44: %d %*s

tool: D111: 
%s


tool: D46: (mc)*(sn^2)+dn^2

tool: D2: sn^2+cn^2

tool: D162: sn

tool: D67: actual

tool: D139: u

tool: D61: mc

tool: D106: %4s %8s %16s %13s %15s %18s


tool: D118: %f %f %f

tool: D88: %5.2f %8.2f %15.5f %15.5f %12.5f %14.5f


tool: D141: 1

assistant: ```c
int main(void) {
  char txt[80];
  int i, nval;
  float em, emmc, uu, val, sn, cn, dn;
  FILE *fp;
  if ((fp = fopen("fncval.dat", "r")) == NULL)
    nrerror("Data file fncval.dat not found\n");
  fgets(txt, 80, fp);
  while (strncmp(txt, "Jacobian Elliptic Function", 26)) {
    fgets(txt, 80, fp);
    if (feof(fp))
      nrerror("Data not found in fncval.dat\n");
  }
  fscanf(fp, "%d %*s", &nval);
  printf("\n%s\n", txt);
  printf("%4s %8s %16s %13s %15s %18s\n", "mc", "u", "actual", "sn", "sn^2+cn^2", "(mc)*(sn^2)+dn^2");
  for (i = 1; i <= nval; i++) {
    fscanf(fp, "%f %f %f", &em, &uu, &val);
    emmc = 1.0 - em;
    sncndn(uu, emmc, &sn, &cn, &dn);
    printf("%5.2f %8.2f %15.5f %15.5f %12.5f %14.5f\n", emmc, uu, val, sn, (sn * sn + cn * cn), (em * sn * sn + dn * dn));
  }
  fclose(fp);
  return 0;
}
```