The internal documentation on NDEigenvalues is a bit unclear. I want to solve two coupled eigenvalue equations in 2D for two coupled eigenfunctions. The problem is of the form:
LinOp1[ u[x,y], v[x,y] ] == eval u[x,y] LinOp2[ u[x,y], v[x,y] ] == eval v[x,y]
LinOp2 are linear combinations of derivatives of the functions
v[x,y]. The documentation implies that this is possible, because it discusses coupled systems with a syntax giving the list of linear differential operators and a list of functions. But it's not clear in the documentation how Mathematica "knows" the eigenfunction on the RHS for the first operator is
u[x,y], and the second one is
v[x,y]. Is it just a matter of the order in the list of arguments? Can anyone point me to a sample of a coupled system solved with NDEigenvalues?