Computational Physics - M. Jensen Episode 2 Part 7 pot

... With respect to pointers this means that matr is pointer-to-a-pointer-to-an-integer which we can write matr. Furthermore matr is a-pointer-to-a-pointer of five integers. This interpretation is important ... 1 (by construction since the diagonals of equal 1) we can use the inverse of to obtain (5 .13 ) which yields the intermediate step (5 .14 ) and multiplying with on both sides we reob...

Ngày tải lên: 07/08/2014, 12:22

... the Gauss-Legendre method. 1 1.305 3.334 2 6 .74 7 7. 473 3 16 .030 10 .954 4 28.330 13 .463 5 42.556 14 .77 6 6 57. 444 14 .77 6 7 71 . 670 13 .463 8 83. 970 10 .954 9 93.253 7. 473 10 98.695 3.334 8.3.6 Applications ... Gauss-Legendre 10 1. 8 210 20 1. 214 025 0 .14 60448 20 0. 912 678 0.6098 97 0. 2 17 80 91 40 0. 478 456 0.333 71 4 0. 219 3834 10 0 0. 273 724 0.2 312...

Ngày tải lên: 07/08/2014, 12:22

... spin [ j ] spin [ j +1] ; } 10 .4. THE METROPOLIS ALGORITHM AND DETAILED BALANCE 18 1 Figure 10 .5: Probability distribution for one walker after 10 , 10 0 and 10 00 steps. 10 .5. PHYSICS PROJECT: SIMULATION ... statistical physics at a given temperature, we need a probability distribution (11 .2) 18 7 10 .3. MICROSCOPIC DERIVATION OF THE DIFFUSION EQUATION 17 7 with the obv...

Ngày tải lên: 07/08/2014, 12:22

... two-dimensional lattice we obtain the following partition function (11 .22 ) and mean energy (11 .23 ) 20 4 CHAPTER 12. QUANTUM MONTE CARLO METHODS which can be rewritten as ( 12. 5) In general, the integrals involved ... the one-dimensional Ising model with spins with free ends (FE) and periodic boundary conditions (PBC). Number spins up Degeneracy Energy (FE) Energy (PBC) Magnetization 2...

Ngày tải lên: 07/08/2014, 12:22

... 7.71 826 E-04 1 .20 000E+00 1.08667E+00 2. 233 89E-01 1.49462E- 03 1 .30 000E+00 1.17168E+00 4.78446E-01 2. 18 734 E- 03 1.40000E+00 1 .26 37 4E+00 8.55 524 E-01 2. 924 93E- 03 1.50000E+00 1 .38 897E+00 1 .30 720 E+00 3. 61553E- 03 ... of 2 was used in order to obtain an acceptance of . 5.00000E-01 2. 06479E+00 5.78 739 E+00 7.60749E- 03 6.00000E-01 1.50495E+00 2. 32 7 82E+00 4....

Ngày tải lên: 07/08/2014, 12:22

... 9.898985E-01 2. 949 052E+00 4. 86 622 3E+00 6.739916E+00 8.56 844 2E+00 100 9.9 748 93E-01 2. 98 744 2E+00 4. 96 727 7E+00 6.936913E+00 8.89 628 2E+00 20 0 9.993715E-01 2. 996864E+00 4. 991877E+00 6.9 843 35E+00 8.9 743 01E+00 40 0 ... the potential is symmetric with nm eV nm [eV] 86 42 0 -2 -4 -6 -8 0 -1 0 -2 0 -3 0 -4 0 -5 0 -6 0 Figure 12. 5: Plot of f...

Ngày tải lên: 07/08/2014, 12:22

... second-order equation goes typically like (14 .2) A well-known second-order equation is Newton’s second law (14.3) where is the force constant. ODE depend only on one variable, whereas 25 5 26 2 CHAPTER ... derivative in our Taylor expansion. We have then (14 .22 ) The second derivative can be rewritten as (14 .23 ) and we can rewrite Eq. (14.14) as (14 .24 ) 25 0 CHAPTER 13. EIGENSYS...

Ngày tải lên: 07/08/2014, 12:22

... errmax ; double yout [ 2] , y_h [ 2] , y_m [2] , y1 [ 2] , y2 [ 2 ] , delt a [ 2] , ys cal [ 2 ]; const double eps =1.0e 6; const double s a fet y =0.9; const double er rcon =6. 0e 4; const double ... phi and v / yout [0]= yin [ 0 ] +1.0 / 6 .0 ( k1 [0] +2 k2 [0] +2 k3 [0]+ k4 [ 0]) ; yout [1]= yin [ 1 ] +1.0 / 6 .0 ( k1 [1] +2 k2 [1] +2 k3 [1]+ k4 [ 1]) ; } 28 0 CHAPTER 1...

Ngày tải lên: 07/08/2014, 12:22

... following condition at (15 .22 ) We can modify this expression by normalizing the function . Then Eq. (15 .22 ) becomes at (15 .23 ) For an arbitary value of the eigenvalue Eq. (15 .22 ) will not be satisfied. ... simplified as (16 .26 ) and (16 . 27 ) The one-dimensional diffusion equation can then be rewritten in its discretized version as (16 .28 ) Defining results in the explicit scheme (16...

Ngày tải lên: 07/08/2014, 12:22

... Possible boundary and initial conditions with are (16. 58) 3 12 CHAPTER 16. PARTIAL DIFFERENTIAL EQUATIONS programs/chap16/program2.f90 ! Program to solve the 2 dim Laplace equation usi ng i t e r a t i ... as (16.47) or in matrix-vector form as (16. 48) where the vector is the same as defined in the implicit case while the matrix is (16.49) 16 .2. 5 Non-linear terms and implementation...

Ngày tải lên: 07/08/2014, 12:22

... Computational Physics (Wiley, 199 7). [7] E.W. Schmid, G. Spitz and W. Lösch, Theoretische Physik mit dem Personal Computer, (Springer, 198 7). [8] J .M. Thijssen, Computational Physics, (Springer, 199 9). [9] ... (Springer, 199 3). [4] S.E. Koonin and D. Meredith, Computational Physics, (Addison Wesley, 199 0). [5] F.S. Acton, Numerical Methods that work, (Harper and Row, 197 0)....

Ngày tải lên: 07/08/2014, 12:22