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

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  mov    %rsi,%rbx
  sub    $0xd8,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  cmp    $0x5,%edi
  jle    L17
L105:
  mov    0x28(%rbx),%rdi
  lea    D43(%rip),%r12
  call   <atof@plt>
  mov    0x30(%rbx),%rdi
  movsd  %xmm0,-0xd8(%rbp)
  call   <atof@plt>
  mov    0x20(%rbx),%rdi
  movsd  %xmm0,-0xe0(%rbp)
  call   <atoi@plt>
  mov    0x8(%rbx),%rdi
  mov    $0x2,%edx
  mov    %r12,%rsi
  call   <openfile_or_exit@plt>
  lea    -0x70(%rbp),%rcx
  lea    -0x68(%rbp),%r8
  mov    %rax,%rdi
  lea    -0x88(%rbp),%rdx
  lea    -0x84(%rbp),%rsi
  call   <read_slap@plt>
  mov    0x10(%rbx),%rdi
  mov    $0x2,%edx
  mov    %r12,%rsi
  call   <openfile_or_exit@plt>
  lea    -0x60(%rbp),%rcx
  lea    -0x58(%rbp),%r8
  mov    %rax,%rdi
  lea    -0x80(%rbp),%rdx
  lea    -0x7c(%rbp),%rsi
  call   <read_slap@plt>
  mov    0x18(%rbx),%rdi
  mov    $0x2,%edx
  mov    %r12,%rsi
  call   <openfile_or_exit@plt>
  lea    -0x50(%rbp),%rcx
  lea    -0x78(%rbp),%rdx
  mov    %rax,%rdi
  lea    -0x74(%rbp),%rsi
  lea    -0x40(%rbp),%r9
  lea    -0x48(%rbp),%r8
  call   <read_slap_w@plt>
  mov    -0x78(%rbp),%edi
  mov    %rsp,%rcx
  lea    0x1(%rdi),%eax
  lea    0xf(,%rax,4),%rdx
  mov    %rdx,%rax
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rsi
  and    $0xfffffffffffff000,%rsi
  sub    %rsi,%rcx
  cmp    %rcx,%rsp
  je     L542
L5:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L5
L542:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L291
L512:
  shr    $0x4,%rdx
  mov    %rsp,%rax
  mov    %rsp,%r15
  shl    $0x4,%rdx
  mov    %rdx,%rcx
  and    $0xfffffffffffff000,%rcx
  sub    %rcx,%rax
L391:
  cmp    %rax,%rsp
  je     L31
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L391
L17:
  mov    (%rsi),%rsi
  lea    D51(%rip),%rdi
  call   <printf@plt>
  mov    $0x1,%edi
  call   <exit@plt>
  jmp    L105
L291:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L512
L31:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  jne    L88
L234:
  mov    -0x88(%rbp),%edx
  mov    -0x68(%rbp),%r10
  mov    %rsp,%r14
  mov    -0x58(%rbp),%r9
  mov    %rdx,%rax
  lea    (%r10,%rdx,4),%rbx
  mov    -0x80(%rbp),%edx
  cmp    %eax,%edi
  cmovbe %edi,%eax
  mov    %rdx,%rsi
  lea    (%r9,%rdx,4),%r11
  mov    %eax,%r8d
  mov    $0x1,%eax
  add    $0x2,%r8
  nopw   0x0(%rax,%rax,1)
L327:
  mov    -0x4(%r10,%rax,4),%edx
  mov    %eax,%ecx
  mov    %edx,-0x4(%r15,%rax,4)
  lea    -0x1(%rax),%edx
  cmp    %edx,%esi
  jb     L370
  mov    -0x4(%r9,%rax,4),%edx
  mov    %edx,-0x4(%r14,%rax,4)
  mov    %eax,%edx
  add    $0x1,%rax
  cmp    %r8,%rax
  jne    L327
L382:
  cmp    %ecx,%edi
  jb     L551
  mov    (%rbx),%r8d
  movslq %edx,%rax
  jmp    L548
  nopl   0x0(%rax)
L230:
  mov    (%r9,%rax,4),%edx
  lea    0x1(%rax),%ecx
  mov    %edx,(%r14,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%edi
  jb     L81
L548:
  mov    %r8d,(%r15,%rax,4)
  cmp    %ecx,%esi
  jae    L230
  mov    (%r11),%edx
  lea    0x1(%rax),%ecx
  mov    %edx,(%r14,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%edi
  jae    L548
L81:
  lea    0xf(,%rdi,8),%rdx
  mov    %rsp,%rcx
  mov    %rdx,%rax
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rsi
  and    $0xfffffffffffff000,%rsi
  sub    %rsi,%rcx
L184:
  cmp    %rcx,%rsp
  je     L239
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L184
  nopw   0x0(%rax,%rax,1)
L370:
  mov    (%r11),%edx
  mov    %edx,-0x4(%r14,%rax,4)
  mov    %eax,%edx
  add    $0x1,%rax
  cmp    %rax,%r8
  jne    L327
  jmp    L382
L239:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L498
  orq    $0x0,-0x8(%rsp,%rax,1)
L498:
  mov    %rdx,%rax
  mov    %rsp,-0xe8(%rbp)
  mov    %rsp,%rcx
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rsi
  and    $0xfffffffffffff000,%rsi
  sub    %rsi,%rcx
L3:
  cmp    %rcx,%rsp
  je     L33
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L3
L551:
  mov    %edi,%eax
  mov    %rsp,%rcx
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rdx
  shr    $0x4,%rdx
  shl    $0x4,%rdx
  mov    %rdx,%rsi
  and    $0xfffffffffffff000,%rsi
  sub    %rsi,%rcx
L278:
  cmp    %rcx,%rsp
  je     L70
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L278
L33:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L317
L349:
  mov    %rdx,%rax
  mov    %rsp,-0xf8(%rbp)
  mov    %rsp,%rdx
  shr    $0x4,%rax
  shl    $0x4,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  sub    %rcx,%rdx
L479:
  cmp    %rdx,%rsp
  je     L439
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L479
L88:
  orq    $0x0,-0x8(%rsp,%rdx,1)
  jmp    L234
L317:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L349
L70:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L455
  orq    $0x0,-0x8(%rsp,%rdx,1)
L455:
  mov    %rax,%rdx
  mov    %rsp,-0xe8(%rbp)
  mov    %rsp,%rcx
  shr    $0x4,%rdx
  shl    $0x4,%rdx
  mov    %rdx,%rsi
  and    $0xfffffffffffff000,%rsi
  sub    %rsi,%rcx
L411:
  cmp    %rcx,%rsp
  je     L298
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L411
L439:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L354
  orq    $0x0,-0x8(%rsp,%rax,1)
L354:
  mov    %rsp,-0xf0(%rbp)
L510:
  movq   $0x0,-0xd0(%rbp)
  pxor   %xmm4,%xmm4
  movsd  %xmm4,-0xa8(%rbp)
  nop
L384:
  mov    -0xd0(%rbp),%rcx
  mov    -0x48(%rbp),%rdx
  lea    0x0(,%rcx,4),%rax
  mov    (%rdx,%rcx,4),%r13d
  lea    0x4(%rax),%rsi
  mov    %rsi,-0xc8(%rbp)
  cmp    0x4(%rdx,%rax,1),%r13d
  jae    L21
  movsd  -0xa8(%rbp),%xmm3
  movslq %r13d,%rbx
  shl    $0x2,%rbx
  movsd  %xmm3,-0xa0(%rbp)
  movsd  %xmm3,-0x98(%rbp)
  nopl   0x0(%rax)
L4:
  mov    -0x40(%rbp),%rax
  pxor   %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  cvttsd2si (%rax,%rbx,2),%eax
  mov    %eax,-0xc0(%rbp)
  mov    -0x50(%rbp),%rax
  mov    (%rax,%rbx,1),%eax
  lea    0x1(%rax),%edx
  mov    (%r15,%rdx,4),%r8d
  mov    (%r14,%rdx,4),%r12d
  sub    (%r15,%rax,4),%r8d
  sub    (%r14,%rax,4),%r12d
  cvtsi2sd %r8d,%xmm0
  lea    (%r8,%r12,1),%eax
  mov    %r8d,-0xb8(%rbp)
  cvtsi2sd %eax,%xmm1
  mov    D13(%rip),%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm1,-0xb0(%rbp)
  cvttsd2si %xmm0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  pxor   %xmm0,%xmm0
  movsd  -0xb0(%rbp),%xmm1
  mov    D13(%rip),%rax
  cvtsi2sd %r12d,%xmm0
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  mov    -0xb8(%rbp),%r8d
  pxor   %xmm2,%xmm2
  test   %r8d,%r8d
  jle    L20
  movsd  -0xd8(%rbp),%xmm0
  mov    %r8d,%edi
  call   <pow@plt>
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
L20:
  test   %r12d,%r12d
  jg     L520
  pxor   %xmm1,%xmm1
  addsd  -0xa8(%rbp),%xmm2
  add    $0x1,%r13d
  add    $0x4,%rbx
  cvtsi2sdl -0xc0(%rbp),%xmm1
  movsd  -0x98(%rbp),%xmm7
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm7
  movapd %xmm1,%xmm0
  movsd  %xmm1,-0xb0(%rbp)
  movsd  %xmm7,-0x98(%rbp)
  call   <log@plt>
  mov    -0x48(%rbp),%rax
  movsd  -0xb0(%rbp),%xmm1
  mov    -0xc8(%rbp),%rsi
  mulsd  %xmm0,%xmm1
  addsd  -0xa0(%rbp),%xmm1
  movsd  %xmm1,-0xa0(%rbp)
  cmp    %r13d,(%rax,%rsi,1)
  ja     L4
L153:
  movapd %xmm1,%xmm6
L295:
  mov    -0xd0(%rbp),%rax
  pxor   %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  movsd  -0x98(%rbp),%xmm5
  mov    -0xe8(%rbp),%rsi
  mov    -0xf8(%rbp),%rcx
  mov    0x4(%r15,%rax,4),%r13d
  mov    0x4(%r14,%rax,4),%r12d
  sub    (%r15,%rax,4),%r13d
  sub    (%r14,%rax,4),%r12d
  movsd  %xmm5,(%rsi,%rax,8)
  movsd  %xmm6,(%rcx,%rax,8)
  cvtsi2sd %r13d,%xmm0
  lea    0x0(%r13,%r12,1),%eax
  cvtsi2sd %eax,%xmm1
  mov    D13(%rip),%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm1,-0x98(%rbp)
  cvttsd2si %xmm0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  pxor   %xmm0,%xmm0
  movsd  -0x98(%rbp),%xmm1
  mov    D13(%rip),%rax
  cvtsi2sd %r12d,%xmm0
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  test   %r13d,%r13d
  jle    L415
  movsd  -0xd8(%rbp),%xmm0
  mov    %r13d,%edi
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
L415:
  test   %r12d,%r12d
  jg     L525
  addsd  -0xa8(%rbp),%xmm1
  mov    -0xd0(%rbp),%rsi
  mov    -0xf0(%rbp),%rax
  movsd  %xmm1,(%rax,%rsi,8)
  mov    -0x78(%rbp),%eax
  add    $0x1,%rsi
  mov    %rsi,-0xd0(%rbp)
  cmp    %esi,%eax
  ja     L384
L92:
  test   %eax,%eax
  je     L217
  movsd  -0xa8(%rbp),%xmm6
  mov    %eax,%r12d
  xor    %ebx,%ebx
  movapd %xmm6,%xmm3
  movapd %xmm6,%xmm4
  nopl   (%rax)
L509:
  mov    -0xe8(%rbp),%rax
  movsd  (%rax,%rbx,8),%xmm0
  mov    -0xf0(%rbp),%rax
  comisd -0xa8(%rbp),%xmm0
  movsd  (%rax,%rbx,8),%xmm1
  movapd %xmm0,%xmm2
  mulsd  %xmm1,%xmm2
  ja     L68
  mov    -0xf8(%rbp),%rax
  addsd  %xmm2,%xmm3
  mulsd  (%rax,%rbx,8),%xmm1
  add    $0x1,%rbx
  addsd  %xmm1,%xmm4
  cmp    %r12,%rbx
  jne    L509
L56:
  subsd  %xmm4,%xmm6
  movsd  -0xe0(%rbp),%xmm2
  movsd  -0xd8(%rbp),%xmm1
  lea    D3(%rip),%rdi
  mov    $0x3,%eax
  divsd  %xmm3,%xmm6
  movapd %xmm6,%xmm0
  call   <printf@plt>
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L299
  lea    -0x28(%rbp),%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
  nopl   0x0(%rax)
L520:
  movsd  -0xe0(%rbp),%xmm0
  mov    %r12d,%edi
  add    $0x1,%r13d
  add    $0x4,%rbx
  movsd  %xmm2,-0xb0(%rbp)
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  movsd  -0xb0(%rbp),%xmm2
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm1
  movsd  -0x98(%rbp),%xmm5
  cvtsi2sdl -0xc0(%rbp),%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm5
  movapd %xmm1,%xmm0
  movsd  %xmm1,-0xb0(%rbp)
  movsd  %xmm5,-0x98(%rbp)
  call   <log@plt>
  mov    -0x48(%rbp),%rax
  movsd  -0xb0(%rbp),%xmm1
  mov    -0xc8(%rbp),%rcx
  mulsd  %xmm0,%xmm1
  addsd  -0xa0(%rbp),%xmm1
  movsd  %xmm1,-0xa0(%rbp)
  cmp    %r13d,(%rax,%rcx,1)
  ja     L4
  jmp    L153
  nop
L68:
  movsd  %xmm2,-0xc0(%rbp)
  movsd  %xmm1,-0xb8(%rbp)
  movsd  %xmm6,-0xb0(%rbp)
  movsd  %xmm3,-0xa0(%rbp)
  movsd  %xmm4,-0x98(%rbp)
  call   <log@plt>
  movsd  -0xc0(%rbp),%xmm2
  mov    -0xf8(%rbp),%rax
  movapd %xmm0,%xmm5
  movsd  -0xb8(%rbp),%xmm1
  movsd  -0xb0(%rbp),%xmm6
  mulsd  (%rax,%rbx,8),%xmm1
  movsd  -0x98(%rbp),%xmm4
  add    $0x1,%rbx
  mulsd  %xmm2,%xmm5
  movsd  -0xa0(%rbp),%xmm3
  addsd  %xmm2,%xmm3
  addsd  %xmm1,%xmm4
  addsd  %xmm5,%xmm6
  cmp    %rbx,%r12
  jne    L509
  jmp    L56
  nopw   0x0(%rax,%rax,1)
L525:
  movsd  -0xe0(%rbp),%xmm0
  mov    %r12d,%edi
  movsd  %xmm1,-0x98(%rbp)
  call   <pow@plt>
  pxor   %xmm0,%xmm0
  movsd  -0x98(%rbp),%xmm1
  mov    -0xd0(%rbp),%rcx
  cvtsi2sd %eax,%xmm0
  mov    -0xf0(%rbp),%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rcx,8)
  mov    -0x78(%rbp),%eax
  add    $0x1,%rcx
  mov    %rcx,-0xd0(%rbp)
  cmp    %ecx,%eax
  ja     L384
  jmp    L92
  nopw   0x0(%rax,%rax,1)
L21:
  movsd  -0xa8(%rbp),%xmm4
  movsd  %xmm4,-0xa0(%rbp)
  movapd %xmm4,%xmm6
  movsd  %xmm4,-0x98(%rbp)
  jmp    L295
L405:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L78
  orq    $0x0,-0x8(%rsp,%rax,1)
L78:
  mov    %rsp,-0xf0(%rbp)
  test   %edi,%edi
  jne    L510
  pxor   %xmm4,%xmm4
  movsd  %xmm4,-0xa8(%rbp)
L217:
  pxor   %xmm3,%xmm3
  movapd %xmm3,%xmm4
  movapd %xmm3,%xmm6
  jmp    L56
L298:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L496
  orq    $0x0,-0x8(%rsp,%rdx,1)
L496:
  shr    $0x4,%rax
  mov    %rsp,-0xf8(%rbp)
  mov    %rsp,%rdx
  shl    $0x4,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  sub    %rcx,%rdx
L321:
  cmp    %rdx,%rsp
  je     L405
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L321
L299:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_label: D43
type: string
name: parse_data
arguments: data_label: D51
type: string
name: parse_data
arguments: data_label: D13
type: f64
name: parse_data
arguments: data_label: D3
type: string

tool: D43: r

tool: D51: Usage: %s <layer1> <layer2> <overlapping network> <N nodes> <b1> <b2>


tool: D13: 2

tool: D3: %f %f %f


assistant: ```c
int main(int argc, char *argv[]) {
  if (argc < 6) {
    printf("Usage: %s <layer1> <layer2> <overlapping network> <N nodes> <b1> <b2>\n", argv[0]);
    exit(1);
  }
  FILE *filein0, *filein1, *filein, *fileout;
  unsigned int N0, K0, N1, K1, N, K;
  unsigned int *J_slap0, *r_slap0, *J_slap1, *r_slap1, *J_slap, *r_slap;
  double *w_slap;
  int i, j;
  double c_i, d_i, f_i, f_i_2;
  double alpha = (atof(argv[5]));
  double beta = (atof(argv[6]));
  int ov;
  int deg0, deg1;
  double degM, part, f_j, f_j_2;
  double degMrid, maxdegM = 200.0;
  double num1, num2, den, h;
  int number_nodes = (atoi(argv[4]));
  double M = 2.0;
  filein0 = openfile_or_exit(argv[1], "r", 2);
  read_slap(filein0, &K0, &N0, &J_slap0, &r_slap0);
  filein1 = openfile_or_exit(argv[2], "r", 2);
  read_slap(filein1, &K1, &N1, &J_slap1, &r_slap1);
  filein = openfile_or_exit(argv[3], "r", 2);
  read_slap_w(filein, &K, &N, &J_slap, &r_slap, &w_slap);
  int r_slap0_n[N + 1], r_slap1_n[N + 1];
  for (i = 0; i <= N; i++) {
    if (i <= N0) {
      r_slap0_n[i] = r_slap0[i];
    } else {
      r_slap0_n[i] = r_slap0[N0];
    }
    if (i <= N1) {
      r_slap1_n[i] = r_slap1[i];
    } else {
      r_slap1_n[i] = r_slap1[N1];
    }
  }
  double c_i_vec[N];
  double d_i_vec[N];
  double f_i_vec[N];
  for (i = 0; i < N; i++) {
    c_i = 0;
    d_i = 0;
    for (j = r_slap[i]; j < r_slap[i + 1]; j++) {
      ov = w_slap[j];
      deg0 = r_slap0_n[J_slap[j] + 1] - r_slap0_n[J_slap[j]];
      deg1 = r_slap1_n[J_slap[j] + 1] - r_slap1_n[J_slap[j]];
      degM = (deg0 + deg1) * 1.0;
      part = (M / (M - 1)) * (1 - (pow((deg0 / degM), 2)) - (pow((deg1 / degM), 2)));
      if (deg0 > 0.0000000001) {
        f_j = pow(deg0, alpha);
      } else {
        f_j = 0;
      }
      if (deg1 > 0.0000000001) {
        f_j_2 = pow(deg1, beta);
      } else {
        f_j_2 = 0;
      }
      c_i += ov * (f_j + f_j_2);
      d_i += ov * (f_j + f_j_2) * log((ov * (f_j + f_j_2)));
    }
    c_i_vec[i] = c_i;
    d_i_vec[i] = d_i;
    deg0 = r_slap0_n[i + 1] - r_slap0_n[i];
    deg1 = r_slap1_n[i + 1] - r_slap1_n[i];
    degM = (deg0 + deg1) * 1.0;
    part = (M / (M - 1)) * (1 - (pow((deg0 / degM), 2)) - (pow((deg1 / degM), 2)));
    if (deg0 > 0.0000000001) {
      f_i = pow(deg0, alpha);
    } else {
      f_i = 0;
    }
    if (deg1 > 0.0000000001) {
      f_i_2 = pow(deg1, beta);
    } else {
      f_i_2 = 0;
    }
    f_i_vec[i] = f_i + f_i_2;
  }
  num1 = 0;
  num2 = 0;
  den = 0;
  for (i = 0; i < N; i++) {
    if (c_i_vec[i] > 0.0) {
      num1 += f_i_vec[i] * c_i_vec[i] * log(c_i_vec[i]);
    }
    num2 = num2 + f_i_vec[i] * d_i_vec[i];
    den = den + f_i_vec[i] * c_i_vec[i];
  }
  h = (num1 - num2) / den;
  printf("%f %f %f\n", h, alpha, beta);
}
```