function a = f_prime(x)
    a = (f(x + 0.001) - f(x))/0.001;
end