Beginning with related basic concepts, the derivation of the blackbody radiation law by MATLAB is presented. It is assumed that the energy of photons is nhv (where n is an integer number, h is Plank constant and v is radiation frequency) in the existing theory. The Generalized Plank law and Stefan-Boltzman law with photon energy of αhv (0<α<1, 2, ...) are given. The possibility of photon splitting is discussed.