Abstract: Hamiltonian of a one-dimensional Bose-Hubbard model is re-formulated by using differential realization of the boson algebra. Energy matrices can then be generated systematically by using a Mathematica package. The output can be taken as the input of other diagonalization codes. As examples, exact energy eigenvalues and the corresponding wavefunctions for some cases are obtained with a Fortran diagonalization code. Phase transition of the model is analyzed.

Key words: insulator, superfluid, phase transition, phase diagram