else if(j==-1 || j==n) b(id) -= w * boundary(i); // constrained coefficient insertCoefficient(id, i-1,j, -1, coefficients, b, boundary); insertCoefficient(id, i+1,j ...