This paper presents a transmission probability method(TPM) to solve the neutron transport equation in three-dimensional triangularz geometry.The source within the mesh is assumed to be spatially uniform and isotropic.At the mesh surface,the constant and the simplified P1 approximation are invoked for the anisotropic angular flux distribution.Based on this model,a code TPMTDT is encoded.It was verified by three 3D Takeda benchmark problems.The results of the present method agree well with those of Monte-Carlo calculation method and Spherical Harmonics method.