P := [5,7]: Q := [6,8]: PQ := evalm(Q-P); len := v -> sqrt(dotprod(v,v)): lenth_of_PQ := len(PQ);