% % H = makeH(coord, [L, a])%% Parameter:% coord  - coordinate of nodes coord(N,2)% L      - a routing matrix L(N,2) where for each link <i>, L(i,1) contains %          the index of the source node and L(i,2) is the index of the destination node. %          Default L is each node talks to itself (hence H(i,j) is fading between nodes i and j).% a      - fading coefficient (a < -2)%% Returns:% H(i,j) - fading between source of link i and destination of link j.%function H = makeH(coord,L,a)if nargin < 2	L = [];endif nargin < 3   a = -4;end    N = size(coord,1);if isempty(L)	L = [(1:N)', (1:N)'];endl = size(L,1);H = zeros(l);for ii=1:l	for jj=1:l		i = L(ii,1);		j = L(jj,2);     	H(ii,jj) = sqrt((coord(i,1)-coord(j,1))^2 + (coord(i,2)-coord(j,2))^2);     	if H(ii,jj) > 0 	     	H(ii,jj) = H(ii,jj)^a;	   else	   	H(ii,jj) = 1e10;	   end	endend