Versão 1.0 Objetivo: Calcula as correções em energia de zero-, primeira- e segunda-ordem, utilizando teoria de pertubacao. Síntaxe: resultado=pertubacao(H0,W,bit), onde H0 é a Hamiltoniana não pertubada (diagonal) e W a pertubação (diagonal ou não, e hermitiana) e «bit» é o controle para agrupar (=1) ou não (=0) os autovalores degenerados. Exemplo: syms a b real;H0=[a 0 0;0 1 0;0 0 b];W=[1 1 0; 1 0 0; 0 0 -1]; resultado=pertubacao(H0,W) resultado = [ a, 1, 1/(a-1)] [ 1, 0, 1/(1-a)] [ b, -1, 0] resultado no formato [E_n(0) E_n(1) E_n(2)]. size(resultado)=[n 3] Observações: (1) Esta rotina serve tanto para H0 com autovalores degenerados e não-degenerados. (2) No caso de autovalores degenerados, estes devem estar agrupados. Exemplo: H0=[1 0 0 0;0 1 0 0;0 0 2 0; 0 0 0 2] %%%%%%%%%%%%%%%%%% Rotina otimizada para MATLAB 7.1 Rotina elaborada por Mario Reis (marior@if.uff.br). Favor comunicar possíveis erros e/ou 'bugs'. Histórico de alterações: Versão 1.0 (julho 2009) *ORIGINAL*