<<VilGebra`

N1={1,2,-1};
N2={2,-1,3};
d=Cross[N1,N2]
P={1,y,z}/.Solve[{x+2y-z==3,2x-y+3z==1}][[1]]/.x->1
EcuacVecToSime[P+t d]