Easter Sunday is the first Sunday after the first full moon of spring. To compute the date, you can use this algorithm, invented by the mathematician Carl Friedrich Gauss in 1800: 1. Let y be the year (such as 1800 or 2001). 2. Divide y by 19 and call the remainder a. Ignore the quotient. 3. Divide y by 100 to get a quotient b and a remainder c. 4. Divide b by 4 to get a quotient d and a remainder e. 5. Divide 8 * b + 13 by 25 to get a quotient g. Ignore the remainder. 6. Divide 19 * a + b - d - g + 15 by 30 to get a remainder h. Ignore the quotient. 7. Divide c by 4 to get a quotient j and a remainder k. 8. Divide a + 11 * h by 319 to get a quotient m. Ignore the remainder. 9. Divide 2 * e + 2 * j - k - h + m + 32 by 7 to get a remainder r. Ignore the quotient. 10. Divide h - m + r + 90 by 25 to get a quotient n. Ignore the remainder. 11. Divide h - m + r + n + 19 by 32 to get a remainder p. Ignore the quotient. Then Easter falls on day p of month n. For example, if y is 2001:

# The user is prompted to enter the year

# the entered year is assigned to y

y = int(input("Enter the year: "))

# the remainder when y is divided by 19 is assigned to a

a = y % 19

# The quotient is assigned to b when y is divided by 100

b = y // 100

# The remainder is assigned to c when y is divided by 100

c = y % 100

# The quotient is assigned to d when b is divided by 4

d = b // 4

# The remainder is assigned to e when b is divided by 4

e = b % 4

# The quotient is assigned to g when 8 * b + 13 is divided by 25

g = (8 * b + 13) // 25

# The remainder is assigned to h when 19 * a + b - d - g + 15 is divided by 30

h = (19 * a + b - d - g + 15) % 30

# The quotient is assigned to j when c is divided by 4

j = c // 4

# The remainder is assigned to k when c is divided by 4

k = c % 4

# The quotient is assigned to m when a + 11 * h is divided by 319

m = (a + 11 * h) // 319

# The remainder is assigned to r when 2 * e + 2 * j - k - h + m + 32 is divided by 7

r = (2 * e + 2 * j  - k - h + m + 32) % 7

# The quotient is assigned to n when h - m + r + 90 is divided by 25

n = (h - m + r + 90) // 25

# The remainder is assigned to p when h - m + r + n + 19 is divided by 32

p = (h - m + r + n + 19) % 32

# The date of easter is displayed

print ("Easter falls on day ", p, "of month ", n);

Explanation:

The code is written in Python and is well commented. A sample image of program output is attached.