A numerical method for second order thin airfoil theory