% % 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 % %----- close all; clc; clear all; addpath('../../BookCode'); e = 0.96727464; M = 4.52759e-3; fn = @(E) ( E - e*sin(E) - M ); fn_prime = @(E) ( 1 - e*cos(E) ); x_root = fnewton( fn, fn_prime, 1, 0.00005 )