close all; clc; clear; 

format long;

% Problem EPage 82
% 
k = 100;
fn = @ (x) ( sin( x ) .* cos( k * x ) );
simpsons_approx = simp1( fn, 0, pi, 1024 )
rhom_approx = romb( fn, 0, pi, 9 )

format short;