




'''simple docstring'''







from math import asin, atan, cos, radians, sin, sqrt, tan

lowerCAmelCase_     						=				6_37_81_37.0
lowerCAmelCase_     						=				6_35_67_52.31_42_45
lowerCAmelCase_     						=				6_3_7_8_1_3_7







def 						_A      ( UpperCAmelCase      ,UpperCAmelCase      ,UpperCAmelCase      ,UpperCAmelCase						):






   '''simple docstring'''







   A__       =     (AXIS_A - AXIS_B) / AXIS_A
   A__       =     atan((1 - flattening) * tan(radians(UpperCAmelCase						)						)						)
   A__       =     atan((1 - flattening) * tan(radians(UpperCAmelCase						)						)						)
   A__       =     radians(UpperCAmelCase						)
   A__       =     radians(UpperCAmelCase						)
   # Equation
   A__       =     sin((phi_a - phi_a) / 2						)
   A__       =     sin((lambda_a - lambda_a) / 2						)
   # Square both values
   sin_sq_phi *= sin_sq_phi
   sin_sq_lambda *= sin_sq_lambda
   A__       =     sqrt(sin_sq_phi + (cos(UpperCAmelCase						) * cos(UpperCAmelCase						) * sin_sq_lambda)						)
   return 2 * RADIUS * asin(UpperCAmelCase						)


if __name__ == "__main__":
    import doctest

    doctest.testmod()








'''simple docstring'''







def 						_A      ( UpperCAmelCase						):






   '''simple docstring'''







   A__       =     0
   while num > 0:
      digit_sum += num % 10
      num //= 10
   return digit_sum







def 						_A      ( UpperCAmelCase = 100						):






   '''simple docstring'''







   A__       =     1
   A__       =     2

   for i in range(2      ,max_n + 1						):
      A__       =     pre_numerator
      A__       =     2 * i // 3 if i % 3 == 0 else 1
      A__       =     cur_numerator
      A__       =     e_cont * pre_numerator + temp

   return sum_digits(UpperCAmelCase						)


if __name__ == "__main__":
    print(f'''{solution() = }''')


