function res = heavi( x )

som = 0;
        PP = P(12);
        l = length (PP);
        for j = 1:l
            for jj = 1:l
                if (x>abs(PP(j)-PP(jj)))
                    som = som + 3^(-2*12);
                end
            end
        end
       
res = 1.5^(2*12)*som;


end

