for each iteration, for each patch i

for each patch j:

calculate the formfactors Fij using hemi-cube at patch i

/*update change since last time patch j shot light */
/* update total radiosity of patch j */
/* reset unshot radiosity for patch i to zero */

initialization:

for all patch i:

if patch i is a light source, then , else