|
Identification Division.
Program-Id. Metropolis.
* from Stauffer/Stanley
* From Newton to Mandelbrot
* Springer 1996
Data Division.
Working-Storage Section.
01 dis pic S9(4)V9(4) occurs 1680.
77 w pic S9(4)V9(4) occurs 9.
77 vy pic S9(4)V9(4).
77 dt pic S9(4)V9(4).
Procedure Division.
move 40 to L
move 2.5 to t
compute L1=L+1
compute Lp=L*L+L
compute Lm=Lp+1
perform varying i from 1 until i>Lm
move 1 to is(i)
end-Perform
perform varying ie from 1 by 2 until ie>9
compute ex=function exp(-2*(ie-5)/t)
compute w(ie)=ex/(1+ex)
end-Perform
perform varying it from 1 until it>100
move 0 to m
perform varying i from L1 until i>Lp
compute ie=5+is(i)*(is(i-1)+is(i+1)+is(i-L)+is(i+L))
if function rnd(i)<w(ie) then move -is(i) to is(i)
add is(i) to m
end-Perform
display it, m
end-Perform
End-Program Metropolis.
¤ Dauer der Verarbeitung: 0.21 Sekunden
(vorverarbeitet)
¤
|
Haftungshinweis
Die Informationen auf dieser Webseite wurden
nach bestem Wissen sorgfältig zusammengestellt. Es wird jedoch weder Vollständigkeit, noch Richtigkeit,
noch Qualität der bereit gestellten Informationen zugesichert.
Bemerkung:
Die farbliche Syntaxdarstellung ist noch experimentell.
|
