function [Q1n,Q2n,etan] = coefQn12(N,delta,c)

    i=1;
    etan = raizFun1(N,h,c,H1);
    
    for i=1:1:N
        Q1n(i) = quadl(@myfun1,0,c,1.e-9,0,c,etan(i),delta);
        Q2n(i) = quadl(@myfun2,0,c,1.e-9,0,c,etan(i),delta);     
    end
    