We construct an ensemble of second-quantized Hamiltonians with two bosonic degrees of freedom, whose members display with probability one GOE or GUE statistics. Nevertheless, these Hamiltonians have a second integral of motion, namely the boson number, and thus are integrable. To construct this ensemble we use some "reverse engineering" starting from the fact that n-bosons in a two-level system with random interactions have an integrable classical limit by the old Heisenberg association of boson operators to actions and angles. By choosing an n-body random interaction and degenerate levels we end up with GOE or GUE Hamiltonians. Ergodicity of these ensembles completes the example.PACS numbers: 05.45. Mt, 05.30.Jp, 31.15.Gy Recently, there has been considerable interest in spinless n-boson systems with k-body random interactions [1,2], both for degenerate [3] and nondegenerate [4, 5] single-particle levels. In particular, anomalous statistics for the two-level system have been understood from the fact that these systems are integrable [6]. The two-level ensemble corresponds to a time-independent two degrees of freedom system, in the sense that creation and annihilation operators for the two single-particle levels are canonical operators for the system. The second integral of motion corresponds to the conservation of the number of bosons. The Hamiltonian is a function of the creation and annihilation operators, whose precise form depends on the question whether we have two-body interactions or more complicated manybody interactions. The number of (random) coefficients in this model increases quadratically with the rank k of the interaction. A classical Hamiltonian can formally be written for any number of particles but in fact only for large particle number the quantum problem reaches the classical limit [7,8]. Thus the boson number plays the role of an action. For fixed particle number, the model can be reduced to a Hamiltonian of one degree of freedom with the number of bosons appearing as parameter.In this paper we shall consider a different large-n limit for the two-level system. Our model consists in choosing the rank of the interaction equal to the particle number, k = n. Hence the Hamiltonian changes for different values of the particle number. This choice, by definition of the ensemble, directly leads to a GOE or GUE matrix [9]. We then consider the limit of large particle number. Clearly, we will have an ensemble of systems which have GOE or GUE statistics with probability one [10], while being integrable in a well-defined limit of large actions. This represents an important caveat concerning the idea that Wigner-Dyson statistics are characteristic of classical chaotic systems. Note that this has no bearing on the fairly-well established quantum chaos conjecture [11,12,13,14,15,16,17], which establishes that classical chaos implies typically Wigner-Dyson statistics.In a similar vein, Wu et al. [18] have performed the following calculation: after having chosen an unfolded fixed spectrum of given length ...