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