Consolidation analysis of soils can be undertaken using empirical and numerical methods. Thernproblems related to empirical analysis of consolidation are the involvement of the determinationrnfor time factor using empirical equation obtained from curve fitting method. The empiricalrncorrelation of the time factor and the load application do not represent the actual situation wherernsoil properties are assumed to be spatially and temporally invariant and thus, average propertyrnvalues are used. In reality, however, these soil parameters usually vary from point to pointrn(heterogeneous) and these parameters may vary in time while a consolidation process isrnoccurring due to an external influence such as surface pressure. Moreover, it is difficult to modelrnand make consolidation analysis for variable loading conditions.rnIn order to solve this problem, numerical analysis using Finite Difference Method is proposed.rnThe finite difference method uses the original governing consolidation equation and provides thernsolution of a continuous problem by representing it in terms of a discrete set of elements suchrnthat there is a number of points in depth and time at which one calculates the field variables, inrnthis case the excess pore water pressure, consolidation settlement and average degree ofrnconsolidation.rnThus, this thesis involves development of user friendly application software for consolidationrnanalysis using finite difference method. In the hope of achieving this goal, a computer programrnnamed ConsolAnalysis has been developed. The analysis capabilities of the software comprisesrnof consolidation analysis for homogeneous and stratified soils subjected to constant loading,rnvariable loading, and abrupt change of loading with constant and variable coefficient ofrnconsolidation.