def magnetisation_deviation(m_t, j):
    """
    Returns the deviation of the `j` component of the
    magnetisation from the starting magnetisation.

    """
    j = ord(j) - 120
    mj_t = m_t[:, j]
    delta_mj_t = mj_t - mj_t[0]
    return delta_mj_t