How to show that the magnetisation M(H) is a concave (convex) function of the external field H on the interval (0, \infty) ( (-\infty, 0) ) in the ferromagnetic Ising model ( number of spins N is finite) ? Perhaps someone could suggest me how to approach this question   or where to find the answer.