function sumr=functionDAPPM(g) global DAPPM_A global DAPPM_L global DAPPM_LIM global rDAPPM global rinf; A=DAPPM_A; L=DAPPM_L; L1=L+1; q=(pi)*g; Pw=sin(q)/q; sumr=rDAPPM(1)-rinf; for d=1:DAPPM_LIM sumr=sumr+(rDAPPM(d+1)-rinf)*2*cos(2.0*q*d); end sumr=Pw*Pw*sumr; end Copyright (C) 2009 Konstantin Kirillov