function l=compute_length(v1,v2,v3)

l(1)=norm(v1);
l(2)=norm(v2);
l(3)=norm(v3);

r=compute_rotation(v1,v2,v3);

for i=1:3
    l(i)=l(i)/sqrt(1-r(3,i)^2);
end