The transonic small disturbance(TSD) theory has been used for calculation of aerodynamic forces in transonic region and a time accurate approximate factorization(AF) algorithm is formulated for the solution of the TSD equations. Steady and unsteady aerodynamic analysises have been performed for two dimensional airfoils and three dimensional wings. The flutter has been analyzed for pitching and plunging coupled system of two dimensional typical section using the coupled time marching method. Entropy effects have been incorporated within the solution procedure to analyze flows with strong shock waves more accurately. An Engquist-Osher(E-O) type-dependent mixed difference is implemented in the AF algorithm to treat supersonic regions of the flow. For more accurate unsteady applications, nonreflecting far-field boundary conditions are used. The developed TSD code gives reasonable results compared with the existing results.