An implementation for the algorithm of Hirota bilinear form of PDE in the Maple system. Based on the bilinear method of R. Hirota [The direct method in soliton theory. (2004; Zbl 1099.35111)], a reliable algorithm for generating a bilinear form of a class of nonlinear partial differential equations (PDEs) is described and a Maple package Bilinearization is established. The package will be helpful for looking for the bilinear form of nonlinear PDEs by a logarithm transformation.