% question5.m
% Solution to Computer Problem 1.2.14, Page 38
% Solution to question 5 of assignment 1

a = 0.0;
b = 1.0;
c = 1.0/sqrt(2.0);
d = 0.25;
e = 1.0; 
for k = 1 : 5 
   a = b; 
   b = (b+c)/2; 
   c = sqrt(c*a); 
   d = d - e*((b-a)^2); 
   e = 2*e; 
   f = (b^2)/d; 
   g = (b+c)^2/(4*d); 
   fprintf('%4d%25.15e%25.15e%25.15e%25.15e\n', ... 
      k, f, abs(f-pi), g, abs(g-pi)); 
end