wolffd@0: function schedule = mk_schedule(init_temp, final_temp, anneal_rate)
wolffd@0: 
wolffd@0: init_temp = 10; final_temp = 1e-2; anneal_rate = 0.8;
wolffd@0: schedule = [];
wolffd@0: temp=init_temp;
wolffd@0: schedule = [schedule temp];
wolffd@0: while temp > final_temp
wolffd@0:   temp = temp * anneal_rate;
wolffd@0:   schedule = [schedule temp];
wolffd@0: end