function [Z,f] = modlineal(N,H,c)
%function [Z,f] = modlineal(N,H,c)
B=0;k=1;
    
    for Z=0:0.1:c
      
            for n=1:1:N
                etan = n*pi/c;
                B = B + (1-(-1)^n)*sin(etan*H)*sin(etan*Z)/(etan^2);
            end
        
            f(k) = 2*B/(H*c);
            B = 0;
            k = k + 1;
    end
    Z=0:0.1:c;
        