Just in case anyone is interested the coefficient of drag varies with aspect ratio as follows: Cd'= Cd + Cl^2/pi( 1/A' - 1/A)
For computation of span-wise load distribution for non-elliptical loading Munk's integrals can be used with the help of computers...