The use of high dose rate remote afterloading system for the treatment of intraluminal lesions necessitates the need for a more accurate of dose distributions around the high intensity brachytherapy sources, doses are often prescribed to a distance of few centimeters from the linear source, and in this range the dose distribution is very difficult to assess. Accurated and optimized dose calculation with stable numerical algorithms by PC level computer was required to treatment intraluminal lesions by high dose rate brachytherapy system. The exposure rate from sources was calculated with Sievert integral and dose rate in tissue was calculated with Meisberger equation. An algorithm for generating a treatment plan with optimized dose distribution was developed for high dose rate intraluminal radiotherapy. The treatment volume becomes the locus of the constrained target surface points that is the specified radial distance from the source dwelling positions. The treatment target volume may be alternately outlined on a x-ray film of the implant dummy sources. The routine used a linear programming formulism to compute which dwell time at each position to irradiate the constrained dose rate at the target surface points whiles minimizing the total volume integrated dose to the patient. The exposure rate and the dose distribution to be confirmed the result of calculation with algorithm were measured with film dosimetry, TLD and small size ion chambers.