Different methods have been proposed for history matching production of shale gas/oil wells which are drilled horizontally and usually hydraulically fractured with multiple stages. These methods are simulation, analytical models, and empirical equations. It has been well known that among the methods listed above, analytical models are more favorable in application to field data for two reasons. First, analytical solutions are faster than simulation, and second, they are more rigorous than empirical equations. Production behavior of horizontally drilled shale gas/oil wells has never been completely matched with the models which are described in this thesis. For shale gas wells, correction due to adsorption is explained with derived equations. The algorithm which is used for history matching and forecasting is explained in detail with a computer program as an implementation of it that is written in Excel's VBA. As an objective of this research, robust method is presented with a computer program which is applied to field data. The method presented in this thesis is applied to analyze the production performance of gas wells from Barnett, Woodford, and Fayetteville shales. It is shown that the method works well to understand reservoir description and predict future performance of shale gas wells. Moreover, synthetic shale oil well also was used to validate application of the method to oil wells. Given the huge unconventional resource potential and increasing energy demand in the world, the method described in this thesis will be the "game changing" technology to understand the reservoir properties and make future predictions in short period of time.