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
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %edi,%ebp
  mov    %fs:0x28,%rax
  mov    %rax,0x18(%rsp)
  xor    %eax,%eax
  call   <prompt_for_artificial@plt>
  cmp    $0x1,%eax
  je     L5
  mov    $0x0,%eax
  call   <read_quaddobl_target_system@plt>
  mov    $0x0,%eax
  call   <copy_quaddobl_target_solutions_to_container@plt>
  mov    $0x0,%eax
  call   <prompt_for_parameter@plt>
  mov    %eax,%ebx
  mov    %eax,%edi
  call   <solcon_quaddobl_drop_coordinate_by_index@plt>
  mov    $0x0,%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
  lea    D11(%rip),%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  test   %ebx,%ebx
  je     L11
  mov    $0x1,%esi
  mov    %ebx,%edi
  call   <padcon_quaddobl_initialize_parameter_homotopy@plt>
  jmp    L6
L5:
  mov    $0x0,%eax
  call   <read_quaddobl_target_system@plt>
  mov    $0x0,%eax
  call   <read_quaddobl_start_system@plt>
  mov    $0x0,%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
  lea    D11(%rip),%rdi
  mov    $0x0,%eax
  call   <printf@plt>
L11:
  mov    %ebp,%esi
  mov    $0x1,%edi
  call   <padcon_quaddobl_initialize_homotopy@plt>
L6:
  mov    $0x1,%ebx
  lea    D35(%rip),%r12
  lea    D3(%rip),%r15
  lea    0xf(%rsp),%r14
  lea    D26(%rip),%r13
  cmpl   $0x0,0x10(%rsp)
  jle    L22
L33:
  mov    %ebx,%esi
  mov    %r12,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  mov    %ebx,%edi
  call   <quaddobl_next_loop@plt>
  cmp    %ebx,0x10(%rsp)
  je     L22
  mov    %r15,%rdi
  mov    $0x0,%eax
  call   <printf@plt>
  mov    $0x0,%eax
  call   <getchar@plt>
  mov    %eax,%ebp
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <scanf@plt>
  cmp    $0x79,%bpl
  jne    L22
  add    $0x1,%ebx
  cmp    %ebx,0x10(%rsp)
  jge    L33
L22:
  mov    0x18(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L58
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L58:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: label: D11
type_guessed: string
name: read_data
arguments: label: D35
type_guessed: string
name: read_data
arguments: label: D3
type_guessed: string
name: read_data
arguments: label: D26
type_guessed: string

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


tool: D35: 
Tracking path %d ...


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

tool: D26: %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;
  }
}
```