We present a method for the direct computation of the Wigner function
by solving a coupled system of linear partial differential equations i
n phase space. Our modified spectral method relies on Chebyshev polyno
mials. Since this approach allows us to include arbitrary high orders
of partial derivatives, our procedure is applicable to arbitrary bindi
ng potentials. We apply our scheme to Wigner functions of the harmonic
oscillator, the Morse oscillator, and an asymmetric double-well poten
tial.