Code covered by the BSD License  

Highlights from
MATLAB en zijn Symbolic Math Toolbox

image thumbnail

MATLAB en zijn Symbolic Math Toolbox

by

Tjibbele Miedema

 

Leergang voor geïnteresseerden in de Computeralgebra - Systemen MATLAB en zijn Symbolic Math Toolbox

PiA.m
%
%               *****
% Script-file : PiA.m
%               *****
%
% Bepaling van pi uit de omtrek van regelmatige in- en omgeschreven 4, 8, ... ,
% veelhoeken
%

i   = 1;
n   = 4;         % initialisatie van een in- en omgeschreven vierkant
s   = sqrt(2.0);
c   = 2.0;
pil = n * s / 2;
piu = n * c / 2;

disp('---------------------------------------------------------------------------')
disp('|  i       n              s                 pil                piu        |')
disp('---------------------------------------------------------------------------')

disp(sprintf('| %2d  %10d  %17.15f  %17.15f  %17.15f |', i, n, s, pil, piu))

for i = 2 : 29
   n  = n * 2;
   
   k  = 1;
   d  = (s / 2.0)^2;
   sa = d;
   
   while(1)
      d  = (2 * k - 1) / (k + 1) * (s / 2.0)^2 / 2.0 * d;
      k  = k + 1;
      sa = sa + d;
      if d < 0.5e-16
         break;
      end % if
   end % while
   
   s   = sqrt(sa);
   c   = 2.0 * s / sqrt(4.0 - s^2);
   pil = n * s / 2;
   piu = n * c / 2;
   
   disp(sprintf('| %2d  %10d  %17.15f  %17.15f  %17.15f |', i, n, s, pil, piu))
   
end % for

disp('---------------------------------------------------------------------------')

%
%               *****
% Script-file : PiA.m
%               *****
%

Contact us