Program Foil CFD :

Solving The Entire Steady Incompressible Subsonic Flow Field About Two Dimensional Lifting And Non-Lifting Surfaces , By Integrating The Ordinary Differential Elliptic Laplace Equation Describing The Velocity Function , Over A Delaunay Triangular Grid Using The Numerical Dormand-Prince Fifth Order Finite Difference Method With Neumann Boundary Conditions