>
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)] );
>
jacobian(A, [u,v,w]);
>
simplify(%);
>