I'm working on a basic 2D-"panelmethod" program. Currently it's "no-lift" but I would like to expand it to lifting bodies. The routines I'm using are loaned from my litterature and the program is working just fine for non-lifting bodies. It uses sources distributed over the panels to calculate influence coefficients Tij,Nij (induced velocities over panel). In the no-lift case I just need to sum these coefficients over all panels and I'm done. The problems arise when expanding to lifting-bodies and one term of vorticity is introduced along with few more Nij and Tij. The question is how i should implement the vorticity and extra influence coefficients to obtain the velocity at a certain panel?