I posted my question at MS but unfortunately it is still without a response, so let me ask it here. We can think about a bounded operator $T\colon c_0\to c_0$ as a double-infinite matrix $[T_{mn}]_{m,n\geq 1}$ which acts on a sequence $a=[a_1, a_2, a_3, \ldots ]\in c_0$ in the same way as usual (finite) matrices act on vectors ($n$-tuples of scalars), i.e. $$ Ta= [T_{mn}][a_n] = \left[ \sum_{n=1}^\infty T_{mn}a_n\right] $$ Suppose $a=[a_1, a_2, a_3, \ldots ]\in \ell^\infty = (c_0)^{**}$. Does the following formula still hold: $$ T^{**}a= [T_{mn}][a_n] = \left[ \sum_{n=1}^\infty T_{mn}a_n\right] $$ Are there any methods of computing the values of second conjugate operators acting on biduals $X^{**}$ of classical (non-reflexive) spaces $X$ at elements from $X$?