function[] = draw(NodalData,Connectivity)

num_elements = size(Connectivity,2);
num_nodes = size(NodalData,2);

for i=1:1:num_elements
    node1 = Connectivity(i).a;
    node2 = Connectivity(i).b;
    
    x = [NodalData(node1).x,NodalData(node2).x];
    y = [NodalData(node1).y,NodalData(node2).y];
    
    if (i==10 | (i==14) | (i==35) | (i==21) | (i==33) | (i==18) | (i==11) | (i==36) | (i==22) | (i==15))
        
        figure(1)
        plot(x,y,'k','LineWidth',4);
        
        %-----------------------------------
        % Use this for Theta = 15 degrees
        xlim([0,35])
        ylim([-10,25])
        
        % Use this for Theta = 30 degrees
%         xlim([0,35])
%         ylim([-6,29])
        %-----------------------------------
        hold on
        %     pause(0.2)
        axis square
    else
        figure(1)
        plot(x,y,'r','LineWidth',3);
        %-----------------------------------
        % Use this for Theta = 15 degrees
        xlim([0,35])
        ylim([-10,25])
        
        % Use this for Theta = 30 degrees
%         xlim([0,35])
%         ylim([-6,29])
        %-----------------------------------
        hold on
        %     pause(0.2)
        axis square
    end
end