% This function finds the length of a vector with a component deltaR(i).
function r = partsep(deltaR);
Rij = 0;
for i = 1:3
  Rij = Rij + deltaR(1,i)*deltaR(1,i);
end;

r = sqrt(Rij);