source('utils.R')

##
## Test the function 'table':
##
test_function_table()


## 2.4-1:
##

## Part (a):
##
xs = c(-10, +10)
ys = xs^2
xbar = 0.0
res = table(xbar, xs, ys, tol=1e-5, verbose=TRUE)
print(abs(res-0))

## Part (b):
##
xs = c(-100, 0, +100)
ys = xs^3
xbar = -50.0
res = table(xbar, xs, ys, tol=1e-5, verbose=TRUE)
print(abs(res-xbar^3) )


## 2.4-4:
##
## Generate the direct data:
##
xs = seq(0.2, 0.8, length.out=20)
ys = sin(xs)

## Find sin(xbar) = 0.6
##
ybar = 0.6
res = table(0.6, ys, xs, verbose=TRUE)
print(sprintf('xbar= %f', res))