dis := (n,D1,D2) -> abs(D2-D1)/len(n): len := vec -> sqrt( innerprod(vec,vec) ): ans := dis([1,1,2],2,4); evalf(%);