function[]= Kreisbogen (Radius, P1,P2,P3,P4)
clc;
clear all;

%Punkte der Geraden
x1 = 4;
x2 = 1;
x3 = 3;
x4 = 1;
y1 = -2;
y2 = 2;
y3 = 6;
y4 = 2;
z1 = 0;
z2 = 0;
z3 = 0;
z4 = 0;

R=2;

P1 = [x1;y1;z1];
P2 = [x2;y2;z2];
P3 = [x3;y3;z3];
P4 = [x4;y4;z4];

%Matrix der Geradengleichungen
%Richtungsvektoren
m = P2-P1
n = P4-P3
b = P3-P1
A = [m -n]

%Lösbarkeit des LGS
if rank(A) == rank([A,b])
if rank(A) == size(A,2)
disp('Schnittpunkt')
x = A\b;
v = x(1)
s = x(2)
else
disp('hmm')
return
end
else
disp('Keinen Schnittpunkt')
return
end
% Schnittpunkt der Hauptgeraden
S = P1+ s*m


%Parallele Hilfsgeraden

%Kreuzprodukt der Normalvektor auf der Ebene h
h = cross(m,n)
%Normieren des Vektoren
h0 = h/sqrt(dot(h,h))
%Kreuzprodukt um Abstandsvektor zu bekommen
j = cross(h0,m)
%Normieren von j
j0 = j/sqrt(dot(j,j))

%Kreuzprodukt um Abstandsvektor zu bekommen
k = cross(h0,n)
%Normieren von k
k0 = k/sqrt(dot(k,k))


%Berechnung des Punktes P5, erster Punkt der erzeugten Hilfsgeraden h1
P5 = P1+j0*R

%Bestimmung des Punktes P6, erster Punkt der erzeugten Hilfsgeraden h2
P6 = P2+k0*R

%Berechnung des Schnittpunktes der Hilfsgeraden

c = P6-P5
w = A\c
% Sowohl Schnittpunkt Hilfsgeraden als auch Kreismittelpunkt 
M = P5 + w(1,1)*m

%Kreisbogenanfang K
e = P5-P3
l = A\e
K = P3 + l(1,1)*n

%Kreisbogenende E 
f = P6-P1
p = A\f

E = P1 + p(1,1)*m

% Die Geradengleichungen
g1 = P1 + s*m
g2 = P2 + v*n
h1 = P5 + s*m
h2 = P6 + v*n

% Grafische Ausgabe

%Das Zeichnen der Geraden

axis auto


x = g1
y = g2
t = 0:0.01:100;

axis([-10 10 -10 10 -10 10])
grid on
plot3(x,y,t,'b','LineWidth',2)

%title('plot3()','FontSize',20)
%xlabel('x_[(t)]','FontSize',20)
%ylabel('y_[(t)]','FontSize',20)
%zlabel('t','FontSize',20)

end