% % Written by: % -- % John L. Weatherwax 2009-04-21 % % email: wax@alum.mit.edu % % Please send comments and especially bug reports to the % above email address. % % Problem on EPage 62; Solutions on EPage 161 % %----- addpath('../../BookCode'); fn = @(x) ( x^1.4 - sqrt(x) + 1/x - 100 ); fn_prime = @(x) ( 1.4 * x^0.4 - 1/(2*sqrt(x)) - 1/x^2 ); fnewton(fn,fn_prime,50,1.e-4)