Stiff differential equations frequently arise in physical problems due to the existence of greatly differing time constants. One way of solving stiff differential equations is taking the step size very small, another way is using an A-stable or stiffly stable implicit method. The purpose of this thesis is to discuss the method of determining initial values fo r stiff systems. The algorithm was motivated by a problem from ionospheric physics in which some of the initial values are unknown. If we give those unknown initial values an arbitrary number and integrated with those values, we have an initial transient. However, if we integrate with the correct initial values then there is no initial transient.