% test_fpi.m
% Try various functions for Fixed Point Iteration
% to root of x^3 + x - 1

g1 = @(x) 1.0 - x^3;
g2 = @(x) (1.0 - x)^(1.0/3.0);
g3 = @(x) (2.0*x^3 + 1)/(3*x^2 + 1);

fpi(g1,0.5,10)