user: 这段汇编代码用C语言表示是什么样的?

```asm
<GPS_Math_Cassini_EN_To_LatLon>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x1f0,%rsp
  movsd  %xmm0,-0x1a8(%rbp)
  movsd  %xmm1,-0x1b0(%rbp)
  mov    %rdi,-0x1b8(%rbp)
  mov    %rsi,-0x1c0(%rbp)
  movsd  %xmm2,-0x1c8(%rbp)
  movsd  %xmm3,-0x1d0(%rbp)
  movsd  %xmm4,-0x1d8(%rbp)
  movsd  %xmm5,-0x1e0(%rbp)
  movsd  %xmm6,-0x1e8(%rbp)
  movsd  %xmm7,-0x1f0(%rbp)
  mov    -0x1d0(%rbp),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1d0(%rbp)
  mov    -0x1c8(%rbp),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1c8(%rbp)
  movsd  D321(%rip),%xmm0
  movsd  %xmm0,-0x198(%rbp)
  movsd  D463(%rip),%xmm0
  movsd  %xmm0,-0x190(%rbp)
  movsd  -0x1e8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x188(%rbp)
  movsd  -0x1f0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x180(%rbp)
  movsd  -0x188(%rbp),%xmm0
  subsd  -0x180(%rbp),%xmm0
  divsd  -0x188(%rbp),%xmm0
  movsd  %xmm0,-0x178(%rbp)
  movsd  -0x178(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x170(%rbp)
  movsd  -0x178(%rbp),%xmm0
  mulsd  -0x170(%rbp),%xmm0
  movsd  %xmm0,-0x168(%rbp)
  movsd  -0x170(%rbp),%xmm1
  movsd  D168(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x160(%rbp)
  movsd  -0x168(%rbp),%xmm1
  movsd  D564(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D7(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x158(%rbp)
  movsd  -0x178(%rbp),%xmm0
  movsd  D30(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D301(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  -0x160(%rbp),%xmm1
  movsd  D374(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x168(%rbp),%xmm2
  movsd  D116(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D293(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x150(%rbp)
  movsd  -0x178(%rbp),%xmm1
  movsd  D168(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D221(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  -0x160(%rbp),%xmm0
  movsd  D443(%rip),%xmm2
  divsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x158(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x148(%rbp)
  movsd  -0x170(%rbp),%xmm1
  movsd  D447(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D293(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x158(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x140(%rbp)
  movsd  -0x168(%rbp),%xmm1
  movsd  D82(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D292(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x138(%rbp)
  movsd  -0x150(%rbp),%xmm0
  mulsd  -0x1c8(%rbp),%xmm0
  movsd  %xmm0,-0x130(%rbp)
  movsd  -0x1c8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x148(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x128(%rbp)
  movsd  -0x1c8(%rbp),%xmm1
  movsd  D30(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x140(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x120(%rbp)
  movsd  -0x1c8(%rbp),%xmm1
  movsd  D251(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x138(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x118(%rbp)
  movsd  -0x130(%rbp),%xmm0
  subsd  -0x128(%rbp),%xmm0
  addsd  -0x120(%rbp),%xmm0
  subsd  -0x118(%rbp),%xmm0
  movsd  -0x1e8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x110(%rbp)
  movsd  D301(%rip),%xmm0
  subsd  -0x178(%rbp),%xmm0
  movsd  %xmm0,-0x108(%rbp)
  movsd  D104(%rip),%xmm0
  mov    -0x108(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x100(%rbp)
  movsd  D301(%rip),%xmm0
  subsd  -0x100(%rbp),%xmm0
  movsd  -0x100(%rbp),%xmm2
  movsd  D301(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf8(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  -0xf0(%rbp),%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  -0xe8(%rbp),%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  movsd  -0xf8(%rbp),%xmm1
  movsd  D168(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D34(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0xe8(%rbp),%xmm2
  movsd  D544(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D443(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  movsd  -0xf0(%rbp),%xmm1
  movsd  D6(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D177(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0xe0(%rbp),%xmm2
  movsd  D499(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D443(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  movsd  -0xe8(%rbp),%xmm1
  movsd  D152(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D445(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movsd  -0xe0(%rbp),%xmm1
  movsd  D366(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D201(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  movsd  D415(%rip),%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  -0x1a8(%rbp),%xmm0
  subsd  -0x1d8(%rbp),%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  -0x1b0(%rbp),%xmm0
  subsd  -0x1e0(%rbp),%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0x110(%rbp),%xmm0
  addsd  -0xa8(%rbp),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0x1e8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x150(%rbp),%xmm1
  movsd  -0xa0(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x98(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xd8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  D30(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xd0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  D251(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xc8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  D221(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xc0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x98(%rbp),%xmm0
  addsd  -0x90(%rbp),%xmm0
  addsd  -0x88(%rbp),%xmm0
  addsd  -0x80(%rbp),%xmm0
  movsd  -0x78(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x190(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0xb8(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L75
  movsd  -0x190(%rbp),%xmm0
  addsd  -0xb8(%rbp),%xmm0
  comisd -0x70(%rbp),%xmm0
  jbe    L75
  mov    -0x1b8(%rbp),%rax
  movsd  -0x190(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c0(%rbp),%rax
  movsd  -0x1d0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L34
L75:
  movsd  -0x190(%rbp),%xmm0
  movq   D149(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  subsd  -0xb8(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L134
  movsd  -0xb8(%rbp),%xmm0
  subsd  -0x190(%rbp),%xmm0
  comisd -0x70(%rbp),%xmm0
  jbe    L134
  movsd  -0x190(%rbp),%xmm0
  movq   D149(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x1b8(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x1c0(%rbp),%rax
  movsd  -0x1d0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L34
L134:
  mov    -0x70(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x68(%rbp)
  mov    -0x70(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x60(%rbp)
  mov    -0x70(%rbp),%rax
  movq   %rax,%xmm0
  call   <tan@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x178(%rbp),%xmm0
  mulsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x68(%rbp),%xmm1
  movsd  D301(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movsd  D104(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x48(%rbp)
  movsd  -0x1e8(%rbp),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x108(%rbp),%xmm1
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0xb0(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm0
  mulsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x30(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x30(%rbp),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D168(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D301(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x40(%rbp),%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm2
  divsd  -0x38(%rbp),%xmm2
  movsd  -0x28(%rbp),%xmm0
  movsd  D34(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x8(%rbp),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  movsd  D214(%rip),%xmm3
  divsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm2
  movapd %xmm2,%xmm1
  movsd  -0x70(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  mov    -0x1b8(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  movsd  D168(%rip),%xmm1
  movapd %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movsd  -0x30(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x8(%rbp),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  D447(%rip),%xmm2
  divsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  divsd  -0x60(%rbp),%xmm0
  addsd  -0x1d0(%rbp),%xmm0
  mov    -0x1c0(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x1b8(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd -0x190(%rbp),%xmm0
  jbe    L116
  mov    -0x1b8(%rbp),%rax
  movsd  -0x190(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L43
L116:
  mov    -0x1b8(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x190(%rbp),%xmm0
  movq   D149(%rip),%xmm2
  xorpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jbe    L43
  movsd  -0x190(%rbp),%xmm0
  movq   D149(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x1b8(%rbp),%rax
  movsd  %xmm0,(%rax)
L43:
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D290(%rip),%xmm0
  jbe    L64
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm0
  subsd  -0x198(%rbp),%xmm0
  mov    -0x1c0(%rbp),%rax
  movsd  %xmm0,(%rax)
L64:
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D105(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L133
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  -0x198(%rbp),%xmm0
  mov    -0x1c0(%rbp),%rax
  movsd  %xmm0,(%rax)
L133:
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D290(%rip),%xmm0
  jbe    L93
  mov    -0x1c0(%rbp),%rax
  movsd  D290(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L34
L93:
  mov    -0x1c0(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D105(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L34
  mov    -0x1c0(%rbp),%rax
  movsd  D105(%rip),%xmm0
  movsd  %xmm0,(%rax)
L34:
  mov    -0x1c0(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movq   %xmm0,%rax
  mov    -0x1c0(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x1b8(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movq   %xmm0,%rax
  mov    -0x1b8(%rbp),%rdx
  mov    %rax,(%rdx)
  nop
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_tag: D321
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D463
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D168
type_guessed: f64
name: disassemble_data
arguments: data_tag: D564
type_guessed: f64
name: disassemble_data
arguments: data_tag: D7
type_guessed: f64
name: disassemble_data
arguments: data_tag: D30
type_guessed: f64
name: disassemble_data
arguments: data_tag: D301
type_guessed: f64
name: disassemble_data
arguments: data_tag: D374
type_guessed: f64
name: disassemble_data
arguments: data_tag: D116
type_guessed: f64
name: disassemble_data
arguments: data_tag: D293
type_guessed: f64
name: disassemble_data
arguments: data_tag: D221
type_guessed: f64
name: disassemble_data
arguments: data_tag: D443
type_guessed: f64
name: disassemble_data
arguments: data_tag: D447
type_guessed: f64
name: disassemble_data
arguments: data_tag: D82
type_guessed: f64
name: disassemble_data
arguments: data_tag: D292
type_guessed: f64
name: disassemble_data
arguments: data_tag: D251
type_guessed: f64
name: disassemble_data
arguments: data_tag: D104
type_guessed: f64
name: disassemble_data
arguments: data_tag: D34
type_guessed: f64
name: disassemble_data
arguments: data_tag: D544
type_guessed: f64
name: disassemble_data
arguments: data_tag: D6
type_guessed: f64
name: disassemble_data
arguments: data_tag: D177
type_guessed: f64
name: disassemble_data
arguments: data_tag: D499
type_guessed: f64
name: disassemble_data
arguments: data_tag: D152
type_guessed: f64
name: disassemble_data
arguments: data_tag: D445
type_guessed: f64
name: disassemble_data
arguments: data_tag: D366
type_guessed: f64
name: disassemble_data
arguments: data_tag: D201
type_guessed: f64
name: disassemble_data
arguments: data_tag: D415
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D149
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D214
type_guessed: f64
name: disassemble_data
arguments: data_tag: D290
type_guessed: f64[2]
name: disassemble_data
arguments: data_tag: D105
type_guessed: byte[8]

tool: {"D321": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: {"D463": [0x1e,0x26,0x44,0x54,0xfb,0x21,0xf9,0x3f]}

tool: D168: 3

tool: D564: 45

tool: D7: 1024

tool: D30: 4

tool: D301: 1

tool: D374: 64

tool: D116: 5

tool: D293: 256

tool: D221: 8

tool: D443: 32

tool: D447: 15

tool: D82: 35

tool: D292: 3072

tool: D251: 6

tool: D104: 0.5

tool: D34: 2

tool: D544: 27

tool: D6: 21

tool: D177: 16

tool: D499: 55

tool: D152: 151

tool: D445: 96

tool: D366: 1097

tool: D201: 512

tool: {"D415": [0xf1,0x68,0xe3,0x88,0xb5,0xf8,0xe4,0x3e]}

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

tool: D214: 24

tool: {"D290": [3.141592653589,-3.141592653589]}

tool: {"D105": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x09,0xc0]}

assistant: ```c
void GPS_Math_Cassini_EN_To_LatLon(double E, double N, double *phi, double *lambda, double phi0, double M0, double E0, double N0, double a, double b) {
  double p2;
  double po2;
  double a2;
  double b2;
  double e2;
  double e4;
  double e6;
  double AM0;
  double c0;
  double c1;
  double c2;
  double c3;
  double om0;
  double A0;
  double A1;
  double A2;
  double A3;
  double j;
  double te4;
  double phi0s2;
  double phi0s4;
  double phi0s6;
  double lat;
  double x;
  double E1;
  double E2;
  double E3;
  double E4;
  double dx;
  double dy;
  double mu;
  double mus2;
  double mus4;
  double mus6;
  double mus8;
  double M1;
  double phi1;
  double phi1s;
  double phi1c;
  double phi1t;
  double T;
  double T1;
  double N1;
  double R1;
  double RD;
  double DD;
  double D2;
  double D3;
  double D4;
  double D5;
  double tol;
  M0 = GPS_Math_Deg_To_Rad(M0);
  phi0 = GPS_Math_Deg_To_Rad(phi0);
  p2 = (double)3.141592653589 * (double)2.;
  po2 = (double)3.141592653589 / (double)2.;
  a2 = a * a;
  b2 = b * b;
  e2 = (a2 - b2) / a2;
  e4 = e2 * e2;
  e6 = e2 * e4;
  te4 = (double)3.0 * e4;
  j = (double)45. * e6 / (double)1024.;
  c0 = (double)1.0 - e2 / (double)4. - te4 / (double)64. - (double)5. * e6 / (double)256.;
  c1 = (double)3. * e2 / (double)8. + te4 / (double)32. + j;
  c2 = (double)15. * e4 / (double)256. + j;
  c3 = (double)35. * e6 / (double)3072.;
  lat = c0 * phi0;
  phi0s2 = c1 * sin((double)2. * phi0);
  phi0s4 = c2 * sin((double)4. * phi0);
  phi0s6 = c3 * sin((double)6. * phi0);
  AM0 = a * (lat - phi0s2 + phi0s4 - phi0s6);
  om0 = (double)1.0 - e2;
  x = pow(om0, (double)0.5);
  E1 = ((double)1.0 - x) / ((double)1.0 + x);
  E2 = E1 * E1;
  E3 = E1 * E2;
  E4 = E1 * E3;
  A0 = (double)3. * E1 / (double)2. - (double)27. * E3 / (double)32.;
  A1 = (double)21. * E2 / (double)16. - (double)55. * E4 / (double)32.;
  A2 = (double)151. * E3 / (double)96.;
  A3 = (double)1097. * E4 / (double)512.;
  tol = (double)1.e-5;
  dx = E - E0;
  dy = N - N0;
  M1 = AM0 + dy;
  mu = M1 / (a * c0);
  mus2 = A0 * sin((double)2. * mu);
  mus4 = A1 * sin((double)4. * mu);
  mus6 = A2 * sin((double)6. * mu);
  mus8 = A3 * sin((double)8. * mu);
  phi1 = mu + mus2 + mus4 + mus6 + mus8;
  if ((((po2 - tol) < phi1) && (phi1 < (po2 + tol)))) {
    *phi = po2;
    *lambda = M0;
  } else if ((((-po2 - tol) < phi1) && (phi1 < (-po2 + tol)))) {
    *phi = -po2;
    *lambda = M0;
  } else {
    phi1s = sin(phi1);
    phi1c = cos(phi1);
    phi1t = tan(phi1);
    T1 = phi1t * phi1t;
    RD = pow((double)1. - e2 * phi1s * phi1s, (double).5);
    N1 = a / RD;
    R1 = N1 * om0 / (RD * RD);
    DD = dx / N1;
    D2 = DD * DD;
    D3 = DD * D2;
    D4 = DD * D3;
    D5 = DD * D4;
    T = (double)1. + (double)3. * T1;
    *phi = phi1 - (N1 * phi1t / R1) * (D2 / (double)2. - T * D4 / (double)24.);
    *lambda = M0 + (DD - T1 * D3 / (double)3. + T * T1 * D5 / (double)15.) / phi1c;
    if (*phi > po2)
      *phi = po2;
    else if (*phi < -po2)
      *phi = -po2;
    if (*lambda > 3.141592653589)
      *lambda -= p2;
    if (*lambda < -3.141592653589)
      *lambda += p2;
    if (*lambda > 3.141592653589)
      *lambda = 3.141592653589;
    else if (*lambda < -3.141592653589)
      *lambda = -3.141592653589;
  }
  *lambda = GPS_Math_Rad_To_Deg(*lambda);
  *phi = GPS_Math_Rad_To_Deg(*phi);
  return;
}
```