%% Kreis und Kartesisches Blatt

clear; clc;

%% Kreis
n   = 20;           % Anzahl der Unterteilungen
rad = 3;            % Radius
M   = [2 1];         % Kreismittelpunkt

phi = linspace(0,2*pi,n+1);

[x y] = pol2cart(phi,rad);

x = M(1)+x;
y = M(2)+y;

plot(x,y); axis equal;
pause()


%% Kartesisches Blatt
clear; clc; clf;

n   = 200;           % Anzahl der Unterteilungen
a   = 5;             % Parameter für Kart. Blatt 
phi = linspace(-pi/4,3*pi/4,n+1); % !! Vorzeichenwechsel in x,y für -pi/4 !!

t = tan(phi);

x = 3*a*t./(1+t.^3);
y = 3*a*t.^2./(1+t.^3);

plot(x,y); grid;axis equal; xlim([-10 10]); ylim([-10 10]);



disp(['Darstellung mit n=',num2str(n)]);