>    with (linalg):

>    A := vector( [u/(u^2+v^2+w^2),v/(u^2+v^2+w^2),w/(u^2+v^2+w^2)] );

A := vector([u/(u^2+v^2+w^2), v/(u^2+v^2+w^2), w/(u^2+v^2+w^2)])

>    jacobian(A, [u,v,w]);

matrix([[1/(u^2+v^2+w^2)-2*u^2/(u^2+v^2+w^2)^2, -2*u/(u^2+v^2+w^2)^2*v, -2*u/(u^2+v^2+w^2)^2*w], [-2*u/(u^2+v^2+w^2)^2*v, 1/(u^2+v^2+w^2)-2*v^2/(u^2+v^2+w^2)^2, -2*v/(u^2+v^2+w^2)^2*w], [-2*u/(u^2+v^2+...

>    simplify(%);

matrix([[(-u^2+v^2+w^2)/(u^2+v^2+w^2)^2, -2*u/(u^2+v^2+w^2)^2*v, -2*u/(u^2+v^2+w^2)^2*w], [-2*u/(u^2+v^2+w^2)^2*v, -(-u^2+v^2-w^2)/(u^2+v^2+w^2)^2, -2*v/(u^2+v^2+w^2)^2*w], [-2*u/(u^2+v^2+w^2)^2*w, -2*...

>