[BEMTHover]
  free  A
  free  dC_P[:]
  free  dC_T[:]
  free  Delta_r[:]
  free  dP[:]
  free  Omega
  free  P
  free  R
  free  r[:]
  free  V_i[:]
  fixed dr_min
  fixed Omega_max
  fixed Omega_min
  fixed R_max
  fixed R_min
  fixed rho
  fixed xi[:]
  cons  A = 3.14⋅R²
  cons  R ≤ R_max
  cons  R ≥ R_min
  cons  Omega ≤ Omega_max
  cons  Omega ≥ Omega_min
  cons  Delta_r[:] ≥ dr_min
  cons  r[0] = Delta_r[0]/2
  cons  r[1] ≥ Delta_r[:1].sum()
  cons  r[2] ≥ Delta_r[:2].sum()
  cons  r[3] ≥ Delta_r[:3].sum()
  cons  r[4] ≥ Delta_r[:4].sum()
  cons  r[1] ≥ r[0] + 0.5⋅Delta_r[0] + 0.5⋅Delta_r[1]
  cons  r[2] ≥ r[1] + 0.5⋅Delta_r[1] + 0.5⋅Delta_r[2]
  cons  r[3] ≥ r[2] + 0.5⋅Delta_r[2] + 0.5⋅Delta_r[3]
  cons  r[4] ≥ r[3] + 0.5⋅Delta_r[3] + 0.5⋅Delta_r[4]
  cons  r[4] ≤ 1
  cons  Delta_r[:].sum() ≤ 1
  cons  xi[:] = rho⋅A⋅(Omega⋅R)²⋅dC_T[:]
  cons  V_i[:]²⋅r[:]⋅Delta_r[:]/(dC_T[:]⋅(Omega⋅R)²) = 0.25
  cons  V_i[:]³⋅r[:]⋅Delta_r[:]/(dC_P[:]⋅(Omega⋅R)³) = 0.25
  cons  dP[:] = rho⋅A⋅(Omega⋅R)³⋅dC_P[:]
  cons  P ≥ dP[:].sum()
