user: 请从给定的汇编指令推导出C源代码:

```asm
<target_phi>:
  endbr64
  push   %rbx
  pxor   %xmm0,%xmm0
  mov    %rsi,%rbx
  sub    $0xb0,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xa8(%rsp)
  xor    %eax,%eax
  movaps %xmm0,0x80(%rsp)
  test   %edx,%edx
  jle    L87
  movss  D256(%rip),%xmm5
  cmp    $0x1,%edx
  je     L80
  movss  %xmm5,0x84(%rsp)
  cmp    $0x2,%edx
  je     L80
  movl   $0x3d638e39,0x88(%rsp)
  cmp    $0x3,%edx
  je     L80
  movl   $0x3d638afd,0x8c(%rsp)
L80:
  movss  0x84(%rsp),%xmm6
  pxor   %xmm10,%xmm10
  movss  0x88(%rsp),%xmm2
  movss  0x8c(%rsp),%xmm0
L84:
  movss  D376(%rip),%xmm1
  movss  (%rdi),%xmm4
  movss  0x4(%rdi),%xmm3
  movss  D139(%rip),%xmm7
  movss  D283(%rip),%xmm14
  subss  %xmm1,%xmm4
  movss  D87(%rip),%xmm15
  movss  D72(%rip),%xmm12
  subss  %xmm1,%xmm3
  movss  0x8(%rdi),%xmm1
  subss  D35(%rip),%xmm1
  mulss  %xmm5,%xmm4
  mulss  %xmm6,%xmm3
  mulss  %xmm2,%xmm1
  mulss  0x8(%rbx),%xmm2
  mulss  (%rbx),%xmm5
  addss  %xmm7,%xmm4
  mulss  0x4(%rbx),%xmm6
  addss  %xmm7,%xmm3
  addss  %xmm7,%xmm1
  mulss  %xmm4,%xmm14
  movaps %xmm2,%xmm13
  movss  0xc(%rdi),%xmm2
  subss  D310(%rip),%xmm2
  mulss  %xmm3,%xmm15
  movaps %xmm4,%xmm9
  movaps %xmm3,%xmm11
  mulss  %xmm4,%xmm9
  movaps %xmm1,%xmm8
  mulss  %xmm1,%xmm8
  mulss  %xmm0,%xmm2
  mulss  0xc(%rbx),%xmm0
  movss  %xmm14,0x14(%rsp)
  movaps %xmm1,%xmm14
  movss  %xmm15,0x10(%rsp)
  mulss  %xmm3,%xmm11
  movss  D457(%rip),%xmm15
  mulss  %xmm4,%xmm12
  mulss  %xmm8,%xmm14
  mulss  D67(%rip),%xmm14
  addss  %xmm7,%xmm2
  movss  %xmm0,0x2c(%rsp)
  movss  D84(%rip),%xmm0
  mulss  %xmm1,%xmm15
  mulss  %xmm4,%xmm0
  movaps %xmm2,%xmm7
  mulss  %xmm2,%xmm7
  movss  %xmm14,0x1c(%rsp)
  movaps %xmm4,%xmm14
  mulss  %xmm9,%xmm14
  movss  %xmm15,0x18(%rsp)
  movss  D258(%rip),%xmm15
  movss  %xmm0,0xc(%rsp)
  movss  D114(%rip),%xmm0
  mulss  %xmm3,%xmm15
  mulss  %xmm14,%xmm0
  movss  D363(%rip),%xmm14
  mulss  %xmm11,%xmm14
  movss  %xmm15,0x30(%rsp)
  mulss  %xmm1,%xmm15
  movss  %xmm0,0x20(%rsp)
  movss  %xmm14,0x38(%rsp)
  mulss  %xmm1,%xmm14
  movss  %xmm15,0x34(%rsp)
  mulss  %xmm2,%xmm15
  movss  %xmm14,0x3c(%rsp)
  movss  D200(%rip),%xmm0
  mulss  %xmm2,%xmm14
  addss  %xmm12,%xmm0
  movss  D405(%rip),%xmm12
  mulss  %xmm1,%xmm12
  addss  %xmm0,%xmm12
  movss  0xc(%rsp),%xmm0
  mulss  %xmm3,%xmm0
  addss  %xmm12,%xmm0
  movss  0x10(%rsp),%xmm12
  mulss  %xmm1,%xmm12
  addss  %xmm0,%xmm12
  movss  0x14(%rsp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm12,%xmm0
  movss  0x18(%rsp),%xmm12
  mulss  %xmm2,%xmm12
  addss  %xmm0,%xmm12
  movss  D416(%rip),%xmm0
  mulss  %xmm3,%xmm0
  addss  %xmm12,%xmm0
  movss  D401(%rip),%xmm12
  mulss  %xmm2,%xmm12
  addss  %xmm0,%xmm12
  movss  D351(%rip),%xmm0
  mulss  %xmm4,%xmm0
  mulss  %xmm7,%xmm0
  addss  %xmm12,%xmm0
  movss  0x1c(%rsp),%xmm12
  mulss  %xmm2,%xmm12
  addss  %xmm0,%xmm12
  movss  D89(%rip),%xmm0
  mulss  %xmm4,%xmm0
  mulss  %xmm11,%xmm0
  addss  %xmm12,%xmm0
  movss  D445(%rip),%xmm12
  mulss  %xmm1,%xmm12
  addss  %xmm15,%xmm0
  movss  D102(%rip),%xmm15
  mulss  %xmm3,%xmm15
  mulss  %xmm7,%xmm12
  movss  %xmm15,0x40(%rsp)
  mulss  %xmm1,%xmm15
  addss  %xmm0,%xmm12
  movss  0x20(%rsp),%xmm0
  mulss  %xmm3,%xmm0
  addss  %xmm12,%xmm0
  movss  D306(%rip),%xmm12
  mulss  %xmm4,%xmm12
  addss  %xmm14,%xmm0
  movaps %xmm12,%xmm14
  mulss  %xmm3,%xmm14
  movss  %xmm14,0x24(%rsp)
  movss  D301(%rip),%xmm14
  mulss  %xmm8,%xmm14
  movss  %xmm14,0x28(%rsp)
  movaps %xmm15,%xmm14
  mulss  %xmm7,%xmm14
  addss  %xmm0,%xmm14
  movss  0x24(%rsp),%xmm0
  mulss  %xmm1,%xmm0
  addss  %xmm14,%xmm0
  movss  0x28(%rsp),%xmm14
  mulss  %xmm2,%xmm14
  addss  %xmm14,%xmm0
  movss  D108(%rip),%xmm14
  mulss  %xmm9,%xmm14
  movaps %xmm14,%xmm15
  mulss  %xmm3,%xmm15
  movss  %xmm15,0x44(%rsp)
  mulss  %xmm1,%xmm15
  addss  %xmm0,%xmm15
  movaps %xmm5,%xmm0
  addss  %xmm10,%xmm0
  addss  %xmm6,%xmm0
  addss  %xmm13,%xmm0
  addss  0x2c(%rsp),%xmm0
  ucomiss %xmm10,%xmm0
  jp     L47
  movaps %xmm10,%xmm0
  jne    L47
L67:
  movss  %xmm0,(%rbx)
  mov    0xa8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L31
  add    $0xb0,%rsp
  movaps %xmm15,%xmm0
  pop    %rbx
  ret
  nopl   (%rax)
L47:
  movss  %xmm6,0x4c(%rsp)
  mulss  %xmm5,%xmm5
  movss  D283(%rip),%xmm6
  movss  D84(%rip),%xmm0
  mulss  D89(%rip),%xmm11
  mov    $0x2,%edi
  movss  %xmm15,0x7c(%rsp)
  mulss  %xmm2,%xmm6
  movss  %xmm8,0x78(%rsp)
  mulss  %xmm3,%xmm0
  movss  %xmm13,0x74(%rsp)
  subss  D45(%rip),%xmm0
  movss  %xmm10,0x60(%rsp)
  mulss  D233(%rip),%xmm9
  movss  %xmm5,0x48(%rsp)
  movss  %xmm14,0x70(%rsp)
  movaps %xmm6,%xmm5
  movaps %xmm4,%xmm6
  movss  %xmm12,0x6c(%rsp)
  addss  %xmm0,%xmm5
  movss  D351(%rip),%xmm0
  mulss  D181(%rip),%xmm6
  movss  %xmm2,0x5c(%rsp)
  mulss  %xmm3,%xmm9
  movss  %xmm7,0x68(%rsp)
  mulss  %xmm7,%xmm0
  movss  %xmm4,0x58(%rsp)
  movss  %xmm3,0x54(%rsp)
  movss  %xmm1,0x50(%rsp)
  addss  %xmm5,%xmm0
  movaps %xmm6,%xmm5
  mulss  %xmm3,%xmm5
  addss  %xmm0,%xmm11
  movss  D306(%rip),%xmm0
  mulss  %xmm3,%xmm0
  addss  %xmm11,%xmm9
  mulss  %xmm1,%xmm5
  mulss  %xmm1,%xmm0
  addss  %xmm9,%xmm0
  addss  %xmm5,%xmm0
  call   <pow@plt>
  movss  0x50(%rsp),%xmm1
  mov    $0x2,%edi
  movss  D87(%rip),%xmm5
  movss  0x60(%rsp),%xmm10
  mulss  0x48(%rsp),%xmm0
  movss  0xc(%rsp),%xmm4
  mulss  %xmm1,%xmm5
  movss  0x54(%rsp),%xmm3
  movss  0x5c(%rsp),%xmm2
  movss  0x68(%rsp),%xmm7
  movss  %xmm10,0x64(%rsp)
  addss  %xmm10,%xmm4
  movss  0x6c(%rsp),%xmm12
  movss  0x70(%rsp),%xmm14
  movss  %xmm3,0x5c(%rsp)
  movss  0x4c(%rsp),%xmm6
  movss  %xmm7,0x54(%rsp)
  movss  %xmm0,0x48(%rsp)
  mulss  %xmm1,%xmm12
  mulss  %xmm1,%xmm14
  movaps %xmm4,%xmm0
  movss  0x58(%rsp),%xmm4
  movss  %xmm2,0x58(%rsp)
  addss  %xmm5,%xmm0
  movss  D431(%rip),%xmm5
  subss  D37(%rip),%xmm0
  mulss  %xmm6,%xmm6
  movss  %xmm4,0x60(%rsp)
  mulss  %xmm3,%xmm5
  movss  %xmm6,0x4c(%rsp)
  mulss  %xmm4,%xmm5
  movss  D258(%rip),%xmm4
  mulss  %xmm1,%xmm4
  addss  %xmm0,%xmm5
  movaps %xmm4,%xmm0
  movaps %xmm3,%xmm4
  mulss  %xmm2,%xmm0
  mulss  D130(%rip),%xmm4
  addss  %xmm5,%xmm0
  movss  0x20(%rsp),%xmm5
  addss  %xmm0,%xmm5
  movaps %xmm4,%xmm0
  mulss  %xmm1,%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm5,%xmm0
  movss  D102(%rip),%xmm5
  mulss  %xmm1,%xmm5
  mulss  %xmm7,%xmm5
  addss  %xmm0,%xmm5
  addss  %xmm5,%xmm12
  addss  %xmm14,%xmm12
  movaps %xmm12,%xmm0
  call   <pow@plt>
  mulss  0x4c(%rsp),%xmm0
  movss  0x48(%rsp),%xmm7
  mov    $0x2,%edi
  movss  0x58(%rsp),%xmm2
  movss  0x78(%rsp),%xmm8
  mulss  D422(%rip),%xmm8
  movss  D405(%rip),%xmm5
  addss  0x10(%rsp),%xmm5
  movss  0x50(%rsp),%xmm1
  movss  0x40(%rsp),%xmm6
  movss  0x74(%rsp),%xmm13
  addss  %xmm0,%xmm7
  movss  D457(%rip),%xmm0
  mulss  %xmm13,%xmm13
  mulss  %xmm2,%xmm0
  mulss  %xmm2,%xmm8
  movss  %xmm7,0xc(%rsp)
  movss  0x54(%rsp),%xmm7
  mulss  %xmm7,%xmm6
  movss  %xmm13,0x20(%rsp)
  addss  %xmm5,%xmm0
  movss  0x38(%rsp),%xmm5
  mulss  %xmm2,%xmm5
  addss  %xmm0,%xmm8
  movss  0x30(%rsp),%xmm0
  movss  %xmm1,0x30(%rsp)
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm8
  movss  D445(%rip),%xmm0
  mulss  %xmm7,%xmm0
  addss  %xmm8,%xmm0
  addss  %xmm0,%xmm5
  movaps %xmm1,%xmm0
  mulss  D44(%rip),%xmm0
  addss  %xmm6,%xmm5
  addss  0x24(%rsp),%xmm5
  movss  %xmm2,0x24(%rsp)
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm5
  movss  0x44(%rsp),%xmm0
  addss  %xmm5,%xmm0
  call   <pow@plt>
  movss  0x24(%rsp),%xmm2
  movss  D149(%rip),%xmm5
  mov    $0x2,%edi
  mulss  0x20(%rsp),%xmm0
  movss  0x64(%rsp),%xmm10
  mulss  %xmm2,%xmm5
  addss  0xc(%rsp),%xmm0
  movss  0x14(%rsp),%xmm7
  movss  0x60(%rsp),%xmm4
  movss  0x5c(%rsp),%xmm3
  addss  %xmm10,%xmm7
  movss  0x30(%rsp),%xmm1
  movss  0x2c(%rsp),%xmm6
  movss  %xmm0,0xc(%rsp)
  movss  0x18(%rsp),%xmm0
  mulss  %xmm6,%xmm6
  mulss  %xmm5,%xmm4
  addss  %xmm7,%xmm0
  addss  D401(%rip),%xmm0
  movss  %xmm6,0x10(%rsp)
  addss  %xmm4,%xmm0
  movss  D111(%rip),%xmm4
  addss  0x1c(%rsp),%xmm0
  addss  0x34(%rsp),%xmm0
  mulss  %xmm2,%xmm4
  mulss  D8(%rip),%xmm2
  mulss  %xmm1,%xmm4
  mulss  %xmm3,%xmm2
  addss  %xmm4,%xmm0
  addss  0x3c(%rsp),%xmm0
  mulss  %xmm1,%xmm2
  addss  %xmm2,%xmm0
  addss  0x28(%rsp),%xmm0
  call   <pow@plt>
  mulss  0x10(%rsp),%xmm0
  addss  0xc(%rsp),%xmm0
  call   <sqrt@plt>
  movss  0x7c(%rsp),%xmm15
  jmp    L67
  nopl   0x0(%rax,%rax,1)
L87:
  pxor   %xmm10,%xmm10
  movaps %xmm10,%xmm0
  movaps %xmm10,%xmm2
  movaps %xmm10,%xmm6
  movaps %xmm10,%xmm5
  jmp    L84
L31:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_label: D256
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D376
type_guessed: f32
name: retrieve_data
arguments: data_label: D139
type_guessed: f32
name: retrieve_data
arguments: data_label: D283
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D87
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D72
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D35
type_guessed: f32
name: retrieve_data
arguments: data_label: D310
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D457
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D67
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D84
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D258
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D114
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D363
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D200
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D405
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D416
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D401
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D351
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D89
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D445
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D102
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D306
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D301
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D108
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D45
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D233
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D181
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D431
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D37
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D130
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D422
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D44
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D149
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D111
type_guessed: byte[4]
name: retrieve_data
arguments: data_label: D8
type_guessed: byte[4]

tool: {"D256": [0x89,0x88,0x88,0x3d]}

tool: D376: 15

tool: D139: 1

tool: {"D283": [0xd3,0xb7,0x9d,0xbb]}

tool: {"D87": [0xdd,0xac,0xaa,0xbb]}

tool: {"D72": [0xa4,0xd6,0x69,0xbd]}

tool: D35: 18

tool: {"D310": [0x79,0xe9,0x8f,0x41]}

tool: {"D457": [0x34,0xf5,0x0e,0x3c]}

tool: {"D67": [0x56,0x1e,0xc0,0xb8]}

tool: {"D84": [0x5a,0x2a,0x3b,0x3b]}

tool: {"D258": [0x11,0x38,0x02,0xb9]}

tool: {"D114": [0x9c,0x88,0x1b,0xb8]}

tool: {"D363": [0x9a,0x1b,0x91,0xb8]}

tool: {"D200": [0x8c,0xca,0x97,0x3a]}

tool: {"D405": [0x89,0xe4,0x8c,0x3d]}

tool: {"D416": [0x70,0x85,0x25,0xb9]}

tool: {"D401": [0x9c,0x97,0x89,0x39]}

tool: {"D351": [0x0c,0xdd,0x72,0xb8]}

tool: {"D89": [0xd4,0xb9,0xf7,0x37]}

tool: {"D445": [0xdd,0x80,0x39,0x39]}

tool: {"D102": [0x81,0x4e,0xbe,0x38]}

tool: {"D306": [0x52,0xbb,0x0e,0x37]}

tool: {"D301": [0x5a,0x3e,0x7e,0xb7]}

tool: {"D108": [0xc4,0x8f,0xbe,0x38]}

tool: {"D45": [0xa4,0xd6,0x69,0x3d]}

tool: {"D233": [0xea,0x4c,0xe9,0xb8]}

tool: {"D181": [0xc4,0x8f,0x3e,0x39]}

tool: {"D431": [0xd4,0xb9,0x77,0x38]}

tool: {"D37": [0x70,0x85,0x25,0x39]}

tool: {"D130": [0x9a,0x1b,0x11,0xb9]}

tool: {"D422": [0xc0,0x16,0x90,0xb9]}

tool: {"D44": [0x5a,0x3e,0xfe,0xb7]}

tool: {"D149": [0x0c,0xdd,0xf2,0xb8]}

tool: {"D111": [0xdd,0x80,0xb9,0x39]}

tool: {"D8": [0x81,0x4e,0x3e,0x39]}

assistant: ```c
float target_phi(float *x, float *dx, int m) {
  int ncoeff = 20;
  float avdat = 0.1158064E-02;
  float xmin[10] = {-0.15000E+02, -0.15000E+02, -0.18000E+02, -0.18013E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00};
  float xmax[10] = {0.15000E+02, 0.15000E+02, 0.18000E+02, 0.17989E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00};
  float scale[10] = {0};
  float coeff[21] = {0.10983568E-07, -0.57089463E-01, 0.68795271E-01, 0.28559179E-02, -0.52085952E-02, -0.48131733E-02, 0.87254532E-02, -0.15785336E-03, 0.26243634E-03, -0.57903162E-04, -0.91609239E-04, 0.29531227E-04, -0.12418652E-03, 0.17690980E-03, -0.37082107E-04, -0.69192800E-04, 0.90745285E-04, 0.85074735E-05, -0.15154097E-04, 0.90866844E-04, 0.};
  int ientry = 0;
  int i;
  if (ientry == 0) {
    ientry = 1;
    for (i = 0; i < m; i++) {
      if (xmin[i] == xmax[i])
        continue;
      scale[i] = 2. / (xmax[i] - xmin[i]);
    }
  }
  float x1 = 1. + (x[0] - xmax[0]) * scale[0];
  float dx1 = dx[0] * scale[0];
  float x2 = 1. + (x[1] - xmax[1]) * scale[1];
  float dx2 = dx[1] * scale[1];
  float x3 = 1. + (x[2] - xmax[2]) * scale[2];
  float dx3 = dx[2] * scale[2];
  float x4 = 1. + (x[3] - xmax[3]) * scale[3];
  float dx4 = dx[3] * scale[3];
  float x11 = x1;
  float x12 = x11 * x1;
  float x13 = x12 * x1;
  float x21 = x2;
  float x22 = x21 * x2;
  float x31 = x3;
  float x32 = x31 * x3;
  float x33 = x32 * x3;
  float x41 = x4;
  float x42 = x41 * x4;
  float v_target_phi = avdat + coeff[0] + coeff[1] * x11 + coeff[2] * x31 + coeff[3] * x11 * x21 + coeff[4] * x21 * x31 + coeff[5] * x11 * x41 + coeff[6] * x31 * x41 + coeff[7] * x21;
  v_target_phi = v_target_phi + coeff[8] * x41 + coeff[9] * x11 * x42 + coeff[10] * x33 * x41 + coeff[11] * x11 * x22 + coeff[12] * x21 * x31 * x41 + coeff[13] * x31 * x42 + coeff[14] * x13 * x21 + coeff[15] * x22 * x31 * x41 + coeff[16] * x21 * x31 * x42;
  v_target_phi = v_target_phi + coeff[17] * x11 * x21 * x31 + coeff[18] * x32 * x41 + coeff[19] * x12 * x21 * x31;
  float dxall = 0 + dx1 + dx2 + dx3 + dx4;
  float dv_target_phi = dxall == 0 ? 0 : sqrt(+dx1 * dx1 * pow(0 + coeff[1] + coeff[3] * x21 + coeff[5] * x41 + coeff[9] * x42 + coeff[11] * x22 + 3 * coeff[14] * x12 * x21 + coeff[17] * x21 * x31 + 2 * coeff[19] * x11 * x21 * x31, 2) + dx2 * dx2 * pow(0 + coeff[3] * x11 + coeff[4] * x31 + coeff[7] + 2 * coeff[11] * x21 * x11 + coeff[12] * x31 * x41 + coeff[14] * x13 + 2 * coeff[15] * x21 * x31 * x41 + coeff[16] * x31 * x42 + coeff[17] * x11 * x31 + coeff[19] * x12 * x31, 2) + dx3 * dx3 * pow(0 + coeff[2] + coeff[4] * x21 + coeff[6] * x41 + 3 * coeff[10] * x32 * x41 + coeff[12] * x21 * x41 + coeff[13] * x42 + coeff[15] * x22 * x41 + coeff[16] * x21 * x42 + coeff[17] * x11 * x21 + 2 * coeff[18] * x31 * x41 + coeff[19] * x12 * x21, 2) + dx4 * dx4 * pow(0 + coeff[5] * x11 + coeff[6] * x31 + coeff[8] + 2 * coeff[9] * x41 * x11 + coeff[10] * x33 + coeff[12] * x21 * x31 + 2 * coeff[13] * x41 * x31 + coeff[15] * x22 * x31 + 2 * coeff[16] * x41 * x21 * x31 + coeff[18] * x32, 2));
  dx[0] = dv_target_phi;
  return v_target_phi;
}
```