In the numerical modelling of the vacuum freeze drying of biomaterials (fodder yeast, lactose, milk biopreparation and sulphamerazine) at combined contact-radiant-microwave heating, the explicit MacCormack finite differences method was employed for solving two-region parabolic plane moving boundary problem (PPLMBP) with heat generation in both frozen and dried regions and unknown a priori sublimation temperature T-s(t) at moving ice sublimation front. Simulation of vacuum freeze drying process of biomaterials at combined contact-radiant-microwave heating was made, and additionally, the influence both steady capacity and relaxation time of internal heat source capacity has been discussed. A method of experimental determination of the mu(i) and tau(si) parameters was also presented.