sub p(y)
return mod((y + int(y/4) - int(y/100) + int(y/400)), 7)
end sub
sub isLongYear(y)
return (p(y) = 4) or (p(y - 1) = 3)
for y = 2000 to 2100
if isLongYear(y) print y
next y
end