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

```asm
<quaddobl_next_step>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %edi,%r12d
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x18(%rsp)
  xor    %eax,%eax
  call   <prompt_for_artificial@plt>
  cmp    $0x1,%eax
  mov    $0x0,%eax
  je     L36
  call   <read_quaddobl_target_system@plt>
  xor    %eax,%eax
  call   <copy_quaddobl_target_solutions_to_container@plt>
  xor    %eax,%eax
  call   <prompt_for_parameter@plt>
  mov    %eax,%edi
  mov    %eax,%ebp
  call   <solcon_quaddobl_drop_coordinate_by_index@plt>
  xor    %eax,%eax
  call   <solcon_quaddobl_set_continuation_parameter@plt>
  lea    0x10(%rsp),%rdi
  call   <solcon_number_of_quaddobl_solutions@plt>
  lea    0x14(%rsp),%rdi
  call   <solcon_dimension_of_quaddobl_solutions@plt>
  mov    0x14(%rsp),%edx
  mov    0x10(%rsp),%esi
  xor    %eax,%eax
  lea    D4(%rip),%rdi
  call   <printf@plt>
  test   %ebp,%ebp
  je     L57
  mov    $0x1,%esi
  mov    %ebp,%edi
  call   <padcon_quaddobl_initialize_parameter_homotopy@plt>
L47:
  mov    0x10(%rsp),%eax
  mov    $0x1,%ebx
  lea    D26(%rip),%r12
  lea    D40(%rip),%r15
  lea    0xf(%rsp),%r14
  lea    D39(%rip),%r13
  test   %eax,%eax
  jg     L10
  jmp    L39
  nopw   0x0(%rax,%rax,1)
L7:
  mov    %r15,%rdi
  xor    %eax,%eax
  call   <printf@plt>
  xor    %eax,%eax
  call   <getchar@plt>
  mov    %r14,%rsi
  mov    %r13,%rdi
  mov    %eax,%ebp
  call   <scanf@plt>
  cmp    $0x79,%bpl
  jne    L39
  add    $0x1,%ebx
  cmp    %ebx,0x10(%rsp)
  jl     L39
L10:
  mov    %r12,%rdi
  xor    %eax,%eax
  mov    %ebx,%esi
  call   <printf@plt>
  mov    %ebx,%edi
  call   <quaddobl_next_loop@plt>
  cmp    %ebx,0x10(%rsp)
  jne    L7
L39:
  mov    0x18(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L72
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L36:
  call   <read_quaddobl_target_system@plt>
  xor    %eax,%eax
  call   <read_quaddobl_start_system@plt>
  xor    %eax,%eax
  call   <copy_quaddobl_start_solutions_to_container@plt>
  lea    0x10(%rsp),%rdi
  call   <solcon_number_of_quaddobl_solutions@plt>
  lea    0x14(%rsp),%rdi
  call   <solcon_dimension_of_quaddobl_solutions@plt>
  mov    0x14(%rsp),%edx
  mov    0x10(%rsp),%esi
  xor    %eax,%eax
  lea    D4(%rip),%rdi
  call   <printf@plt>
L57:
  mov    %r12d,%esi
  mov    $0x1,%edi
  call   <padcon_quaddobl_initialize_homotopy@plt>
  jmp    L47
L72:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: data_tag: D4
type_infered: string
name: disassemble_data
arguments: data_tag: D26
type_infered: string
name: disassemble_data
arguments: data_tag: D40
type_infered: string
name: disassemble_data
arguments: data_tag: D39
type_infered: string

tool: D4: Read %d start solutions of dimension %d.


tool: D26: 
Tracking path %d ...


tool: D40: Continue with the next solution ? (y/n) 

tool: D39: %c

assistant: ```c
void quaddobl_next_step(int homo) {
  int fail, length, index, failed, strlensol, dim, arthom;
  int idx = 0;
  char contstep, nlsb;
  double frp, re_pole, im_pole, tval, step;
  arthom = prompt_for_artificial();
  if (arthom == 1) {
    fail = read_quaddobl_target_system();
    fail = read_quaddobl_start_system();
    fail = copy_quaddobl_start_solutions_to_container();
  } else {
    fail = read_quaddobl_target_system();
    fail = copy_quaddobl_target_solutions_to_container();
    idx = prompt_for_parameter();
    fail = solcon_quaddobl_drop_coordinate_by_index(idx);
    fail = solcon_quaddobl_set_continuation_parameter();
  }
  fail = solcon_number_of_quaddobl_solutions(&length);
  fail = solcon_dimension_of_quaddobl_solutions(&dim);
  printf("Read %d start solutions of dimension %d.\n", length, dim);
  if (idx == 0)
    fail = padcon_quaddobl_initialize_homotopy(1, homo);
  else
    fail = padcon_quaddobl_initialize_parameter_homotopy(idx, 1);
  for (index = 1; index <= length; index++) {
    printf("\nTracking path %d ...\n", index);
    quaddobl_next_loop(index);
    if (index == length)
      break;
    printf("Continue with the next solution ? (y/n) ");
    contstep = getchar();
    scanf("%c", &nlsb);
    if (contstep != 'y')
      break;
  }
}
```