For the first time, we compute the quantum discord in bipartite systems containing up to nine qubits. An analytical expression is obtained for the discord in a bipartite system with three qubits. The dependence of the discord on the temperature and the structural parameter of the model is studied.
