Abstract: Based on the Biot’s consolidation theory, the foundation is regarded as transversely isotropic and the soil parameters are considered as random variables. Combined the finite layer theory with gradient optimization schemes, a formula of stochastic finite layer theory is deduced. A program is programmed based on the theory. The quadric-layered technique is used in the program, which makes the calculation more reasonable. An interface is designed using GUI in Matlab. Based on the interface technology of Matlab and Visual C++, a visual program is programmed using hybrid programming of VC++ and Matlab. Finally, the program is testet using an example. It is shown that the method is feasible; it is worth studying reliability analysis of soft foundation settlement.