computation.cpp setBoundaryX
loads of comments: iMax vs. iMax+1 etc.
(TB,AH)