function dist_double_power_wts_plot(k1,k2,b);

%DIST_POWER_WTS computes power weights of the form
% 
% a(dij) = [1 - (dij/b)^k1]^k2 , dij <= b
%        = 0  ,  dij > b
% 
% for any integer k > 1, where b > Dmin = bandwidth. These weights are then row normalized as
%     
% wij = s(dij)/sum[a(dih): h = 1,..n]%


del = b/100; %increment

d = [0:del:b];

N = length(d);

a = zeros(N,1);

for i = 1:N
    
    a(i) = [1 - (d(i)./b).^k1].^k2;   
     
end     
 
r = [0:del:b]; %plot range

M = length(r);

A = zeros(1,M);

A(1:N) = a;

plot(r,A,'k','LineWidth',3);

xlabel('Distance','Fontsize',12, 'FontWeight', 'bold');

ylabel('Influence Weight','Fontsize',12, 'FontWeight', 'bold');

title(['DOUBLE-POWER KERNEL (k1 = ',num2str(k1),', k2 = ',num2str(k2),', b = ',num2str(b),')'],'Fontsize',12, 'FontWeight', 'bold');