wolffd@0: function temp_schedule = compute_cooling_schedule(init_temp, final_temp, anneal_rate)
wolffd@0: 
wolffd@0: temp_schedule = [];
wolffd@0: i = 1;
wolffd@0: temp_schedule(i)=init_temp;
wolffd@0: while temp_schedule(i) > final_temp
wolffd@0:   i = i + 1;
wolffd@0:   temp_schedule(i)=temp_schedule(i-1)*anneal_rate;
wolffd@0: end
wolffd@0: 
wolffd@0:   
wolffd@0: 
wolffd@0: