A polynomial function having eight constants and a non-polynomial analytical function having six constants are fitted to the observed mean vertical profile of the monsoonal zonal wind over Indian region. The constants are related to the observed position and magnitude of the upper level easterly and the lower level westerly jet streams, the level of the wind reversal and the surface wind speed. When the constants are realistically chosen, then the analytical profiles agree very well with the observed profiles, in detail. Further, an iterative scheme is suggested for improving the non-polynomial representation.