f4(x) = x^(-2) - sin(x);
df4dx = x -> -2*x^(-3) - cos(x);
plot(f4, 0.5, 10)
inits = 1.:7.
x = zeros(length(inits))
for (i,y) in enumerate(inits)
f = x -> f4(x) - y
dfdx = x -> df4dx(x)
r = FNC.newton(f,dfdx,y)
x[i] = r[end]
end
pretty_table([inits x])