Fatiha

Hi, Please how I can improve my loops bellow:
n=184;
c=0; comp=0
D=zeros(n*(n-1)/2+n_obs*n,1);
for i=1:n
for j=i+1:n
c=c+1;
p=norm(Q(:,j)-Q(:,i));
D(c)=p-R(i)-R(j);
% possibilite de contact
if D(c)< d_risk_pp
comp=comp+1;
contact(comp,1)=i;
contact(comp,2)=j;
p=norm(Q(:,j)-Q(:,i));
E(:,comp)=(Q(:,j)-Q(:,i))/p;
end
end
end

Professional Interests: Mathematic models