user: Please show me the C implementation that matches the assembly code listed below:

```asm
<gauss_integral>:
  endbr64
  comisd D71(%rip),%xmm0
  ja     L5
  movapd %xmm0,%xmm1
  movsd  D135(%rip),%xmm0
  pxor   %xmm2,%xmm2
  comisd %xmm1,%xmm0
  jbe    L19
  movapd %xmm2,%xmm0
  ret
  nopl   (%rax)
L19:
  sub    $0x18,%rsp
  movapd %xmm1,%xmm0
  andpd  D98(%rip),%xmm0
  movsd  D67(%rip),%xmm3
  mulsd  D24(%rip),%xmm0
  movapd %xmm1,%xmm4
  movsd  %xmm1,(%rsp)
  mulsd  D79(%rip),%xmm4
  movapd %xmm3,%xmm2
  addsd  %xmm3,%xmm0
  divsd  %xmm0,%xmm2
  movapd %xmm1,%xmm0
  xorpd  D80(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  movsd  %xmm2,0x8(%rsp)
  call   <exp@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  (%rsp),%xmm1
  movapd %xmm0,%xmm4
  movsd  D70(%rip),%xmm0
  mov    D67(%rip),%rax
  mulsd  D16(%rip),%xmm4
  mulsd  %xmm2,%xmm0
  movq   %rax,%xmm3
  subsd  D42(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D68(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D142(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D17(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  movq   %rax,%xmm2
  mulsd  %xmm4,%xmm0
  subsd  %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L6
  movapd %xmm2,%xmm0
  add    $0x18,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L5:
  movsd  D67(%rip),%xmm2
  movapd %xmm2,%xmm0
  ret
  nopl   (%rax)
L6:
  subsd  %xmm2,%xmm3
  add    $0x18,%rsp
  movapd %xmm3,%xmm2
  movapd %xmm2,%xmm0
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D71
datatype: f64[2]
name: disassemble_data
arguments: data_label: D135
datatype: f32[2]
name: disassemble_data
arguments: data_label: D98
datatype: qword[2]
name: disassemble_data
arguments: data_label: D67
datatype: f64
name: disassemble_data
arguments: data_label: D24
datatype: f64
name: disassemble_data
arguments: data_label: D79
datatype: f32[2]
name: disassemble_data
arguments: data_label: D80
datatype: f32[4]
name: disassemble_data
arguments: data_label: D70
datatype: f64
name: disassemble_data
arguments: data_label: D16
datatype: f64
name: disassemble_data
arguments: data_label: D42
datatype: f64
name: disassemble_data
arguments: data_label: D68
datatype: f64
name: disassemble_data
arguments: data_label: D142
datatype: f64
name: disassemble_data
arguments: data_label: D17
datatype: f64

tool: {"D71": [6.0,-6.0]}

tool: {"D135": [0.0,-2.375]}

tool: {"D98": [0x7fffffffffffffff,0x0000000000000000]}

tool: D67: 1

tool: D24: 0.2316419

tool: {"D79": [0.0,1.75]}

tool: {"D80": [0.0,-0.0,0.0,0.0]}

tool: D70: 1.330274429

tool: D16: 0.3989423

tool: D42: 1.821255978

tool: D68: 1.781477937

tool: D142: 0.356563782

tool: D17: 0.31938153

assistant: ```c
double gauss_integral(double z) {
  double b1;
  double b2;
  double b3;
  double b4;
  double b5;
  double p;
  double c2;
  double a;
  double t;
  double b;
  double n;
  if (z > 6.0)
    return 1.0;
  if (z < -6.0)
    return 0.0;
  b1 = 0.31938153;
  b2 = -0.356563782;
  b3 = 1.781477937;
  b4 = -1.821255978;
  b5 = 1.330274429;
  p = 0.2316419;
  c2 = 0.3989423;
  a = fabs(z);
  t = 1.0 / (1.0 + a * p);
  b = c2 * exp((-z) * (z / 2.0));
  n = ((((b5 * t + b4) * t + b3) * t + b2) * t + b1) * t;
  n = 1.0 - b * n;
  if (z < 0.0)
    n = 1.0 - n;
  return n;
}
```