Well, the fictive nodes are only necessary at the bondaries (upper ones), so yeah, when the laplace equation steered into a boundary point, it transformed the $u_{i+1,j}$ and $u_{i,j+1}$ variables to specific $u_{i,M+1}$ and $ u_{N+1,j}$ fictive points, which have their own special formulas...