The thickness distribution is given by the following equation ahead of the maximum thickness:
y = A0 sqrt(x) + A1 x + A2 x**2 + A3 x**3
where ** denotes exponentiation, (t/c) is the maximum thickness to chord ratio of the airfoil, x is the position as fraction of chord, and y is the half-thickness as fraction of chord.
and by the following equation from maximum thickness to trailing edge:
y = D0 + D1(1-x) + D2 (1-x)**2 + D3 (1-x)**3
The constants A0, A1, A2, A3, D0, D1, D2, D3 are calculated from the values of maximum thickness, position of maximum thickness, and leading edge radius that are specified by the user.
The airfoil must satisfy the following constraints:
D1 is the negative of the trailing edge slope.
These conditions are sufficient to determine all of the A and D terms in the polynomial equations.