function A = apfel(NX, NY, maxit)

% Berechnet ein in der Matrix A kodiertes Bild des Apfelmaennchens
% im Bildbereich [-2.5,1.5]+i[-1.5,1.5].
% Eingabewerte:   NX        Anzahl der Gitterpunkte in x-Richtung
%                 NY        Anzahl der Gitterpunkte in y-Richtung
%                 maxit     maximale Iterationszahl
% Ausgabewert:    A         Matrix mit Farbwerten der entsprechenden Gitterpunkte



A=zeros(NY,NX);           % Initialisierung der Matrix A
x=linspace(-2.5,1.5,NX);  % Vektor der x-Koordinaten
y=linspace(-1.5,1.5,NY);  % Vektor der y-Koordinaten


for k=1:NX
   for j=1:NY          % Schleife ueber die Gitterpunkte c
       c=x(k)+i*y(j);
       z=c;            % Startwert der Iteration
       for m=1:maxit
           z=z*z+c;    % Iteriere
           if (norm(z)>2)  % Abbruch der Iteration, falls |z|>2
               A(j,k)=m;   % Setze Farbwert des entspr. Gitterpunktes
               break;
           end
       end
   end
end 