function [Qn,etan] = coefQn(N,delta,c)

    i=1;
    
    for i=1:1:N
        Qn(i) = quadl(@myfun,0,c,1.e-9,0,c,i,delta);
        etan(i) = i*pi/c;
    end
    