Ya se mencionó en un comentario anterior sobre un tema en physics.SE, que puede estar relacionado con el problema de dispersión inversa. Sólo puedo añadir, que hay un método preciso de construcción del operador de Shroedinger $p^2 + V(q)$ con espectro finito arbitrario para el caso unidimensional. Está muy bien desarrollado debido a la aplicación a la teoría de solitones. Los operadores polinómicos $V(q)$ tienen algunos problemas debido a los infinitos para $q \to \pm \infty$ pero para $V(q)$ con rápida convergencia a cero hay cientos de trabajos. Una construcción especialmente sencilla puede encontrarse, por ejemplo, en M. Ablowitz, H. Segur, Solitons and the Inverse Scattering Transform, SIAM, Philadelphia, 1981.
La idea, que si tienes algo de $L_1 = p^2 + u(q)$ con un espectro dado, existe un método elegante para añadir otro valor propio $\zeta$ . Debería considerar la ecuación $L_1 g = \zeta g$ , encontrar la solución $g(q)$ . Entonces $L_2 = p^2 + u + 2 (\ln (g))''$ guarda todos los valores propios de $L_1$ pero también tiene un nuevo valor propio $\zeta$ con la función propia $1/g(q)$ . Así que puede empezar con $u=0$ y añadir los valores propios uno tras otro.
Sin embargo, como mencioné aquí hay un problema con el caso polinómico, es unidimensional y puede ser lo único útil es la idea de buscar el truco como $g \to 1/g$ con la adición de valores propios...
[ADD] (1). Creo que también es posible considerar los productos $L_k(q_k)$ etc. (2). No es único, por ejemplo, la deformación isoespectral de $L$ se describe mediante la ecuación de KdV (sin embargo, también se mencionó en una referencia en el hilo de Physics.SE)