Si en un instante inicial (t=0) la función de ondas es una superposición de estados estacionarios: Phi0 = suman [cn·phin] En cualquier otro instante t Phi_t = suman [cn·e-iEn·t/h·phin] donde nos tomamos la licencia de escribir h en lugar de h barra por cuestiones de edición. En el caso de un pozo de potencial, las funciones propias son funciones armónicas. Tomemos dos de estas funciones: phi1[i] = Math.sin(n1*Math.PI*x[i]/2); phi2[i] = Math.sin(n2*Math.PI*x[i]/2);

          psir[i] = c1*phi1[i] + c2*phi2[i] *Math.cos(n2*n2*t - n1*n1*t);
          psii[i] = - c2*phi2[i] *Math.sin(n2*n2*t - n1*n1*t);
          psipsi[i] = psir[i]*psir[i] + psii[i]*psii[i];

En el siguiente enlace se puede ver un ejemplo de una Simulación realizada con EJS.