user: What would the C source look like for the assembly code provided here:

```asm
<nbodyC>:
  endbr64
  push   %r15
  pxor   %xmm0,%xmm0
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdx,0x30(%rsp)
  mov    %ecx,0x4c(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  movaps %xmm0,0x50(%rsp)
  movq   $0x0,0x60(%rsp)
  test   %ecx,%ecx
  jle    L39
  lea    0x50(%rsp),%rax
  mov    %rdi,0x38(%rsp)
  mov    %rsi,%rbx
  xor    %r15d,%r15d
  mov    %rax,0x40(%rsp)
  movslq %ecx,%rax
  mov    %rax,0x8(%rsp)
L26:
  mov    0x40(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  xor    %r13d,%r13d
  call   <memset@plt>
  lea    (%r15,%r15,2),%eax
  shl    $0x3,%rax
  lea    0x8(%rbx,%rax,1),%rcx
  lea    (%rbx,%rax,1),%r14
  mov    %rcx,0x28(%rsp)
  lea    0x10(%rbx,%rax,1),%r12
  jmp    L80
  nopl   0x0(%rax,%rax,1)
L44:
  add    $0x1,%r13
  cmp    0x8(%rsp),%r13
  je     L9
L80:
  cmp    %r13d,%r15d
  je     L44
  lea    0x0(%r13,%r13,2),%edx
  mov    $0x2,%edi
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rbp
  lea    (%rbx,%rbp,1),%r8
  movsd  (%r8),%xmm0
  subsd  (%r14),%xmm0
  mov    %r8,0x20(%rsp)
  call   <pow@plt>
  mov    0x28(%rsp),%rax
  mov    $0x2,%edi
  movsd  %xmm0,0x10(%rsp)
  movsd  0x8(%rbx,%rbp,1),%xmm0
  subsd  (%rax),%xmm0
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm4
  lea    0x10(%rbx,%rbp,1),%rdx
  mov    $0x2,%edi
  mov    %rdx,0x18(%rsp)
  addsd  %xmm0,%xmm4
  movsd  (%rdx),%xmm0
  subsd  (%r12),%xmm0
  movsd  %xmm4,0x10(%rsp)
  call   <pow@plt>
  addsd  0x10(%rsp),%xmm0
  mov    $0x1,%edi
  call   <pow@plt>
  pxor   %xmm6,%xmm6
  mov    0x18(%rsp),%rdx
  mov    0x20(%rsp),%r8
  ucomisd %xmm6,%xmm0
  movapd %xmm0,%xmm3
  unpcklpd %xmm3,%xmm3
  jp     L67
  je     L63
L67:
  movupd (%r8),%xmm1
  movupd (%r14),%xmm7
  mov    0x30(%rsp),%rax
  subpd  %xmm7,%xmm1
  movsd  (%rax,%r13,8),%xmm2
  add    $0x1,%r13
  divpd  %xmm3,%xmm1
  movapd %xmm2,%xmm3
  unpcklpd %xmm3,%xmm3
  mulpd  %xmm3,%xmm1
  addpd  0x50(%rsp),%xmm1
  movaps %xmm1,0x50(%rsp)
  movsd  (%rdx),%xmm1
  subsd  (%r12),%xmm1
  divsd  %xmm0,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  0x60(%rsp),%xmm1
  movsd  %xmm1,0x60(%rsp)
  cmp    0x8(%rsp),%r13
  jne    L80
L9:
  mov    0x38(%rsp),%rax
  movapd 0x50(%rsp),%xmm7
  add    $0x1,%r15d
  movsd  0x60(%rsp),%xmm0
  movups %xmm7,(%rax)
  add    $0x18,%rax
  movsd  %xmm0,-0x8(%rax)
  mov    %rax,0x38(%rsp)
  cmp    %r15d,0x4c(%rsp)
  jne    L26
L39:
  xor    %eax,%eax
  jmp    L2
  xchg   %ax,%ax
L63:
  mov    $0xffffffff,%eax
L2:
  mov    0x68(%rsp),%rdx
  sub    %fs:0x28,%rdx
  jne    L73
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L73:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
int nbodyC(double dx[], double x[], double mus[], int n) {
  double temp[3] = {0.0};
  double rnorm3;
  int j;
  int k;
  int i;
  for (j = 0; j < n; j++) {
    memset(temp, 0, sizeof(temp));
    for (k = 0; k < n; k++) {
      if (k == j) {
        continue;
      }
      rnorm3 = pow(pow(x[k * 3] - x[j * 3], 2.) + pow(x[k * 3 + 1] - x[j * 3 + 1], 2.) + pow(x[k * 3 + 2] - x[j * 3 + 2], 2.), 3. / 2.);
      if (rnorm3 == 0.0) {
        return -1;
      }
      for (i = 0; i < 3; i++) {
        temp[i] += (x[k * 3 + i] - x[j * 3 + i]) / rnorm3 * mus[k];
      }
    }
    for (i = 0; i < 3; i++) {
      dx[j * 3 + i] = temp[i];
    }
  }
  return 0;
}
```