%
% 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.
%
%-----

% EPage 103
% 

close all; clc; clear all; 

addpath('../../BookCode');

intdiv = 8;

% (i): 
fprime = @(t,x) ( 3 * x / t );
sol_i = rombergx( fprime, [1,20], intdiv, 1 );

% (ii): 
fprime = @(t,x) ( 2 * t * x );
sol_ii = rombergx( fprime, [0,2], intdiv, 2 );

fprintf('sol_i= %10.6f, sol_ii= %10.6f\n', sol_i, sol_ii);