We present numerical methods to solve the Israel-Stewart (IS) equations of causal relativistic dissipative fluid dynamics with bulk and shear viscosities. We then test these methods studying the Riemann problem in (1+1)-and (2+1)-dimensional geometry. The numerical schemes investigated here are applicable to realistic (3+1)-dimensional modeling of a relativistic dissipative fluid.