238
CHAPITRE 7. MODELE D'ISING
{
- 0
B···· 0.1
-
--- 0.01
-
0.001
0,5
'•,,
-0,5
2
4
6
8
T
FIGURE 7.5. Solutions de !'equation (7.14) pour T successivement B = 0 (trait epais), 0.1, 0.01 et 0,001. La ligne pointillee est la
solution instable.
champ magnetique) en-dessous de la temperature de Curie. Elle est egalement
verifiee en simulation numerique : pour T < T,, un tout petit champ B > O
conduit a un m > 0 et inversement pour B < 0. Nous donnons ci-apres un
programme en Fortran permettant de calculer la magnetisation m.
real "'8 J, q,m,mt9, T, B, ctr, Tmin, Tmax
integer i
B=l.d-6
q=4.d(9
J=l. d(S)
m=l.dt9
dT=l.d-3
Tmin=dT
Tmax=J'''q-dT
do T=Tmin,Tmax,dT
i=t9
1
continue
i=i+l
if (i.gt.1(9(9(9(9(9) stop 'diverge'