user: Demonstrate the C code that generates the assembly output below:

```asm
<target_phi>:
  endbr64
  push   %r14
  push   %rbx
  sub    $0x108,%rsp
  mov    %rdi,%rcx
  mov    %rsi,%rbx
  mov    %fs:0x28,%rax
  mov    %rax,0xf8(%rsp)
  xor    %eax,%eax
  movss  D142(%rip),%xmm0
  movss  %xmm0,0x70(%rsp)
  movss  %xmm0,0x74(%rsp)
  movl   $0xc1900000,0x78(%rsp)
  movl   $0xc1901aa0,0x7c(%rsp)
  movl   $0x0,0x80(%rsp)
  movl   $0x0,0x84(%rsp)
  movl   $0x0,0x88(%rsp)
  movl   $0x0,0x8c(%rsp)
  movl   $0x0,0x90(%rsp)
  movl   $0x0,0x94(%rsp)
  movss  D360(%rip),%xmm0
  movss  %xmm0,0xa0(%rsp)
  movss  %xmm0,0xa4(%rsp)
  movl   $0x41900000,0xa8(%rsp)
  movl   $0x418fe979,0xac(%rsp)
  movl   $0x0,0xb0(%rsp)
  movl   $0x0,0xb4(%rsp)
  movl   $0x0,0xb8(%rsp)
  movl   $0x0,0xbc(%rsp)
  movl   $0x0,0xc0(%rsp)
  movl   $0x0,0xc4(%rsp)
  movq   $0x0,0xd0(%rsp)
  movq   $0x0,0xd8(%rsp)
  movq   $0x0,0xe0(%rsp)
  movq   $0x0,0xe8(%rsp)
  movq   $0x0,0xf0(%rsp)
  test   %edx,%edx
  jle    L32
  mov    %edx,%edx
  movss  D83(%rip),%xmm2
  jmp    L85
L49:
  subss  %xmm1,%xmm0
  movaps %xmm2,%xmm1
  divss  %xmm0,%xmm1
  movss  %xmm1,0xd0(%rsp,%rax,4)
L93:
  add    $0x1,%rax
  cmp    %rdx,%rax
  je     L32
L85:
  movss  0x70(%rsp,%rax,4),%xmm1
  movss  0xa0(%rsp,%rax,4),%xmm0
  ucomiss %xmm0,%xmm1
  jp     L49
  jne    L49
  jmp    L93
L32:
  movss  0xd0(%rsp),%xmm1
  movss  (%rcx),%xmm0
  movss  D360(%rip),%xmm4
  subss  %xmm4,%xmm0
  movaps %xmm0,%xmm2
  mulss  %xmm1,%xmm2
  movss  D93(%rip),%xmm0
  movaps %xmm2,%xmm6
  addss  %xmm0,%xmm6
  movaps %xmm6,%xmm7
  movss  %xmm6,0x30(%rsp)
  mulss  (%rbx),%xmm1
  movss  0xd4(%rsp),%xmm3
  movss  0x4(%rcx),%xmm2
  subss  %xmm4,%xmm2
  mulss  %xmm3,%xmm2
  movaps %xmm2,%xmm10
  addss  %xmm0,%xmm10
  movss  %xmm10,0xc(%rsp)
  mulss  0x4(%rbx),%xmm3
  movss  %xmm3,0x24(%rsp)
  movss  0xd8(%rsp),%xmm3
  movss  0x8(%rcx),%xmm2
  subss  D139(%rip),%xmm2
  mulss  %xmm3,%xmm2
  movaps %xmm2,%xmm11
  addss  %xmm0,%xmm11
  movaps %xmm11,%xmm13
  movss  %xmm11,0x10(%rsp)
  mulss  0x8(%rbx),%xmm3
  movss  %xmm3,0x28(%rsp)
  movss  0xdc(%rsp),%xmm3
  movss  0xc(%rcx),%xmm2
  subss  D377(%rip),%xmm2
  mulss  %xmm3,%xmm2
  movaps %xmm2,%xmm15
  addss  %xmm0,%xmm15
  movss  %xmm15,0x14(%rsp)
  movaps %xmm3,%xmm2
  mulss  0xc(%rbx),%xmm2
  movss  %xmm2,0x2c(%rsp)
  movaps %xmm6,%xmm2
  mulss  %xmm6,%xmm2
  movaps %xmm10,%xmm3
  mulss  %xmm10,%xmm3
  movaps %xmm11,%xmm5
  movaps %xmm11,%xmm12
  mulss  %xmm11,%xmm5
  movaps %xmm5,%xmm11
  movss  %xmm5,0x38(%rsp)
  movaps %xmm15,%xmm14
  mulss  %xmm15,%xmm14
  movss  %xmm14,0x34(%rsp)
  mulss  D508(%rip),%xmm6
  movss  %xmm6,0x3c(%rsp)
  movaps %xmm10,%xmm9
  mulss  D123(%rip),%xmm9
  movss  %xmm9,0x40(%rsp)
  movaps %xmm7,%xmm5
  mulss  D206(%rip),%xmm5
  movss  %xmm5,0x44(%rsp)
  mulss  D246(%rip),%xmm12
  movss  %xmm12,0x48(%rsp)
  movaps %xmm13,%xmm0
  mulss  %xmm11,%xmm0
  mulss  D210(%rip),%xmm0
  movss  %xmm0,0x18(%rsp)
  movaps %xmm10,%xmm8
  mulss  D30(%rip),%xmm8
  movss  %xmm8,0x4c(%rsp)
  mulss  %xmm13,%xmm8
  movss  %xmm8,0x50(%rsp)
  movaps %xmm7,%xmm0
  mulss  %xmm2,%xmm0
  movaps %xmm0,%xmm4
  mulss  D341(%rip),%xmm4
  movss  %xmm4,0x1c(%rsp)
  movaps %xmm3,%xmm4
  mulss  D117(%rip),%xmm4
  movss  %xmm4,0x54(%rsp)
  mulss  %xmm13,%xmm4
  movss  %xmm4,0x20(%rsp)
  movaps %xmm7,%xmm0
  mulss  D207(%rip),%xmm0
  addss  D416(%rip),%xmm0
  movaps %xmm13,%xmm4
  mulss  D358(%rip),%xmm4
  addss  %xmm4,%xmm0
  mulss  %xmm10,%xmm6
  addss  %xmm6,%xmm0
  movaps %xmm9,%xmm4
  mulss  %xmm13,%xmm4
  addss  %xmm4,%xmm0
  mulss  %xmm15,%xmm5
  addss  %xmm5,%xmm0
  movaps %xmm12,%xmm4
  mulss  %xmm15,%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm10,%xmm4
  mulss  D120(%rip),%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm15,%xmm4
  mulss  D233(%rip),%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm7,%xmm4
  mulss  D471(%rip),%xmm4
  mulss  %xmm14,%xmm4
  addss  %xmm4,%xmm0
  movss  0x18(%rsp),%xmm5
  mulss  %xmm15,%xmm5
  addss  %xmm5,%xmm0
  movaps %xmm7,%xmm4
  mulss  D80(%rip),%xmm4
  mulss  %xmm3,%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm8,%xmm4
  mulss  %xmm15,%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm13,%xmm4
  mulss  D23(%rip),%xmm4
  mulss  %xmm14,%xmm4
  addss  %xmm4,%xmm0
  movss  0x1c(%rsp),%xmm5
  mulss  %xmm10,%xmm5
  addss  %xmm5,%xmm0
  movss  0x20(%rsp),%xmm4
  mulss  %xmm15,%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm10,%xmm5
  mulss  D450(%rip),%xmm5
  movss  %xmm5,0x58(%rsp)
  mulss  D278(%rip),%xmm7
  movss  %xmm7,0x5c(%rsp)
  mulss  %xmm10,%xmm7
  movss  %xmm7,0x60(%rsp)
  mulss  %xmm13,%xmm5
  movaps %xmm5,%xmm4
  mulss  %xmm14,%xmm4
  addss  %xmm4,%xmm0
  mulss  %xmm13,%xmm7
  addss  %xmm7,%xmm0
  movaps %xmm10,%xmm4
  mulss  D9(%rip),%xmm4
  mulss  %xmm11,%xmm4
  addss  %xmm4,%xmm0
  movaps %xmm2,%xmm5
  mulss  D309(%rip),%xmm5
  movss  %xmm5,0x64(%rsp)
  mulss  %xmm10,%xmm5
  movss  %xmm5,0x68(%rsp)
  movaps %xmm5,%xmm4
  mulss  %xmm13,%xmm4
  addss  %xmm0,%xmm4
  movd   %xmm4,%r14d
  pxor   %xmm4,%xmm4
  movaps %xmm1,%xmm0
  addss  %xmm4,%xmm0
  addss  0x24(%rsp),%xmm0
  addss  0x28(%rsp),%xmm0
  addss  0x2c(%rsp),%xmm0
  ucomiss %xmm4,%xmm0
  jp     L3
  movaps %xmm4,%xmm0
  jne    L3
L40:
  movss  %xmm0,(%rbx)
  mov    0xf8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L67
  movd   %r14d,%xmm0
  add    $0x108,%rsp
  pop    %rbx
  pop    %r14
  ret
L3:
  mulss  %xmm1,%xmm1
  movss  %xmm1,0x6c(%rsp)
  movss  0xc(%rsp),%xmm6
  movaps %xmm6,%xmm0
  mulss  D508(%rip),%xmm0
  subss  D335(%rip),%xmm0
  movss  0x14(%rsp),%xmm1
  mulss  D206(%rip),%xmm1
  addss  %xmm1,%xmm0
  movss  0x34(%rsp),%xmm7
  mulss  D471(%rip),%xmm7
  addss  %xmm7,%xmm0
  mulss  D80(%rip),%xmm3
  addss  %xmm3,%xmm0
  mulss  D407(%rip),%xmm2
  mulss  %xmm6,%xmm2
  addss  %xmm2,%xmm0
  movaps %xmm6,%xmm2
  mulss  D278(%rip),%xmm2
  movaps %xmm2,%xmm1
  movss  0x10(%rsp),%xmm5
  mulss  %xmm5,%xmm1
  addss  %xmm1,%xmm0
  movss  0x30(%rsp),%xmm1
  mulss  D111(%rip),%xmm1
  mulss  %xmm6,%xmm1
  mulss  %xmm5,%xmm1
  addss  %xmm1,%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  mulss  0x6c(%rsp),%xmm0
  movss  %xmm0,0x6c(%rsp)
  movss  0x24(%rsp),%xmm0
  mulss  %xmm0,%xmm0
  movss  %xmm0,0x24(%rsp)
  pxor   %xmm12,%xmm12
  addss  0x3c(%rsp),%xmm12
  movaps %xmm12,%xmm0
  movss  0x10(%rsp),%xmm5
  movaps %xmm5,%xmm2
  mulss  D123(%rip),%xmm2
  addss  %xmm2,%xmm0
  subss  D352(%rip),%xmm0
  movss  0xc(%rsp),%xmm6
  movaps %xmm6,%xmm2
  mulss  D156(%rip),%xmm2
  movaps %xmm2,%xmm1
  mulss  0x30(%rsp),%xmm1
  addss  %xmm1,%xmm0
  movaps %xmm5,%xmm2
  mulss  D30(%rip),%xmm2
  movaps %xmm2,%xmm1
  movss  0x14(%rsp),%xmm3
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  addss  0x1c(%rsp),%xmm0
  movaps %xmm6,%xmm1
  mulss  D470(%rip),%xmm1
  mulss  %xmm5,%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  movaps %xmm5,%xmm1
  mulss  D450(%rip),%xmm1
  mulss  0x34(%rsp),%xmm1
  addss  %xmm1,%xmm0
  movss  0x5c(%rsp),%xmm12
  mulss  %xmm5,%xmm12
  addss  %xmm12,%xmm0
  movss  0x38(%rsp),%xmm12
  mulss  D9(%rip),%xmm12
  addss  %xmm12,%xmm0
  movss  0x64(%rsp),%xmm14
  mulss  %xmm5,%xmm14
  addss  %xmm14,%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  mulss  0x24(%rsp),%xmm0
  addss  0x6c(%rsp),%xmm0
  movss  %xmm0,0x1c(%rsp)
  movss  0x28(%rsp),%xmm0
  mulss  %xmm0,%xmm0
  movss  %xmm0,0x24(%rsp)
  movss  0x40(%rsp),%xmm8
  addss  D358(%rip),%xmm8
  movaps %xmm8,%xmm0
  movss  0x14(%rsp),%xmm3
  movaps %xmm3,%xmm1
  mulss  D246(%rip),%xmm1
  addss  %xmm1,%xmm0
  movss  0x38(%rsp),%xmm12
  mulss  D167(%rip),%xmm12
  movaps %xmm12,%xmm1
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  movss  0x4c(%rsp),%xmm12
  mulss  %xmm3,%xmm12
  addss  %xmm12,%xmm0
  movss  0x34(%rsp),%xmm2
  movaps %xmm2,%xmm7
  mulss  D23(%rip),%xmm7
  addss  %xmm7,%xmm0
  movss  0x54(%rsp),%xmm8
  mulss  %xmm3,%xmm8
  addss  %xmm8,%xmm0
  movss  0x58(%rsp),%xmm12
  mulss  %xmm2,%xmm12
  addss  %xmm12,%xmm0
  addss  0x60(%rsp),%xmm0
  movss  0x10(%rsp),%xmm1
  mulss  D274(%rip),%xmm1
  mulss  0xc(%rsp),%xmm1
  addss  %xmm1,%xmm0
  addss  0x68(%rsp),%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  mulss  0x24(%rsp),%xmm0
  movaps %xmm0,%xmm2
  addss  0x1c(%rsp),%xmm2
  movss  %xmm2,0x1c(%rsp)
  movss  0x2c(%rsp),%xmm2
  mulss  %xmm2,%xmm2
  movss  %xmm2,0x24(%rsp)
  pxor   %xmm0,%xmm0
  addss  0x44(%rsp),%xmm0
  addss  0x48(%rsp),%xmm0
  addss  D233(%rip),%xmm0
  movss  0x14(%rsp),%xmm3
  movaps %xmm3,%xmm1
  mulss  D235(%rip),%xmm1
  mulss  0x30(%rsp),%xmm1
  addss  %xmm1,%xmm0
  addss  0x18(%rsp),%xmm0
  addss  0x50(%rsp),%xmm0
  movaps %xmm3,%xmm1
  mulss  D62(%rip),%xmm1
  movss  0x10(%rsp),%xmm5
  mulss  %xmm5,%xmm1
  addss  %xmm1,%xmm0
  addss  0x20(%rsp),%xmm0
  mulss  D42(%rip),%xmm3
  movaps %xmm3,%xmm1
  mulss  0xc(%rsp),%xmm1
  mulss  %xmm5,%xmm1
  addss  %xmm1,%xmm0
  mov    $0x2,%edi
  call   <pow@plt>
  mulss  0x24(%rsp),%xmm0
  addss  0x1c(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L40
L67:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: tag: D142
type: byte[4]
name: read_data
arguments: tag: D360
type: f32
name: read_data
arguments: tag: D83
type: f32
name: read_data
arguments: tag: D93
type: f32
name: read_data
arguments: tag: D139
type: f32
name: read_data
arguments: tag: D377
type: byte[4]
name: read_data
arguments: tag: D508
type: byte[4]
name: read_data
arguments: tag: D123
type: byte[4]
name: read_data
arguments: tag: D206
type: byte[4]
name: read_data
arguments: tag: D246
type: byte[4]
name: read_data
arguments: tag: D210
type: byte[4]
name: read_data
arguments: tag: D30
type: byte[4]
name: read_data
arguments: tag: D341
type: byte[4]
name: read_data
arguments: tag: D117
type: byte[4]
name: read_data
arguments: tag: D207
type: byte[4]
name: read_data
arguments: tag: D416
type: byte[4]
name: read_data
arguments: tag: D358
type: byte[4]
name: read_data
arguments: tag: D120
type: byte[4]
name: read_data
arguments: tag: D233
type: byte[4]
name: read_data
arguments: tag: D471
type: byte[4]
name: read_data
arguments: tag: D80
type: byte[4]
name: read_data
arguments: tag: D23
type: byte[4]
name: read_data
arguments: tag: D450
type: byte[4]
name: read_data
arguments: tag: D278
type: byte[4]
name: read_data
arguments: tag: D9
type: byte[4]
name: read_data
arguments: tag: D309
type: byte[4]
name: read_data
arguments: tag: D335
type: byte[4]
name: read_data
arguments: tag: D407
type: byte[4]
name: read_data
arguments: tag: D111
type: byte[4]
name: read_data
arguments: tag: D352
type: byte[4]
name: read_data
arguments: tag: D156
type: byte[4]
name: read_data
arguments: tag: D470
type: byte[4]
name: read_data
arguments: tag: D167
type: byte[4]
name: read_data
arguments: tag: D274
type: byte[4]
name: read_data
arguments: tag: D235
type: byte[4]
name: read_data
arguments: tag: D62
type: byte[4]
name: read_data
arguments: tag: D42
type: byte[4]

tool: {"D142": [0x00,0x00,0x70,0xc1]}

tool: D360: 15

tool: D83: 2

tool: D93: 1

tool: D139: 18

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

tool: {"D508": [0x55,0x3a,0x3b,0x3b]}

tool: {"D123": [0xd0,0xb5,0xaa,0xbb]}

tool: {"D206": [0xe1,0xbf,0x9d,0xbb]}

tool: {"D246": [0x99,0xf9,0x0e,0x3c]}

tool: {"D210": [0xa7,0xc7,0xbf,0xb8]}

tool: {"D30": [0xe8,0xe3,0x01,0xb9]}

tool: {"D341": [0xb9,0x92,0x1b,0xb8]}

tool: {"D117": [0x9a,0x51,0x90,0xb8]}

tool: {"D207": [0x0c,0xd8,0x69,0xbd]}

tool: {"D416": [0x6c,0xcb,0x97,0x3a]}

tool: {"D358": [0x62,0xe5,0x8c,0x3d]}

tool: {"D120": [0xbd,0x47,0x21,0xb9]}

tool: {"D233": [0x8a,0x0b,0x87,0x39]}

tool: {"D471": [0x87,0x55,0x72,0xb8]}

tool: {"D80": [0x2c,0x22,0xf7,0x37]}

tool: {"D23": [0xe0,0x0b,0x39,0x39]}

tool: {"D450": [0x42,0x42,0xbd,0x38]}

tool: {"D278": [0x2d,0x36,0x11,0x37]}

tool: {"D9": [0x4e,0xc5,0x56,0xb7]}

tool: {"D309": [0xb1,0x5f,0xbe,0x38]}

tool: {"D335": [0x0c,0xd8,0x69,0x3d]}

tool: {"D407": [0x16,0x5c,0xe9,0xb8]}

tool: {"D111": [0xb1,0x5f,0x3e,0x39]}

tool: {"D352": [0xbd,0x47,0x21,0x39]}

tool: {"D156": [0x2c,0x22,0x77,0x38]}

tool: {"D470": [0x9a,0x51,0x10,0xb9]}

tool: {"D167": [0xbd,0xd5,0x8f,0xb9]}

tool: {"D274": [0x4e,0xc5,0xd6,0xb7]}

tool: {"D235": [0x87,0x55,0xf2,0xb8]}

tool: {"D62": [0xe0,0x0b,0xb9,0x39]}

tool: {"D42": [0x42,0x42,0x3d,0x39]}

assistant: ```c
float target_phi(float *x, float *dx, int m) {
  int ncoeff = 20;
  float avdat = 0.1158090E-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.11092651E-07, -0.57090804E-01, 0.68796888E-01, 0.28568704E-02, -0.52096620E-02, -0.48141335E-02, 0.87265009E-02, -0.15380881E-03, 0.25757804E-03, -0.57776950E-04, -0.91447779E-04, 0.29460607E-04, -0.12387300E-03, 0.17647399E-03, -0.37091526E-04, -0.68816546E-04, 0.90245638E-04, 0.86552873E-05, -0.12801333E-04, 0.90777299E-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] * x21 * x32 + 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[18] * x32 + 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 * x21 + 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, 2));
  dx[0] = dv_target_phi;
  return v_target_phi;
}
```