We build and analytically calculate the Generalised Gibbs Ensemble
partition function of the integrable Soft Neumann Model. This is the
model of a classical particle which is constrained to move, on average
over the initial conditions, on an NN
dimensional sphere, and feels the effect of anisotropic harmonic
potentials. We derive all relevant averaged static observables in the
(thermodynamic) N\to\inftyN→∞
limit. We compare them to their long-term dynamic averages finding
excellent agreement in all phases of a non-trivial phase diagram
determined by the characteristics of the initial conditions and the
amount of energy injected or extracted in an instantaneous quench. We
discuss the implications of our results for the proper Neumann model in
which the spherical constraint is imposed strictly.