[BEMTHover]
  free  A
  free  Delta_r[:]
  free  Omega
  free  P
  free  R
  free  V_i[:]
  free  dC_P[:]
  free  dC_T[:]
  free  dP[:]
  free  r[:]
  fixed Omega_max
  fixed Omega_min
  fixed R_max
  fixed R_min
  fixed dr_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[0] ≥ dr_min
  cons  Delta_r[1] ≥ dr_min
  cons  Delta_r[2] ≥ dr_min
  cons  Delta_r[3] ≥ dr_min
  cons  Delta_r[4] ≥ 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[0] = rho·A·(Omega·R)²·dC_T[0]
  cons  xi[1] = rho·A·(Omega·R)²·dC_T[1]
  cons  xi[2] = rho·A·(Omega·R)²·dC_T[2]
  cons  xi[3] = rho·A·(Omega·R)²·dC_T[3]
  cons  xi[4] = rho·A·(Omega·R)²·dC_T[4]
  cons  V_i[0]²·r[0]·Delta_r[0]/(dC_T[0]·(Omega·R)²) = 0.25
  cons  V_i[1]²·r[1]·Delta_r[1]/(dC_T[1]·(Omega·R)²) = 0.25
  cons  V_i[2]²·r[2]·Delta_r[2]/(dC_T[2]·(Omega·R)²) = 0.25
  cons  V_i[3]²·r[3]·Delta_r[3]/(dC_T[3]·(Omega·R)²) = 0.25
  cons  V_i[4]²·r[4]·Delta_r[4]/(dC_T[4]·(Omega·R)²) = 0.25
  cons  V_i[0]³·r[0]·Delta_r[0]/(dC_P[0]·(Omega·R)³) = 0.25
  cons  V_i[1]³·r[1]·Delta_r[1]/(dC_P[1]·(Omega·R)³) = 0.25
  cons  V_i[2]³·r[2]·Delta_r[2]/(dC_P[2]·(Omega·R)³) = 0.25
  cons  V_i[3]³·r[3]·Delta_r[3]/(dC_P[3]·(Omega·R)³) = 0.25
  cons  V_i[4]³·r[4]·Delta_r[4]/(dC_P[4]·(Omega·R)³) = 0.25
  cons  dP[0] = rho·A·(Omega·R)³·dC_P[0]
  cons  dP[1] = rho·A·(Omega·R)³·dC_P[1]
  cons  dP[2] = rho·A·(Omega·R)³·dC_P[2]
  cons  dP[3] = rho·A·(Omega·R)³·dC_P[3]
  cons  dP[4] = rho·A·(Omega·R)³·dC_P[4]
  cons  P ≥ dP[:].sum()
