def calculate_dc_coefficients(contour):
dxy = np.diff(contour, axis=<NUM_LIT:0>)<EOL>dt = np.sqrt((dxy ** <NUM_LIT:2>).sum(axis=<NUM_LIT:1>))<EOL>t = np.concatenate([([<NUM_LIT:0.>, ]), np.cumsum(dt)])<EOL>T = t[-<NUM_LIT:1>]<EOL>xi = np.cumsum(dxy[:, <NUM_LIT:0>]) - (dxy[:, <NUM_LIT:0>] / dt) * t[<NUM_LIT:1>:]<EOL>A0 = (<NUM_LIT:1> / T) * np.sum(((dxy[:, <NUM_LIT:0>] / (<NUM_LIT:2> * dt)) * np.diff(t ** <NUM_LIT:2>)) + xi * dt)<EOL>delta = np.cumsum(dxy[:, <NUM_LIT:1>]) - (dxy[:, <NUM_LIT:1>] / dt) * t[<NUM_LIT:1>:]<EOL>C0 = (<NUM_LIT:1> / T) * np.sum(((dxy[:, <NUM_LIT:1>] / (<NUM_LIT:2> * dt)) * np.diff(t ** <NUM_LIT:2>)) + delta * dt)<EOL>return contour[<NUM_LIT:0>, <NUM_LIT:0>] + A0, contour[<NUM_LIT:0>, <NUM_LIT:1>] + C0<EOL>
Calculate the :math:`A_0` and :math:`C_0` coefficients of the elliptic Fourier series.

    :param numpy.ndarray contour: A contour array of size ``[M x 2]``.
    :return: The :math:`A_0` and :math:`C_0` coefficients.
    :rtype: tuple