One of the major subsystems of each airplane is landing gear system which must be capable of tolerating extreme forces applied to the airplane during landing. Using conservative techniques to find landing loads result in overestimation and unnecessary extra structural weight. New commercial softwares can simulate real landing conditions with acceptable accuracy if detailed mechanical data about landing gear system subparts are provided. Although these softwares work well but due to lack of detailed information about the subparts at the conceptual design phase, complexity and time consuming of modeling, expensive license price, etc. they do not seem to be the best choice for design purpose. In this study, in order to calculate landing loads more precisely than the estimating conservative methods, flight dynamic differential equations of an airplane during landing phase are derived and through numeric and state space techniques are solved for different initial conditions including, three point landing, two point landing and one wheel landing. Each landing gear of the airplane is modeled as a two-degree of freedom mass-spring-damper set. Time history of the airplane center of gravity, pitch and roll angle, vertical landing loads of each landing gear and their spin-up loads for different landing types (different initial conditions) are obtained to show capabilities of this new, fast and accurate landing simulation code, generated.