direct injection at center point
ef[2*i-1, 2*j-1] = unc[i,j]