%% Zeitpunkt t0
fid = fopen('QuaderT0.txt');
S = textscan(fid, '%s %f32 %f32 %f32 %f32 %f32 %f32 %f32 %f32 ');
fclose(fid);

figure('color','white');
hold on
grid off

for I = 2:9
    plot3(S{I}(1),S{I}(2),S{I}(3),'g*')
end

for I=2:5
    J = mod(I, 5) + 1;
    if (mod(I,5) == 0)
        J = 2;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=6:9
    J = mod(I, 9) + 1;
    if (mod(I, 9) == 0)
        J = 6;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=2:5 
    J = mod(I,5)  + 4; 
    if (mod(I,5) == 0)
        J = 9;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end
        

axis([0,300,0,300,0,300])
title('Quader t0')
xlabel('x')
ylabel('y')
zlabel('z')
view([30,30])
saveas(gcf,'QuaderT0.png')




%% Zeitpunkt t1
fid = fopen('QuaderT1.txt');
S = textscan(fid, '%s %f32 %f32 %f32 %f32 %f32 %f32 %f32 %f32 ');
fclose(fid);

figure('color','white');
hold on
grid off

for I = 2:9
    plot3(S{I}(1),S{I}(2),S{I}(3),'g*')
end

for I=2:5
    J = mod(I, 5) + 1;
    if (mod(I,5) == 0)
        J = 2;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=6:9
    J = mod(I, 9) + 1;
    if (mod(I, 9) == 0)
        J = 6;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=2:5 
    J = mod(I,5)  + 4; 
    if (mod(I,5) == 0)
        J = 9;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end
        

axis([0,300,0,300,0,300])
title('Quader t1')
xlabel('x')
ylabel('y')
zlabel('z')
view([30,30])
saveas(gcf,'QuaderT1.png')
%% Zeitpunkt t2
fid = fopen('QuaderT2.txt');
S = textscan(fid, '%s %f32 %f32 %f32 %f32 %f32 %f32 %f32 %f32 ');
fclose(fid);

figure('color','white');
hold on
grid off

for I = 2:9
    plot3(S{I}(1),S{I}(2),S{I}(3),'g*')
end

for I=2:5
    J = mod(I, 5) + 1;
    if (mod(I,5) == 0)
        J = 2;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=6:9
    J = mod(I, 9) + 1;
    if (mod(I, 9) == 0)
        J = 6;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end

for I=2:5 
    J = mod(I,5)  + 4; 
    if (mod(I,5) == 0)
        J = 9;
    end
    plot3([S{I}(1) S{J}(1)], [S{I}(2) S{J}(2)], [S{I}(3) S{J}(3)])
end
        

axis([0,300,0,300,0,300])
title('Quader t2')
xlabel('x')
ylabel('y')
zlabel('z')
view([30,30])
saveas(gcf,'QuaderT2.png')