%Team Blue Genes
%Jiayi Liang, Steven Hahn, William Schuppert, Runkun Jiang
%BioMOD competition 2011
%This file takes the fluctuations of the synthetic DNA structures and
%outputs a 3D color-coded representation of the figure.

function graphPointsNoText
% Specific DNA data of flexibility and coordinates
fileName = 'FILENAME.xlsx';

fprintf('Loading...\n');

intPointsMat = xlsread(fileName);
size_intPointsMat = size(intPointsMat);
rows_intPointsMat = size_intPointsMat(1, 1);

for r = 1 : rows_intPointsMat
    flux = intPointsMat(r, 5);
    if flux < 0.1
        X1(1, r) = intPointsMat(r, 2);
        Y1(1, r) = intPointsMat(r, 3);
        Z1(1, r) = intPointsMat(r, 4);
    elseif flux >= 0.1 && flux < 0.2
        X2(1, r) = intPointsMat(r, 2);
        Y2(1, r) = intPointsMat(r, 3);
        Z2(1, r) = intPointsMat(r, 4);
    elseif flux >= 0.2 && flux < 0.4
        X3(1, r) = intPointsMat(r, 2);
        Y3(1, r) = intPointsMat(r, 3);
        Z3(1, r) = intPointsMat(r, 4);
    elseif flux >= 0.4 && flux < 0.6
        X4(1, r) = intPointsMat(r, 2);
        Y4(1, r) = intPointsMat(r, 3);
        Z4(1, r) = intPointsMat(r, 4);
    else
        X5(1, r) = intPointsMat(r, 2);
        Y5(1, r) = intPointsMat(r, 3);
        Z5(1, r) = intPointsMat(r, 4);
    end
end

plot3(X1, Y1, Z1, 'b.', X2, Y2, Z2, 'c.', X3, Y3, Z3, 'g.', ...
    X4, Y4, Z4, 'y.', X5, Y5, Z5, 'r.'); 
xlabel('X');
ylabel('Y');
zlabel('Z');
legend('<0.1','0.1<F<0.2','0.2<F<0.4','0.4<F<0.6','>0.6','Location','EastOutside');

end