/**
   * @param wingAr
   * @param verticalAr
   * @param surfaceWing
   * @param surfaceVertical
   * @param sweepC2vertical
   * @param clAlphaVertical
   * @param mach
   * @param kFv
   * @param kWv
   * @param kHv
   * @return
   */
  public static double calcCNbetaVerticalTail(
      double wingAr,
      double verticalAr,
      double armVertical,
      double wingSpan,
      double surfaceWing,
      double surfaceVertical,
      double sweepC2vertical,
      double clAlphaVertical,
      double mach,
      double kFv,
      double kWv,
      double kHv) {

    return calcCNbetaVerticalTail(
        LiftCalc.calculateCLalphaHelmboldDiederich(
            verticalAr, clAlphaVertical, sweepC2vertical, mach),
        kFv,
        kWv,
        kHv,
        armVertical,
        surfaceVertical,
        surfaceWing,
        wingSpan);
  }
  /**
   * @author Lorenzo Attanasio
   * @param mach
   * @param ar
   * @param semispan
   * @param sweepHalfChordEq
   * @param acToCgDistance
   * @param wingMAC
   * @param alpha0L
   * @param cLalpha
   * @param yStations
   * @param clAlphaVsY
   * @param chordsVsY
   * @return
   */
  public static double calcCM0LS(
      double mach,
      double ar,
      double semispan,
      double sweepHalfChordEq,
      double acToCgDistance,
      double wingMAC,
      double alpha0L,
      double cLalpha,
      double[] yStations,
      double[] clAlphaVsY,
      double[] chordsVsY) {

    return cLalpha
        + LiftCalc.calculateLiftCoefficientAtAlpha0(alpha0L, cLalpha) * (acToCgDistance / wingMAC);
  }