% Session p. 174

session_p173

u1 = inline('richardson(tol,order,nmin,u,t)-1','t',...
                'tol','order','nmin','u');
options = optimset('TolX',tol);
t = fzero(u1,t,options,tol,order,nmin,u)

[val,err,ampl] = richardson(tol,order,nmin,u,t);
val, err = max(err,ampl*eps)