%% Wuerfel in 5 Tetraeder zerlegt
clear, clc, clf
% Koordinaten der Eckpunkte
KK = [-500 -500 -500;... 
       500 -500 -500;...
      -500 500 -500;...
       500 500 -500;...
      -500 -500 500;...
       500 -500 500;...
       -500 500 500;...
       500 500 500];

%  Zuordnung der Eckpunkte zu den Tetraedern (C-Indexing)
Tets = [ 0 1 2 4; 2 1 3 7; 1 4 5 7; 1 2 4 7; 4 2 6 7];
%Tets = [ 0 1 2 4; 2 1 3 7; 1 4 5 7; 1 2 4 7];

%  --> Matlab Indexing
Tets = Tets+1;

nelem = size(Tets,1);   % Anzahl der Tetraeder
nnode = size(KK,2);     % Anzahl der Eckpunkte

%%  Alle Tetraeder zeichnen 
figure(1)
tetramesh(Tets,KK,'FaceAlpha',0.3,'Visible','on')
xlabel('x'),ylabel('y'),zlabel('z')
title([num2str(nnode),' Punkte fuer ',num2str(nelem),' Tetraeder'])
view([-14,42])
% alpha(0.7)
axis equal

%% Tetraeder in Explosionsdarstellung
figure(2)
tetraeder_explosiv(Tets,KK, 0.3)
xlabel('x'),ylabel('y'),zlabel('z')
title(['Explosionsdarstellung: ',num2str(nnode),' Punkte fuer ',num2str(nelem),' Tetraeder'])
view([-14,42])
alpha(0.7)
axis equal

%% Tetraeder animiert in Explosionsdarstellung
figure(3)
tetraeder_animated(Tets,KK,1, 100)


