f_sin_x_over_xsqp1.py

# 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)