A numerical iterative vortex lattice method is developed to study flow past wing(s) at high angles of attack where the separated flow is modelled using NY nascent vortex filaments. The wing itself is modelled using NX × NY bound vortex rings, where NX and NY are the number of sections along the chord and span of the wing respectively. The strength and position of the nascent vortex along the chord corresponding to the local effective angle of attack are evaluated from the residuals in viscous and potential flow, i.e. (Cl)visc - (Cl)pot and (Cm)visc - (Cm)pot. Hence, the 2D airfoil viscous Cl - α and Cm - α is required as input (from experiment, numerical analysis or CFD). Aerodynamic characteristics and section distribution along span are predicted for 3D wings at a high angle of attack. Effect of initial conditions and existence of multiple solutions in the post-stall region is studied. Results are validated with experiment. © Published under licence by IOP Publishing Ltd.