function [energy, stability] = energy_vs_stability( possible_energy, stability, index )
%ENERGY_VS_STABILITY Summary of this function goes here
%   Detailed explanation goes here

    n = size(stability,2);
    c = 1:n;
    
    for i=1:n
        fixed_stability(i,1) = stability(i);
    end
    
    for i = 1:n
        energy(i,1) = possible_energy(index(i));
    end
    
    figure(n+1)
    textCell = arrayfun(@(x,y) sprintf('(%3.2f, %3.2f)',x,y),energy,fixed_stability,'un',0);
    scatter(energy,fixed_stability,25,c,'filled')
    xlabel('Energy')
    ylabel('Stability')
    
%     for i = 1:n
%         text(energy(i), fixed_stability(i), textCell{i}, 'FontSize',8)
%     end
    
end