A method is outlined to obtain, directly in phase space, a complete set of eigenfunctions of the Liouville operator. The method is applied to quadratic lagrangians and the results are compared with the Wigner equivalents of the corresponding basic elements of the density matrix.