a:=vector(3): b:=vector(3): u:= vector(3): v:=vector(3): eq1 := evalm((a+s*u) - (b+t*v)); sol1 := solve({eq1[1],eq1[2]},{s,t}); tt := op(2,sol1[1]); ss := op(2,sol1[2]); interlines := proc(a,u,b,v) eq := evalm(a+s*u-b-t*v); end: interlines([1,0,1],[-2,0,3],[1,0,2],[1,0,-1]); st := interlines([1,0,1],[-2,0,3],[1,0,2],[1,0,-1]): st[1]; st[2];