user: 提供的汇编代码用C语言表示会是什么样子:

```asm
<r_mvndu>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  sub    $0xc18,%rsp
  mov    0x1c50(%rsp),%rax
  mov    (%r8),%ebx
  mov    (%rdi),%ebp
  mov    %rax,0x50(%rsp)
  mov    0x1c58(%rsp),%rax
  mov    %r8,0x40(%rsp)
  mov    %rax,0x58(%rsp)
  mov    0x1c60(%rsp),%rax
  mov    %r9,0x48(%rsp)
  mov    %rax,0x60(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x1c08(%rsp)
  xor    %eax,%eax
  mov    %ebx,%eax
  neg    %eax
  cmovs  %ebx,%eax
  mov    %eax,0x8(%rsp)
  test   %ebp,%ebp
  jle    L252
  mov    %ebp,%r9d
  mov    %ebp,%r12d
  mov    %ebp,%r15d
  mov    %rsi,%r10
  mov    %rdx,%r13
  mov    %rcx,%rdi
  lea    0x1(%rbp),%r14d
  shr    %r9d
  and    $0xfffffffe,%r12d
  or     $0x1,%r15d
  xor    %r11d,%r11d
  mov    $0x1,%esi
  lea    0x3a8(%rsp),%r8
  xchg   %ax,%ax
L124:
  cmp    $0x1,%ebp
  je     L44
  movslq %r11d,%rax
  lea    (%rdi,%rax,8),%rax
  movupd (%rax),%xmm5
  movups %xmm5,(%r8)
  cmp    $0x1,%r9d
  jbe    L241
  movupd 0x10(%rax),%xmm6
  movups %xmm6,0x10(%r8)
  cmp    $0x2,%r9d
  je     L241
  movupd 0x20(%rax),%xmm7
  movups %xmm7,0x20(%r8)
  cmp    $0x3,%r9d
  je     L241
  movupd 0x30(%rax),%xmm7
  movups %xmm7,0x30(%r8)
  cmp    $0x4,%r9d
  je     L241
  movupd 0x40(%rax),%xmm6
  movups %xmm6,0x40(%r8)
  cmp    $0x5,%r9d
  je     L241
  movupd 0x50(%rax),%xmm5
  movups %xmm5,0x50(%r8)
  cmp    $0x6,%r9d
  je     L241
  movupd 0x60(%rax),%xmm7
  movups %xmm7,0x60(%r8)
  cmp    $0x7,%r9d
  je     L241
  movupd 0x70(%rax),%xmm6
  movups %xmm6,0x70(%r8)
  cmp    $0x8,%r9d
  je     L241
  movupd 0x80(%rax),%xmm5
  movups %xmm5,0x80(%r8)
L241:
  lea    (%r11,%r12,1),%eax
  movslq %r15d,%rdx
  cmp    %ebp,%r12d
  je     L73
L410:
  movslq %esi,%rcx
  cltq
  movsd  (%rdi,%rax,8),%xmm0
  lea    (%rcx,%rcx,4),%rcx
  lea    (%rdx,%rcx,4),%rdx
  movsd  %xmm0,0x300(%rsp,%rdx,8)
L73:
  add    $0x1,%esi
  add    %ebp,%r11d
  add    $0xa0,%r8
  cmp    %r14d,%esi
  jne    L124
  movslq %ebp,%r12
  lea    0x128(%rsp),%rdi
  mov    %r10,%rsi
  lea    0x0(,%r12,8),%r14
  mov    %r14,%rdx
  call   <memcpy@plt>
  lea    0x1c8(%rsp),%rdi
  mov    %r14,%rdx
  mov    %r13,%rsi
  call   <memcpy@plt>
  cmp    0x8(%rsp),%ebp
  jle    L384
  lea    0x300(%rsp),%rax
  lea    -0x1(%rbp),%r8d
  mov    %rax,0x10(%rsp)
  add    %rax,%r14
  lea    (%r8,%r8,4),%rax
  lea    (%r12,%rax,4),%rax
  mov    %r8d,0x4(%rsp)
  lea    0x3a0(%rsp,%rax,8),%rdx
  movslq 0x8(%rsp),%rax
  mov    %rax,%rcx
  sub    %r12,%rcx
  xchg   %ax,%ax
L129:
  movsd  0xa0(%r14),%xmm0
  movsd  0xa0(%r14,%rcx,8),%xmm1
  movsd  %xmm1,0xa0(%r14)
  movsd  %xmm0,0xa0(%r14,%rcx,8)
  add    $0xa0,%r14
  cmp    %rdx,%r14
  jne    L129
  lea    (%r12,%r12,4),%rdi
  lea    (%rax,%rax,4),%rsi
  shl    $0x5,%rdi
  shl    $0x5,%rsi
  lea    0x18(%rdi),%r9
  lea    0x8(%rsi),%rdx
  cmp    %rdx,%r9
  lea    0x8(%rdi),%rcx
  lea    0x18(%rsi),%r9
  setle  %r10b
  cmp    %rcx,%r9
  setle  %r9b
  or     %r10b,%r9b
  je     L260
  cmpl   $0x1,0x4(%rsp)
  jbe    L260
  mov    0x10(%rsp),%rdi
  mov    %ebp,%esi
  shr    %esi
  add    %rdi,%rcx
  add    %rdi,%rdx
  movupd (%rcx),%xmm0
  movupd (%rdx),%xmm6
  movups %xmm6,(%rcx)
  movups %xmm0,(%rdx)
  cmp    $0x1,%esi
  je     L654
  movupd 0x10(%rcx),%xmm0
  movupd 0x10(%rdx),%xmm7
  movups %xmm7,0x10(%rcx)
  movups %xmm0,0x10(%rdx)
  cmp    $0x2,%esi
  je     L654
  movupd 0x20(%rcx),%xmm0
  movupd 0x20(%rdx),%xmm6
  movups %xmm6,0x20(%rcx)
  movups %xmm0,0x20(%rdx)
  cmp    $0x3,%esi
  je     L654
  movupd 0x30(%rcx),%xmm0
  movupd 0x30(%rdx),%xmm5
  movups %xmm5,0x30(%rcx)
  movups %xmm0,0x30(%rdx)
  cmp    $0x4,%esi
  je     L654
  movupd 0x40(%rcx),%xmm0
  movupd 0x40(%rdx),%xmm7
  movups %xmm7,0x40(%rcx)
  movups %xmm0,0x40(%rdx)
  cmp    $0x5,%esi
  je     L654
  movupd 0x50(%rcx),%xmm0
  movupd 0x50(%rdx),%xmm5
  movups %xmm5,0x50(%rcx)
  movups %xmm0,0x50(%rdx)
  cmp    $0x6,%esi
  je     L654
  movupd 0x60(%rcx),%xmm0
  movupd 0x60(%rdx),%xmm6
  movups %xmm6,0x60(%rcx)
  movups %xmm0,0x60(%rdx)
  cmp    $0x7,%esi
  je     L654
  movupd 0x70(%rcx),%xmm0
  movupd 0x70(%rdx),%xmm5
  movups %xmm5,0x70(%rcx)
  movups %xmm0,0x70(%rdx)
  cmp    $0x8,%esi
  je     L654
  movupd 0x80(%rcx),%xmm0
  movupd 0x80(%rdx),%xmm6
  movups %xmm6,0x80(%rcx)
  movups %xmm0,0x80(%rdx)
L654:
  mov    %ebp,%edx
  or     $0x1,%edx
  test   $0x1,%bpl
  je     L210
  movslq %edx,%rdx
  lea    (%r12,%r12,4),%rcx
  lea    (%rax,%rax,4),%rsi
  lea    (%rdx,%rcx,4),%rcx
  lea    (%rdx,%rsi,4),%rdx
  movsd  0x300(%rsp,%rcx,8),%xmm0
  movsd  0x300(%rsp,%rdx,8),%xmm1
  movsd  %xmm1,0x300(%rsp,%rcx,8)
  movsd  %xmm0,0x300(%rsp,%rdx,8)
L210:
  movsd  0x120(%rsp,%r12,8),%xmm0
  movsd  0x120(%rsp,%rax,8),%xmm1
  movsd  %xmm1,0x120(%rsp,%r12,8)
  movsd  0x1c0(%rsp,%rax,8),%xmm1
  movsd  %xmm0,0x120(%rsp,%rax,8)
  movsd  0x1c0(%rsp,%r12,8),%xmm0
  movsd  %xmm1,0x1c0(%rsp,%r12,8)
  movsd  %xmm0,0x1c0(%rsp,%rax,8)
  test   %ebx,%ebx
  jg     L448
L3:
  movsd  0x120(%rsp,%r12,8),%xmm2
  mov    0x4(%rsp),%esi
  test   %esi,%esi
  jg     L599
L12:
  lea    0x120(%rsp),%rsi
  lea    0x1c0(%rsp),%rdx
  nopl   (%rax)
L127:
  lea    0xf80(%rsp),%rcx
  movsd  %xmm2,0x8(%rsp)
  push   0x58(%rsp)
  lea    0x80(%rsp),%rax
  push   %rax
  mov    0x60(%rsp),%rax
  mov    0x14(%rsp),%edi
  movsd  (%rax),%xmm0
  mov    0x58(%rsp),%rax
  lea    0x80(%rsp),%r9
  mov    (%rax),%r8d
  call   <mvn1@plt>
  movsd  0x18(%rsp),%xmm2
  movapd %xmm2,%xmm0
  call   <phi@plt>
  mulsd  0x80(%rsp),%xmm0
  pop    %rax
  pop    %rdx
  mov    0x40(%rsp),%rax
  mov    (%rax),%ecx
  test   %ecx,%ecx
  jns    L76
  xorpd  D9(%rip),%xmm0
L76:
  mov    0x60(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x1c08(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L64
  add    $0x1c18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L44:
  mov    %r11d,%eax
  mov    $0x1,%edx
  jmp    L410
L384:
  lea    -0x1(%rbp),%eax
  mov    %eax,0x4(%rsp)
L264:
  test   %ebx,%ebx
  jle    L3
L448:
  movsd  0x1c0(%rsp,%r12,8),%xmm2
  mov    0x4(%rsp),%esi
  test   %esi,%esi
  jle    L12
L599:
  mov    0x4(%rsp),%ecx
  lea    -0x2(%rbp),%edi
  lea    0x300(%rsp),%r11
  mov    $0xa0,%r8d
  lea    0x2(%rdi),%rbx
  sub    %r11,%r8
  mov    %r11,0x10(%rsp)
  lea    (%r11,%r12,8),%rax
  mov    %rcx,%r15
  mov    %rbx,0x8(%rsp)
  lea    0x0(,%rcx,8),%rbx
  mov    $0xa8,%r9d
  imul   $0xffffffffffffff60,%r12,%rcx
  mov    %r15d,%esi
  mov    %rbx,0x18(%rsp)
  mov    %r15d,%ebx
  or     $0x1,%r15d
  mov    %r8,0x20(%rsp)
  mov    %r9,%r8
  and    $0xfffffffe,%esi
  lea    0x8(,%rdi,8),%rdi
  mov    %r15d,0x6c(%rsp)
  sub    %r11,%r8
  shr    %ebx
  sub    %r11,%rcx
  mov    %rdi,0x38(%rsp)
  movslq %r15d,%rdi
  lea    0x3a0(%rsp),%rdx
  mov    %rcx,0x30(%rsp)
  lea    (%r12,%r12,4),%rcx
  lea    0x80(%rsp),%r12
  mov    %rcx,%r13
  mov    %esi,0x68(%rsp)
  mov    $0x1,%esi
  shl    $0x5,%r13
  mov    %r8,0x28(%rsp)
  lea    (%rdi,%rcx,4),%r8
  lea    0x8(%r11,%r13,1),%r10
  nopw   0x0(%rax,%rax,1)
L274:
  mov    0x18(%rsp),%r15
  movslq %esi,%rcx
  lea    (%r15,%r9,1),%r14
  mov    0x20(%rsp),%r15
  lea    (%r15,%rax,1),%r11
  mov    0x28(%rsp),%r15
  cmp    %r11,%r14
  lea    (%r15,%rax,1),%r11
  setle  %r14b
  cmp    %r9,%r11
  setle  %r11b
  or     %r14b,%r11b
  je     L119
  mov    0x30(%rsp),%r15
  lea    (%r15,%rdx,1),%r11
  cmp    $0x8,%r11
  setne  %r14b
  cmp    $0x2,%ebp
  setne  %r11b
  test   %r11b,%r14b
  je     L119
  movsd  0xa0(%rax),%xmm0
  movupd (%r10),%xmm3
  movupd 0x8(%rdx),%xmm1
  unpcklpd %xmm0,%xmm0
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x8(%rdx)
  cmp    $0x1,%ebx
  jbe    L532
  movupd 0x10(%r10),%xmm3
  movupd 0x18(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x18(%rdx)
  cmp    $0x2,%ebx
  je     L532
  movupd 0x20(%r10),%xmm3
  movupd 0x28(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x28(%rdx)
  cmp    $0x3,%ebx
  je     L532
  movupd 0x30(%r10),%xmm3
  movupd 0x38(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x38(%rdx)
  cmp    $0x4,%ebx
  je     L532
  movupd 0x40(%r10),%xmm3
  movupd 0x48(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x48(%rdx)
  cmp    $0x5,%ebx
  je     L532
  movupd 0x50(%r10),%xmm3
  movupd 0x58(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x58(%rdx)
  cmp    $0x6,%ebx
  je     L532
  movupd 0x60(%r10),%xmm3
  movupd 0x68(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x68(%rdx)
  cmp    $0x7,%ebx
  je     L532
  movupd 0x70(%r10),%xmm3
  movupd 0x78(%rdx),%xmm1
  mulpd  %xmm0,%xmm3
  subpd  %xmm3,%xmm1
  movups %xmm1,0x78(%rdx)
  cmp    $0x8,%ebx
  je     L532
  movupd 0x80(%r10),%xmm1
  mulpd  %xmm0,%xmm1
  movupd 0x88(%rdx),%xmm0
  subpd  %xmm1,%xmm0
  movups %xmm0,0x88(%rdx)
L532:
  mov    0x68(%rsp),%r11d
  cmp    %r11d,0x4(%rsp)
  je     L170
  movsd  0xa0(%rax),%xmm1
  lea    (%rcx,%rcx,4),%rcx
  mulsd  0x300(%rsp,%r8,8),%xmm1
  lea    (%rdi,%rcx,4),%rcx
  movsd  0x300(%rsp,%rcx,8),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x300(%rsp,%rcx,8)
L170:
  movsd  0xa0(%rax),%xmm0
  add    $0xa0,%r9
  add    $0xa0,%rax
  add    $0xa0,%rdx
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,(%r12,%rsi,8)
  add    $0x1,%rsi
  cmp    0x8(%rsp),%rsi
  jne    L274
  mov    0x8(%rsp),%rax
  mov    $0x8,%r14d
  pxor   %xmm4,%xmm4
  lea    0x260(%rsp),%r13
  lea    0x120(%rsp),%rsi
  lea    0x1c0(%rsp),%rdx
  lea    0x0(,%rax,8),%r15
  nopl   0x0(%rax)
L604:
  mov    0x10(%rsp),%rdi
  lea    (%r14,%r14,4),%rax
  lea    (%r14,%rax,4),%rax
  movsd  (%rdi,%rax,1),%xmm0
  ucomisd %xmm0,%xmm4
  ja     L637
  movsd  (%r12,%r14,1),%xmm3
  movsd  (%rsi,%r14,1),%xmm1
  sqrtsd %xmm0,%xmm0
  subsd  %xmm3,%xmm1
  movsd  %xmm0,0x0(%r13,%r14,1)
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsi,%r14,1)
  movsd  (%rdx,%r14,1),%xmm1
  subsd  %xmm3,%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdx,%r14,1)
  add    $0x8,%r14
  cmp    %r15,%r14
  jne    L604
L247:
  mov    $0x1,%r9d
  mov    0x68(%rsp),%r11d
  mov    0x6c(%rsp),%r14d
  lea    0x3a8(%rsp),%r8
  mov    0x8(%rsp),%r10
  mov    0x4(%rsp),%r12d
  movslq %r9d,%rcx
  lea    0x1028(%rsp),%rdi
  movsd  0x0(%r13,%r9,8),%xmm1
  cmp    $0x2,%ebp
  je     L437
  nopl   0x0(%rax,%rax,1)
L151:
  movupd 0x268(%rsp),%xmm4
  movapd %xmm1,%xmm0
  movupd (%r8),%xmm3
  unpcklpd %xmm0,%xmm0
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,(%rdi)
  cmp    $0x1,%ebx
  jbe    L105
  movupd 0x278(%rsp),%xmm4
  movupd 0x10(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x10(%rdi)
  cmp    $0x2,%ebx
  je     L105
  movupd 0x288(%rsp),%xmm4
  movupd 0x20(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x20(%rdi)
  cmp    $0x3,%ebx
  je     L105
  movupd 0x298(%rsp),%xmm4
  movupd 0x30(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x30(%rdi)
  cmp    $0x4,%ebx
  je     L105
  movupd 0x2a8(%rsp),%xmm4
  movupd 0x40(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x40(%rdi)
  cmp    $0x5,%ebx
  je     L105
  movupd 0x2b8(%rsp),%xmm4
  movupd 0x50(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x50(%rdi)
  cmp    $0x6,%ebx
  je     L105
  movupd 0x2c8(%rsp),%xmm4
  movupd 0x60(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x60(%rdi)
  cmp    $0x7,%ebx
  je     L105
  movupd 0x2d8(%rsp),%xmm4
  movupd 0x70(%r8),%xmm3
  mulpd  %xmm0,%xmm4
  divpd  %xmm4,%xmm3
  movups %xmm3,0x70(%rdi)
  cmp    $0x8,%ebx
  je     L105
  movupd 0x2e8(%rsp),%xmm3
  mulpd  %xmm0,%xmm3
  movupd 0x80(%r8),%xmm0
  divpd  %xmm3,%xmm0
  movups %xmm0,0x80(%rdi)
L105:
  movslq %r14d,%rax
  cmp    %r11d,%r12d
  je     L608
L111:
  mulsd  0x260(%rsp,%rax,8),%xmm1
  lea    (%rcx,%rcx,4),%rcx
  lea    (%rax,%rcx,4),%rcx
  movsd  0x300(%rsp,%rcx,8),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0xf80(%rsp,%rcx,8)
L608:
  add    $0x1,%r9
  add    $0xa0,%r8
  add    $0xa0,%rdi
  cmp    %r10,%r9
  je     L127
  movsd  0x0(%r13,%r9,8),%xmm1
  movslq %r9d,%rcx
  cmp    $0x2,%ebp
  jne    L151
L437:
  mov    $0x1,%eax
  jmp    L111
  nop
L119:
  mov    0x38(%rsp),%rcx
  mov    $0x8,%r14d
  sub    %r9,%r14
  lea    (%rcx,%rdx,1),%r15
  mov    %rdx,%rcx
  nopl   (%rax)
L651:
  lea    (%r14,%rcx,1),%r11
  movsd  0x8(%rcx),%xmm0
  add    $0x8,%rcx
  movsd  0x8(%r11,%r13,1),%xmm1
  mulsd  0xa0(%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx)
  cmp    %r15,%rcx
  jne    L651
  jmp    L170
L252:
  lea    -0x1(%rbp),%eax
  movslq %ebp,%r12
  mov    %eax,0x4(%rsp)
  jmp    L264
L260:
  mov    0x10(%rsp),%r10
  lea    (%r12,%r12,4),%rcx
  lea    (%r8,%rcx,4),%rcx
  lea    (%r10,%rdi,1),%rdx
  lea    0x8(%r10,%rcx,8),%r8
  nopl   0x0(%rax,%rax,1)
L580:
  mov    %rdx,%rcx
  movsd  0x8(%rdx),%xmm0
  add    $0x8,%rdx
  sub    %rdi,%rcx
  movsd  0x8(%rcx,%rsi,1),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  %xmm0,0x8(%rcx,%rsi,1)
  cmp    %rdx,%r8
  jne    L580
  jmp    L210
L637:
  mov    %rsi,0x28(%rsp)
  mov    %rdx,0x20(%rsp)
  movsd  %xmm2,0x18(%rsp)
  call   <sqrt@plt>
  mov    0x28(%rsp),%rsi
  mov    0x20(%rsp),%rdx
  pxor   %xmm4,%xmm4
  movsd  (%r12,%r14,1),%xmm3
  movapd %xmm0,%xmm1
  movsd  %xmm0,0x0(%r13,%r14,1)
  movsd  0x18(%rsp),%xmm2
  movsd  (%rsi,%r14,1),%xmm0
  subsd  %xmm3,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rsi,%r14,1)
  movsd  (%rdx,%r14,1),%xmm0
  subsd  %xmm3,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx,%r14,1)
  add    $0x8,%r14
  cmp    %r15,%r14
  jne    L604
  jmp    L247
L64:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D9
type_infered: byte[16]

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
void r_mvndu(int *m0, double *w, double *x, double *rhmat, int *ksign, int *nsim, double *eps, int *ifail, double *der0) {
  double der, tem, bb, pr, sd;
  double rr[20][20], mu[20], a[20], b[20], rc[20][20], dd[20];
  int i, j, m1, k, m, ii;
  void mvn1(int, double[], double[], double[][20], int, double, double *, double *, int *);
  double phi(double);
  m = *m0;
  k = abs(*ksign);
  for (i = 1, ii = 0; i <= m; i++) {
    for (j = 1; j <= m; j++) {
      rr[i][j] = rhmat[ii];
      ii++;
    }
    a[i] = w[i - 1];
    b[i] = x[i - 1];
  }
  if (k < m) {
    for (i = 1; i <= m; i++) {
      tem = rr[i][m];
      rr[i][m] = rr[i][k];
      rr[i][k] = tem;
    }
    for (i = 1; i <= m; i++) {
      tem = rr[m][i];
      rr[m][i] = rr[k][i];
      rr[k][i] = tem;
    }
    tem = a[m];
    a[m] = a[k];
    a[k] = tem;
    tem = b[m];
    b[m] = b[k];
    b[k] = tem;
  }
  m1 = m - 1;
  if ((*ksign) > 0)
    bb = b[m];
  else
    bb = a[m];
  for (i = 1; i <= m1; i++) {
    for (j = 1; j <= m1; j++)
      rr[i][j] -= rr[i][m] * rr[m][j];
    mu[i] = rr[i][m] * bb;
  }
  for (i = 1; i <= m1; i++) {
    dd[i] = sqrt(rr[i][i]);
    a[i] = (a[i] - mu[i]) / dd[i];
    b[i] = (b[i] - mu[i]) / dd[i];
  }
  for (i = 1; i <= m1; i++) {
    for (j = 1; j <= m1; j++)
      rc[i][j] = rr[i][j] / (dd[i] * dd[j]);
  }
  mvn1(m1, a, b, rc, *nsim, *eps, &pr, &sd, ifail);
  der = phi(bb) * pr;
  if (*ksign < 0)
    der = -der;
  *der0 = der;
}
```