def fill_jacobian_column(self, jaccol, coordinates):
q, g = self.icfn(coordinates[list(self.indexes)], <NUM_LIT:1>)<EOL>for i, j in enumerate(self.indexes):<EOL><INDENT>jaccol[<NUM_LIT:3>*j:<NUM_LIT:3>*j+<NUM_LIT:3>] += g[i]<EOL><DEDENT>return jaccol<EOL>
Fill in a column of the Jacobian.

           Arguments:
            | ``jaccol`` -- The column of Jacobian to which the result must be
                            added.
            | ``coordinates`` -- A numpy array with Cartesian coordinates,
                                 shape=(N,3)