Computation of plurigenera of a canonical threefold