Tissue reactions to cryosurgical procedures depend on temperature variations and on rates of temperature variations induced by freezing probes. In this paper thermal responses of biological systems undergoing freezing are obtained through the application of the finite element method to the solution of the nonlinear bio-equation. This model allows realistic predictions of isotherm fields and of rates of freezing in practically any cryosurgical procedure.