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 integrator , with neumann boundary conditions