%Lab 9
a=[12,16];
bracketing = fzero(inline('(9.81*68.1/x)*(1-exp(-x*10/68.1))-40'),a)
guess = fzero(inline('(9.81*68.1/x)*(1-exp(-x*10/68.1))-40'),5)
 
 
 
>> example1
bracketing =
 14.8011
guess = 
 14.8011
 
 
function [hi=h,Myarea] = trap(func,p1,p2,n)
 
if nargin < 4 
 n = abs(p2-p1); 
end
i=1;
x=linspace(p1,p2,n+1)
h=(p2-p1)/n
areatemp=0;
while i>0
 i=i+1;
 areatemp = areatemp + h* ( func(x(i-1)) + func(x(i)) ) /2;
 
 if i == n+1
 break
 end
end
hi=h;
Myarea=areatemp;