Jak wyznaczyć dzień tygodnia dla podanej daty z uwzględnieniem lat przestępnych?

(d+=m<3?y--:y-2,23*m/9+d+4+y/4-y/100+y/400)%7
0 - niedziela; 6-sobota
d - dzień
m - miesiąć od 1-12
y - rok

dzień tygodna = ([23m/9] + d + 4 + y + [z/4] - [z/100] + [z/400] - 2 (if m >= 3) ) mod 7

    gdzie
    [ ] oznacza dzielenie bez reszty z zaokrągleniem w dół
    mod - funkcja modulo
    m - numer miesiąca (od stycznia = 1 do grudnia = 12)
    d - numer dnia miesiąca
    z = y - 1 jeśli m < 3,
    z = y w pozostałych przypadkach