



'''simple docstring'''





import math
from datetime import datetime, timedelta

def 			a_	(				_lowerCAmelCase			) -> List[Any]:
  __lowerCamelCase			:			Any       			=					year % 19
  __lowerCamelCase			:			Dict       			=					year % 4
  __lowerCamelCase			:			Any       			=					year % 7
  __lowerCamelCase			:			List[Any]       			=					math.floor(year / 100			)
  __lowerCamelCase			:			List[str]       			=					math.floor((13 + 8 * leap_day_inhibits) / 25			)
  __lowerCamelCase			:			Optional[int]       			=					leap_day_inhibits / 4
  __lowerCamelCase			:			List[Any]       			=					(
      15 - lunar_orbit_correction + leap_day_inhibits - leap_day_reinstall_number
  ) % 30
  __lowerCamelCase			:			Any       			=					(4 + leap_day_inhibits - leap_day_reinstall_number) % 7

  # days to be added to March 21
  __lowerCamelCase			:			Union[str, Any]       			=					(19 * metonic_cycle + secular_moon_shift) % 30

  # PHM -> Paschal Full Moon
  __lowerCamelCase			:			Optional[int]       			=					(
      2 * julian_leap_year
      + 4 * non_leap_year
      + 6 * days_to_add
      + century_starting_point
  ) % 7

  if days_to_add == 29 and days_from_phm_to_sunday == 6:
    return datetime(__lowerCAmelCase			,4			,19			)
  elif days_to_add == 28 and days_from_phm_to_sunday == 6:
    return datetime(__lowerCAmelCase			,4			,18			)
  else:
    return datetime(__lowerCAmelCase			,3			,22			) + timedelta(
        days=int(days_to_add + days_from_phm_to_sunday			)			)


if __name__ == "__main__":
   for year in (1994, 2000, 2010, 2021, 2023):
      _UpperCamelCase											=		"""will be""" if year > datetime.now().year else """was"""
      print(f'''Easter in {year} {tense} {gauss_easter(year)}''')












import numpy as np
def 		a_							( __lowerCAmelCase   , __lowerCAmelCase   , __lowerCAmelCase   , __lowerCAmelCase   , __lowerCAmelCase				):
  lowerCAmelCase__										=      int(np.ceil((x_end - xa) / h				)				)
  lowerCAmelCase__										=      np.zeros((n + 1,)				)
  lowerCAmelCase__										=      ya
  lowerCAmelCase__										=      xa

  for k in range(__lowerCAmelCase				):
    lowerCAmelCase__										=      f(__lowerCAmelCase   , y[k]				)
    lowerCAmelCase__										=      f(x + 0.5 * h   , y[k] + 0.5 * h * ka				)
    lowerCAmelCase__										=      f(x + 0.5 * h   , y[k] + 0.5 * h * ka				)
    lowerCAmelCase__										=      f(x + h   , y[k] + h * ka				)
    lowerCAmelCase__										=      y[k] + (1 / 6) * h * (ka + 2 * ka + 2 * ka + ka)
    x += h

  return y


if __name__ == "__main__":
      import doctest

      doctest.testmod()

