function  distversusikegami

d=linspace(0.02,0.5,200);

% distance dependence interference

intf1=1./(2.*((1.5-d).^(-4))+((1.5-(d./2)).^(-4))+((1.5+(d./2)).^(-4))+((1.5+d).^(-4))+(1.5.^(-4)));
intf1=intf1.*d.^(-4);
intf1=10.*log10(intf1)

% cost231 for interference
intf2=10.^(cost231wi(1.5+d)./10)+10.^(cost231wi(1.5+d./2)./10)+10.^(cost231wi(1.5)./10)+10.^(cost231wi(1.5-d)./10)+10.^(cost231wi(1.5-d)./10)+10.^(cost231wi(1.5-d./2)./10);
intf2=10.*log10(intf2);
power=cost231wi(d);
intf2=power-intf2
% end

plot (d,intf1,'k-. ');
hold on;
plot (d,intf2,'k');
title('CIR');
legend('distance-dependence','cost231');
ylabel('CIR in dB');
xlabel('distance in km');