function nummath(x)
%
% Berechnung der y Werte
%
[m, n] = size(x);
y = zeros(m, n);

% oder kurz
y = zeros( size(x) );

for i=1:m           % alternativ :   for i=1:size(y, 1)
    for j=1:n       % alternativ :   for j=1:size(y, 2)
        y(i,j) = sqrt(1+x(i,j))-1;
    end
end

%
% Berechnung der z Werte, koennte man wie bei y machen; oder kurz
% mit
z = x ./ (sqrt(1+x)+1);

%
% Plot ausgeben
%
loglog(x, y, 'rx');  % alternativen : plot, semilogy, semilogx
%hold on;
%loglog(x, z, 'b-');

%abschaetzung = eps*3/2;
%loglog(x, z+abschaetzung, 'k--', x, z-abschaetzung, 'k--');
%hold off;