Hi. I have written my own 3d panel method code in MATLAB. it could calculate correct pressure coefficients of a NACA0012 airfoil at zero AoA, but fails to get correct result at 5 degree of AOA, pressure coefficients seem to diverge near trailing point. I think this may be caused by wrong calculation of source/doublet influence coefficients (problem with arc tangents functions).

I added figures of both results. could anyone help me with this problem?

Similar questions and discussions