Only for keep working, check in the attached yellow and red regions. I think that the error in the discretization was from considerating a matrix Vectorincog2 with dimensions nE,nZ,nq. In mathcad you can't do that, mathcad only have 2 dimemnsional matrices. After that you assign FefleEHK1D to this function and with this use it in GK1D but with dimensions nE,nq. So, there are an inconsistence somewhere there.
Hope attached can be useful in some things, but can be better that another collab take this for a more physical interpretation, not only the mathcad implementation.
Regards. Alvaro.
Solo para seguir trabajando, verifica en el archivo adjunto las regiones amarillas y rojas. Creo que el error en la discretización proviene de considerar la matriz Vectorincog2 con dimensiones nE,nZ,nq. No puedes hacer eso en mathcad, mathcad s{olo tiene matrices bidimensionales. Luego de eso asignas FefleEHK1 a esa función y lo usas en GK1D, pero con dimensiones nE,nq. Es decir, hay alguna inconsistencia por ahí.
Espero que el archivo adjunto sea útil para algunos de las cuestiones, pero sería mejor que otro colaborador tome este para una interpretación más física, y no sólo la implementación en mathcad.
PD:
In other sort of things, we have a discussion here about units, your file handles the probability but in cartesian form. It can be done also in polar notation; sometimes math can simplify in this way. There are two paths where the units can be extended -as Philip propose- with complexes. One is with the electromagnetic field, and the other is assignin units to complexes in such way that the amplitude of the probabilities fails in the exponentianl of the polar notation of the complexes.
