A pre-computed importance function providing ex ante assessment of optimal placement of a single radioactive source is used for sequential placement of multiple radioactive sources. As each source is placed, the importance function may be modified to reflect the dose implemented by that source and all sources already placed. The greedy algorithm employed by this technique makes possible optimization of additional dimensions of treatment planning including, using directional and isotropic sources, using different source isotopes, using optimum source strengths, using source locations that are not on the regular grid and using a single movable source for variable times at different locations.