I've run some more tests and have found the answer to questions 2 and 3. Mathematica does seem to be renormalising with each value of $x$ and $y$. Here's a graph of the norm of one of the eigenvectors:
Quiet[ Plot3D[ Sqrt[Abs[Eigenvectors[m][[2]][[1]]]^2+Abs[Eigenvectors[m][[2]][[2]]]^2]-1, {x,x1,x2},{y,y1,y2}, Exclusions->None, PlotLabel->Row[{"|",Subscript[v,2],"|"}], PlotRangePadding->0 ] ] Given that (as shown on my original post) the norm of one of the components of this vector is close to 1 in a large region of the domain, I'll ascribe the to machine precision (i.e. the norm of this eigenvector is 1 throughout the domain). This seems to imply that Mathematica is not only renormalising the eigenvectors at each point, but it is actually giving me the eigenvectors normalised to 1 (contrary to what the documentation says).
Of course, this assumes that the order of the eigenvalues and eigenvectors is preserved regardless of which eigenvalue has a larger absolute value at each point. If it isn't, my whole interpretation of the above graph automatically goes out the window.
Edit:
I've found the answer to question 1 by testing a little further by building a matrix with eigenvalues whose ordering should change in a very obvious way:
M=DiagonalMatrix[{x-2,2-x2}]; Graphing the eigenvalues themselves (which are obviously $x-2$ and $2-x$) yields the following (naturally):
Plot[{2-x,x-2},{x,0,4}] However, plotting the eigenvalues as eigenvalues yields the following:
Plot[{Eigenvalues[M][[1]],Eigenvalues[M][[2]]},{x,0,4}] It is as I feared: Eigenvalues are reordered at each point. Clearly, then, eigenvectors are also reordered at each point, since the eigenvalue and eigenvector ordering must be consistent at all times (at least according to Mathematica).
Other than having explicit expressions for the eigenvector components (which I do not have, as explained at the beginning of my original post), is there a simple way around this reordering?


