One of cancer's therapeutic methods is applying ultrasound induced hyperthermia independently or in combination with other therapeutic methods like radiotherapy and chemotherapy. Because of improper distribution of sonication in tissue under treatment and not achieving to a suitable treatment planning during sonication, tissue temperature do not regulate correctly. So in this study to achieve ultrasonic treatment planning protocol, Rayleigh integral has been used for modeling propagation of sound. In this study, Rayleigh integral for two 1 MHz ultrasonic transducers have been solved in cylindrical coordinates for the purpose of estimating effects of 0.8 and 5.0 cm² effective radiation area and different intensities. Ultrasonic pressure distribution in axial, radial and angular directions were extracted. For validation of modeling results, acoustic pressure value in aqueous medium and at specific distances from transducer face has been measured using piston hydrophone. The results show that there are good agreement between modeling results and experimental measurements with 0.82 correlation of coefficient (p-value<0.05). The mean relative error is 7.34 percent for 5.0 cm² effective radiation area with 1.5 to 14.2 percent interval. It is concluded that this model is capable of modeling sound propagation in tissue and performing first step in hyperthermia treatment planning with acceptable accuracy.