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