The main aim of this work is to contrast three various probabilistic computational techniques, namely analytical, simulation and perturbation-based, in a solution of the transient heat transfer problem in specific axisymmetric problem with Gaussian uncertainty in physical parameters. It is done thanks to a common application of the Finite Element Method program ABAQUS (for the deterministic part) and symbolic algebra system MAPLE, where all probabilistic procedures have been programmed. We determine up to the fourth order probabilistic characteristics of the resulting temperatures, i.e. expectations, coefficients of variation, skewness and kurtosis together with the histograms -all as the functions of the input coefficient of variation of random heat conductivity coefficient. Stochastic perturbation technique is implemented here using the tenth order Taylor series expansion and traditional Least Squares Method released with polynomial basis whose final order is a subject of the separate statistical optimization. Probabilistic results computed show almost perfect agreement of all the probabilistic characteristics under consideration, which means that the traditional simulation method may be replaced due to the time and computer scale savings with the stochastic perturbation method.