Everything is finally consistent now. To use original programm, use solve_poisson.m_original and poisson solving method 1.