 ;;(integral-abstracted (hash-map :corrected-max-camber 0 :corrected-position-camber 0 :corrected-thickness 0.45) upper-surface-y-function left-rule 100)
 
 ;;(upper-surface-y-function (hash-map :corrected-max-camber 0 :corrected-position-camber 0.5 :corrected-thickness 0.45) 0)
+
+(* (integral-production (NACA-design 2 4 12) gradient-forward-polar-function 0 1.369 (/ 1.369 200) 200) 1.369)
+
+(integral-abstracted (NACA-design 2 4 12) gradient-forward-polar-function middle-rule 0 (theta-switch (NACA-design 2 4 12)) 200)
+
+(trapazoid-integral (NACA-design 2 4 12) gradient-forward-polar-function 0 (theta-switch (NACA-design 2 4 12)) 10000)
+
+(simpson-integral (NACA-design 2 4 12) gradient-forward-polar-function 0 (theta-switch (NACA-design 2 4 12)) 1000)
+
+(integral-abstracted (NACA-design 2 4 12) gradient-aft-polar-function left-rule (theta-switch (NACA-design 2 4 12)) (Math/PI) 200)
+
+(simpson-integral (NACA-design 2 4 12) gradient-aft-polar-function (theta-switch (NACA-design 2 4 12)) (Math/PI) 200)
+
+
+

  The integral multiplies by the derivative of the x function because the functions are
  paremetric."
  [individual]
  (- (simpson-integral individual adjusted-upper-function 90)
     (simpson-integral individual adjusted-lower-function 90)))

;;;---------------- Coefficantof lift ----------------
(defn symetric-coefficient-of-lift

        theta-switch-constant (theta-switch individual)]
    (- attack-angle
       (* (/ 1 Math/PI)
          (+ (integral-abstracted individual gradient-forward-polar-function middle-rule 0 theta-switch-constant 600)
             (integral-abstracted individual gradient-aft-polar-function left-rule theta-switch-constant Math/PI 600))))))

(defn A-1
  [individual run-constants]
  (let [theta-switch-constant (theta-switch individual)]
    (* (/ 2 Math/PI)
       (+ (integral-abstracted individual GFPCos middle-rule 0 theta-switch-constant 850)
          (integral-abstracted individual GAPCos left-rule theta-switch-constant Math/PI 850)))))


(defn cambered-coefficient-of-lift
