1. V1n = V1 ./ repmat(sqrt(sum(V1.^2,2)),[1 3]);
2.
for iRay=1:size(Vray,2) Vray(:,iRay) = Vray(:,iRay)./norm(Vray(:,iRay));end
