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
    