% Script to test Romberg.m
% Example 8.9 on p.370 of Heath
% Romberg integration of sin(x) from 0 to pi.
% Exact result is given by 2.
% RombergTest.m
% K. Ming Leung, 02/09/03

clear all; format long;
exact = 2;
IRomberg = Romberg('sin',0,pi,6),
error = IRomberg(end,end)-exact,