function [a, der, int] = coeficientes() 
    N = 51;
    x = [1:N];
    
    a = 2.^(-x);
    der = a(2:N).*x(1:N-1);
    int = [0, a./[1:N]];
end