A program for the application described in subsection REF 
including several extra features (graphical output etc.) has been realized in
 3400 lines of Python 3 using the NumPy[1]} and
Matplotlib[2]} packages. The algorithms presented above
contain several loops which could be parallelized by using the
multiprocessing package: the first \(j\)  loop and the \(i\)  loop in Algorithm
REF , the \(j\)  loop in Algorithm REF , the \(i\)  loop in
Algorithm REF  and the \(c\)  loop in
REF . We keep this optimization for future versions of
this implementation.
