
import numpy as np





def 	SCREAMING_SNAKE_CASE__    (      __lowerCAmelCase    ,		__lowerCAmelCase    ,		__lowerCAmelCase    ,		__lowerCAmelCase    ,		__lowerCAmelCase							):
    snake_case__						=				int(np.ceil((x_end - xa) / h							)							)
    snake_case__						=				np.zeros((n + 1,)							)
    snake_case__						=				ya
    snake_case__						=				xa

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

    return y


if __name__ == "__main__":
     import doctest

     doctest.testmod()





def 	SCREAMING_SNAKE_CASE__    (      __lowerCAmelCase    ,		__lowerCAmelCase							):
    return price * (1 + tax_rate)


if __name__ == "__main__":
     print(F'''{price_plus_tax(100, 0.2_5) = }''')
     print(F'''{price_plus_tax(1_2_5.5_0, 0.0_5) = }''')



