% Vektor der Winkelparameter

phi_vec = 0:pi/200:2*pi;

% parameter
a = 1; h = 1.5;


% Vektor der Radius-werte
r_vec = 2*a*cos(phi_vec) + h;

% Vektoren der x- und y-Koordinaten
x_vec = r_vec .* cos(phi_vec);
y_vec = r_vec .* sin(phi_vec);


% Kurve Zeichnen

figure(1)

clf

axis equal
axis([-(2*a+h),(2*a+h),-(2*a+h),(2*a+h)])

hold on

p = line(x_vec(1),y_vec(1),'Color','blue','Erasemode','xor','LineStyle','none','Marker','+');

for k = 2:length(phi_vec)
    
    plot([x_vec(k-1),x_vec(k)],[y_vec(k-1),y_vec(k)],'Color','blue','Erasemode','none','Linewidth',2);

    set(p,'Xdata',x_vec(k),'Ydata',y_vec(k));
    
    drawnow

end