# function [y] = f_sin_x_over_xsqp1(x) # Accompanying program for the text # # Classical and Modern Numerical Analysis: # Theory, Methods and Practice # by Azmy S. Ackleh, Edward J. Allen, # R. Baker Kearfott, and Padmanabhan Seshaiyer # # (Taylor and Francis / CRC Press, 2009) import numpy as np def f_sin_x_over_xsqp1(x): return np.sin(x) / (x**2 + 1)