Thanx for quick replay,
That somewhat clarifies my problem.
However I have particles in continuum, in the presence of some potential, that are described by Coulomb (Coulomb-like to be more precise) wave function that does depend on both, wave vector and position at the same time...
The wave function for fermions has to be anti-symmetric with respect to exchange of positions of electrons, but what if it depends on wave vector as well. Does they have to be exchanged as well, in other words, for two-electron system what is correct
Ψ(r1,k1,r2,k2) = - Ψ(r2,k1,r1,k2)
or...
You should use "eqnarray" instead of "equation". Eg.
\begin{eqnarray}
blah blah blah 1 \nonumber\\
blah blah blah 2 \nonumber\\
...
blah blah blah end
\end{eqnarray}
As it's written in book that you referred to, you have to install (if you don't have already) wasysym package and include at start of your document (\usepackage{wasysym}).
You should avoid to letting Mathematica works with lots of symbols. Just absorb exponential factor in constant b and redefine it later.
Clear[a,b,w,g,t]
ρt={{a,b},{b*,1-a}};
{e1,e2}=Eigenvalues[ρt];
{v1,v2}=Eigenvectors[ρt];
Here you return exponential factor since it is conugated...
You have to generate factors in reverse order. Just change arguments in Table command:
term[n_Integer,m_Integer]:=Apply[Dot,Table[matg[[i]],{i,n,m+1,-1}]].hl[[m]]/;(n>=m);
Subscript[lay,s_]:=Sum[term[s,i],{i,s}];
Remove that "Table" from import, use something like this:
plotdata = Import["test2 10 10.csv"]
z = GatherBy[plotdata, First][[All, All, 3]]
ArrayPlot[z, ColorFunction -> "Rainbow"]
Here you have it explained:
http://www.phy.pmf.unizg.hr/~npaar/teaching/compphys.pdf
page 292, "13.4 More on finite difference methods, Runge-Kutta methods"