data centes;
   year=1600;
   do until(year>20000);
      leap=mdy(2,29,year);
      output;
      year+100;
   end;
run;

proc print data=centes;
   title 'Centesimal Leap Years through the Year 20,000';
run;