for i = 5;

    for j = 6:8

        colors = {'b' 'b' 'b' 'b' 'b' 'r' 'g' 'b'};

        plot(Bz,Em(i,:)-Em(j,:),colors{j}); hold on;

    end

end