.

> simplify(subs({v[1]=-b/sqrt(a^2+b^2),v[2]=a/sqrt(a^2+b^2),-v[2]*x[0]+v[1]*y[0]=c/sqrt(a^2+b^2)},evalm(m)));

matrix([[-(a^2-b^2)/(a^2+b^2), -2*b*a/(a^2+b^2), 0]...