%
%  Berechnung von Mittelwert und Standardabweichung  
%    für übergebene Messreihe g 
%   >> [mw,sa,sf] = statistik( [8,9,10,11,12] ) 
%
function [mw,sa,sf] = statistik( g )

   n = length( g );   
   fprintf( 'Das Feld g hat %d Elemente. \n', n );
   
   % Summe der Elemente
   s = 0;
   for( k=1:n )
      s = s + g(k); 
   end
   
   mw = s/n;   % Mittelwert
   
   % Summe der quadratischen Abweichungen
   sq = 0;
   for( k=1:n )
      sq = sq + ( g(k) - mw )^2; 
   end
   
   sa = sqrt( sq/(n-1) );  % Standardabweichung 
   sf = sa / sqrt(n);      % Standardfehler 
   