function [Rn,etan] = coefRn(N,delta,c) i=1; for i=1:1:N Rn(i) = sqrt(pi*delta)*exp(-delta*((i*pi/c)^2))*cos(i*pi/(2*c)); etan(i) = i*pi/c; end