Versão 1.0

Objetivo:
Calcula as componentes da magnetização (Mx, My, Mz) e o tensor
susceptibilidade Xvu=lim{Bv->0} dMu/dBv a partir do espectro de energia
fornecido.

Síntaxe:
[mag sus]=MX(ene), onde ene é o espectro de
energia do sistema de interesse. O campo magnético deve sempre
estar escrito como Bx, By e Bz.

Resultado no formato:
magnetizacao um vetor 3x1 com as três componentes da magnetização
(Mx;My;Mz).
susceptibilidade uma celula 3x3:
X={[Xxx] [Xxy] [Xxz]
   [Xyx] [Xyy] [Xyz]
   [Xzx] [Xzy] [Xzz]}

A variável t que aparecera é a temperatura vezes a constante de Boltzmann
t=k_B*T e tem unidade de energia.

%%%%%%%%%%%%%%%%%%
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 (agosto 2009) *ORIGINAL*