By applying an improved  method of the Earth's crust classification, we develop an algorithm and build an earthquake prediction program using a combination of geological, geophysical, geomorphological and seismic data. This program includes a system of multiple windows with different functions, which can divide fault zones into the different segments by the maximum magnitude values Mmax. Using the constructed program, we carried out an earthquake prediction test for the Northwest Vietnam by a combination of geological, geophysical, geomorphological and seismic data. According to the received results, zoning maps of maximum earthquake prediction for the researched region has been established. The results show that the   areas, capable of generating earthquakes with Mmax = 6.5 - 6.8  are primarily concentrated along some major fault zones such as Lai Chau-Dien Bien, Son La, Song Ma, Song Da, Tuan Giao or near the intersection of these fault zones. The received results show a good accordance with the actual seismotectonic characteristics of the researched region.