user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<circlefit>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x130,%rsp
  mov    %rdi,-0x108(%rbp)
  mov    %esi,-0x10c(%rbp)
  mov    %rdx,-0x118(%rbp)
  mov    %rcx,-0x120(%rbp)
  mov    %r8,-0x128(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0xf8(%rbp)
  jmp    L61
L44:
  movsd  -0x60(%rbp),%xmm1
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x58(%rbp),%xmm1
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  0x8(%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  addl   $0x1,-0xf8(%rbp)
L61:
  mov    -0xf8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L44
  movsd  -0x60(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x58(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movl   $0x0,-0xf8(%rbp)
  jmp    L85
L29:
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x60(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  0x8(%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x20(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm2
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x20(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xf0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  -0x18(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xe8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0x10(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xe0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xd0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  movsd  -0x18(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xc8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  addl   $0x1,-0xf8(%rbp)
L85:
  mov    -0xf8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L29
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xe8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xe0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xd0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  movsd  -0xc8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movsd  -0xf0(%rbp),%xmm0
  addsd  -0xe8(%rbp),%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  movsd  -0xf0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xe8(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xe8(%rbp),%xmm1
  movsd  -0xc8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0xf0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0xe0(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0xb8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0xd0(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0xc8(%rbp),%xmm0
  mulsd  -0xd8(%rbp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0xc0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0xb8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  -0xe0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xc0(%rbp),%xmm1
  movsd  -0xb8(%rbp),%xmm2
  movsd  D21(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  -0xc0(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0xc8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0xc0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0xc0(%rbp),%xmm2
  subsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0xc0(%rbp),%xmm1
  movsd  D96(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0xc0(%rbp),%xmm0
  subsd  -0xe0(%rbp),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0xc0(%rbp),%xmm1
  movsd  D21(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  D2(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  D37(%rip),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  D33(%rip),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  D49(%rip),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movl   $0x0,-0xf4(%rbp)
L63:
  movsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x50(%rbp),%xmm2
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  addsd  -0xa0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  -0xa8(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  -0xb0(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x48(%rbp),%xmm0
  movq   D53(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm1
  movq   D53(%rip),%xmm2
  andpd  %xmm2,%xmm1
  comisd %xmm1,%xmm0
  jbe    L36
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  lea    D4(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  jmp    L42
L36:
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  addsd  -0x90(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0xa8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movsd  -0x48(%rbp),%xmm1
  divsd  -0x70(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  jbe    L8
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
L8:
  movsd  -0x50(%rbp),%xmm0
  movsd  -0x40(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x50(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movq   D53(%rip),%xmm1
  andpd  %xmm1,%xmm0
  comisd -0x80(%rbp),%xmm0
  jbe    L42
  addl   $0x1,-0xf4(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xf4(%rbp),%xmm1
  movsd  -0x78(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L63
L42:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xf4(%rbp),%xmm0
  ucomisd -0x78(%rbp),%xmm0
  jp     L66
  ucomisd -0x78(%rbp),%xmm0
  jne    L66
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  lea    D109(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
L66:
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x50(%rbp),%xmm1
  mulsd  -0xc0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xb8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0xe8(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  mulsd  -0xd0(%rbp),%xmm0
  movsd  -0xc8(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0xd8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  mov    -0x120(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  mulsd  -0xc8(%rbp),%xmm0
  movsd  -0xd0(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0xd8(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x120(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x120(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm2
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  addsd  -0xc0(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    -0x118(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x120(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x60(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x120(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x128(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0xf8(%rbp)
  jmp    L68
L104:
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x120(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xf8(%rbp),%eax
  cltq
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  0x8(%rax),%xmm0
  mov    -0x120(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm1
  movsd  -0x30(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  -0x28(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    -0x118(%rbp),%rdx
  movsd  (%rdx),%xmm1
  movq   %rax,%xmm0
  subsd  %xmm1,%xmm0
  movq   D53(%rip),%xmm1
  andpd  %xmm0,%xmm1
  mov    -0x128(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x128(%rbp),%rax
  movsd  %xmm0,(%rax)
  addl   $0x1,-0xf8(%rbp)
L68:
  mov    -0xf8(%rbp),%eax
  cmp    -0x10c(%rbp),%eax
  jl     L104
  mov    -0x128(%rbp),%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10c(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x128(%rbp),%rax
  movsd  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L111
  call   <__stack_chk_fail@plt>
L111:
  leave
  ret
```

tool_calls: name: read_data
arguments: label: D21
datatype: f32[2]
name: read_data
arguments: label: D96
datatype: f32[2]
name: read_data
arguments: label: D2
datatype: f64
name: read_data
arguments: label: D37
datatype: byte[8]
name: read_data
arguments: label: D33
datatype: byte[8]
name: read_data
arguments: label: D49
datatype: f32[2]
name: read_data
arguments: label: D53
datatype: qword
name: read_data
arguments: label: D4
datatype: string
name: read_data
arguments: label: D109
datatype: string

tool: {"D21": [0.0,2.25]}

tool: {"D96": [0.0,-2.125]}

tool: D2: 3

tool: {"D37": [0x40,0x8c,0xb5,0x78,0x1d,0xaf,0x15,0x44]}

tool: {"D33": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

tool: {"D49": [0.0,2.8125]}

tool: D53: 9223372036854775807

tool: D4: WARNING: Newton-Taubin goes in wrong direction.


tool: D109: Newton-Taubin don't converge.


assistant: ```c
void circlefit(double pts[][2], int npoints, double *radius, double *center, double *detcoeff) {
  int point, iter;
  double centroid[2];
  double Mxx, Myy, Mzz, Mxy, Mxz, Myz;
  double Mz, Cov_xy;
  double A0, A1, A2, A3, A22, A33;
  double norm_pt[3];
  double new_pt[2], old_pt[2], epsilon, maxIter, Dy;
  double Det;
  double diff[2];
  for (point = 0; point < npoints; point++) {
    centroid[0] += pts[point][0];
    centroid[1] += pts[point][1];
  }
  centroid[0] /= (double)npoints;
  centroid[1] /= (double)npoints;
  Mxx = 0.0;
  Myy = 0.0;
  Mzz = 0.0;
  Mxy = 0.0;
  Mxz = 0.0;
  Myz = 0.0;
  for (point = 0; point < npoints; point++) {
    norm_pt[0] = pts[point][0] - centroid[0];
    norm_pt[1] = pts[point][1] - centroid[1];
    norm_pt[2] = norm_pt[0] * norm_pt[0] + norm_pt[1] * norm_pt[1];
    Mxx += norm_pt[0] * norm_pt[0];
    Myy += norm_pt[1] * norm_pt[1];
    Mzz += norm_pt[2] * norm_pt[2];
    Mxy += norm_pt[0] * norm_pt[1];
    Mxz += norm_pt[0] * norm_pt[2];
    Myz += norm_pt[1] * norm_pt[2];
  }
  Mxx /= (double)npoints;
  Myy /= (double)npoints;
  Mzz /= (double)npoints;
  Mxy /= (double)npoints;
  Mxz /= (double)npoints;
  Myz /= (double)npoints;
  Mz = Mxx + Myy;
  Cov_xy = Mxx * Myy - Mxy * Mxy;
  A0 = Mxz * Mxz * Myy + Myz * Myz * Mxx - Mzz * Cov_xy - 2.0 * Mxz * Myz * Mxy + Mz * Mz * Cov_xy;
  A1 = Mzz * Mz + 4 * Cov_xy * Mz - Mxz * Mxz - Myz * Myz - Mz * Mz * Mz;
  A2 = -3 * Mz * Mz - Mzz;
  A3 = 4 * Mz;
  A22 = 2.0 * A2;
  A33 = 3.0 * A3;
  new_pt[0] = 0.0;
  new_pt[1] = 1e+20;
  epsilon = 1e-12;
  maxIter = 20;
  iter = 0;
  do {
    old_pt[1] = new_pt[1];
    new_pt[1] = A0 + new_pt[0] * (A1 + new_pt[0] * (A2 + new_pt[0] * A3));
    if (fabs(new_pt[1]) > fabs(old_pt[1])) {
      new_pt[0] = 0.0;
      printf("WARNING: Newton-Taubin goes in wrong direction.\n");
      break;
    }
    Dy = A1 + new_pt[0] * (A22 + new_pt[0] * A33);
    old_pt[0] = new_pt[0];
    new_pt[0] = old_pt[0] - new_pt[1] / Dy;
    if (new_pt[0] < 0.0) {
      new_pt[0] = 0.0;
    }
  } while ((fabs((new_pt[0] - old_pt[0]) / new_pt[0]) > epsilon) && (++iter < maxIter));
  if (iter == maxIter) {
    new_pt[0] = 0.0;
    printf("Newton-Taubin don't converge.\n");
  }
  Det = new_pt[0] * new_pt[0] - new_pt[0] * Mz + Cov_xy;
  center[0] = (Mxz * (Myy - new_pt[0]) - Myz * Mxy) / (2.0 * Det);
  center[1] = (Myz * (Mxx - new_pt[0]) - Mxz * Mxy) / (2.0 * Det);
  *radius = sqrt(((center[0] * center[0]) + (center[1] * center[1])) + Mz);
  center[0] += centroid[0];
  center[1] += centroid[1];
  *detcoeff = 0.0;
  for (point = 0; point < npoints; point++) {
    diff[0] = pts[point][0] - center[0];
    diff[1] = pts[point][1] - center[1];
    *detcoeff += fabs((sqrt(diff[0] * diff[0] + diff[1] * diff[1])) - *radius);
  }
  *detcoeff /= npoints;
}
```