Computational Physics - M Jensen Episode 1 Part 6 doc

Computational Physics - M. Jensen Episode 1 Part 6 doc

Computational Physics - M. Jensen Episode 1 Part 6 doc

... second-order polynomial . The first and second derivatives are given by (8.8) and (8.9) 90 CHAPTER 6. NON-LINEAR EQUATIONS AND ROOTS OF POLYNOMIALS 1. (6 .15 ) and 2. (6 . 16 ) 3. and a maximum number ... transferring 6. 6 Roots of polynomials in preparation 6. 6 .1 Polynomials division in preparation 6. 6.2 Root finding by Newton-Raphson’s method in preparation 6. 6.3 Root finding b...

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

20 250 0
Computational Physics - M. Jensen Episode 1 Part 3 doc

Computational Physics - M. Jensen Episode 1 Part 3 doc

... . sigma1 =SQRT( ( sumsq2 sum xbar ) / 1 26. ) ! second method to ev aluat e the st andard d e vi a tio n sumsq2 =0. DO i =1 , 12 7 sumsq2=sumsq2 +(x ( i ) xbar ) 2 ENDDO sigma2 =SQRT( sumsq2 / 1 2 6. ) WRITE( ... mechanical prob- lems from both atomic physics and nuclear physics. In atomic and nuclear physics the single- particle degrees are represented by quantum numbers such orb...

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

20 239 0
Computational Physics - M. Jensen Episode 1 Part 7 docx

Computational Physics - M. Jensen Episode 1 Part 7 docx

... [0 ,10 0] with using the Gauss-Legendre method. 1 1.305 3.334 2 6. 747 7.473 3 16 .030 10 .954 4 28.330 13 . 463 5 42.5 56 14 .7 76 6 57.444 14 .7 76 7 71. 67 0 13 . 463 8 83.970 10 .954 9 93.253 7.473 10 98 .69 5 ... different methods as functions of the num- ber of mesh points . Trapez Simpson Gauss-Legendre 10 1. 8 210 20 1. 214 025 0 .14 60 448 20 0. 912 67 8 0 .60 9...

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

20 202 0
Computational Physics - M. Jensen Episode 1 Part 1 ppsx

Computational Physics - M. Jensen Episode 1 Part 1 ppsx

... . 319 17 .1. 1 The histogram method . . . . . . . . . . . . . . . . . . . . . . . . . . . 319 17 .1. 2 Multi-histogram method . . . . . . . . . . . . . . . . . . . . . . . . . . 319 17 .2 Renormalization ... . . 215 12 .2.5 The helium atom . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 16 12 .2 .6 Program example for atomic systems . . . . . . . . . . . . . . . . . . . 2...

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

20 206 0
Computational Physics - M. Jensen Episode 1 Part 2 pot

Computational Physics - M. Jensen Episode 1 Part 2 pot

... 0.935 762 30E -1 3 -0 .3 066 811 1E-04 10 0 40.000000 0.42483543E -1 7 -0 . 3 16 57 319 E+ 01 127 50.000000 0 .19 287498E- 21 0 .11 072933E+05 15 5 60 .000000 0.87 56 510 8E- 26 -0 .33 5 16 811 E+09 18 2 70.000000 0.39754497E-30 -0 .3297 960 5E +14 ... 0.487 460 E-08 72 30.0 0.935 762 E -1 3 -0 .34 213 4E-04 10 0 40.0 0.424835E -1 7 -0 .2 210 33E+ 01 127 50.0 0 .19...

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

20 231 0
Computational Physics - M. Jensen Episode 1 Part 4 pdf

Computational Physics - M. Jensen Episode 1 Part 4 pdf

... 20.085537 20.250 467 20.085537 4.0 54 .64 366 4 54.59 860 5 54.59 815 5 54.59 815 1 54. 711 789 54.59 815 0 5.0 14 8.5 368 78 14 8. 414 3 96 14 8. 413 172 14 8. 413 1 61 150 .63 50 56 14 8. 413 159 Table 3 .1: Result for numerically ... are MyVector o bj ec t s ; want to s et a ( j ) = v ( i +1) ; 3.2. NUMERICAL DIFFERENTIATION 55 2. 010 00833 61 1 160 70 1. 000833 61 1 160 7230 2.00...

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

20 258 0
Computational Physics - M. Jensen Episode 1 Part 5 pptx

Computational Physics - M. Jensen Episode 1 Part 5 pptx

... 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 when ... a matrix. Again the elements start with zero, matr[0][0], matr[0] [1] , , matr[0][4], matr [1] [0], . This sequence of elements also shows how data are stored in memory. For example, the element matr...

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

20 324 0
Computational Physics - M. Jensen Episode 1 Part 8 pdf

Computational Physics - M. Jensen Episode 1 Part 8 pdf

... 999 10 30 993 0. 3-0 .4 939 960 10 23 937 0. 4-0 .5 10 38 10 01 1002 992 0. 5-0 .6 10 37 10 47 10 09 10 09 0. 6- 0 .7 10 05 989 10 03 989 0. 7-0 .8 9 86 962 985 954 0. 8-0 .9 10 00 10 27 10 09 10 23 0. 9 -1 .0 9 91 1 015 9 61 10 26 0.4997 ... samples. Table 9 .1: Results for as function of number of Monte Carlo samples . The exact answer is with 6 digits. 10...

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

20 242 0
Computational Physics - M. Jensen Episode 1 Part 9 pptx

Computational Physics - M. Jensen Episode 1 Part 9 pptx

... rewriting (9 .10 9) since (9 .11 0) Perform then a Monte Carlo sampling for (9 .11 1) with , 9.5. IMPROVED MONTE CARLO INTEGRATION 15 3 where is a random number in the interval [0 ,1] . The algorithm for the ... ran0(&idum ) ; } fx=brute_force_MC (x ) ; int_mc += fx ; sum_sigma += fx fx ; } int_mc = int_mc / ( ( double ) n ) ; sum_sigma = sum_sigma / ( ( double ) n ) ; va rianc e =sum_si...

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

20 313 0
Computational Physics - M. Jensen Episode 2 Part 6 pptx

Computational Physics - M. Jensen Episode 2 Part 6 pptx

... Four different methods have been used to solve the equations, Euler’s method from Eq. (14 .17 ), the half-step method, Euler-Richardson’s method in Eqs. (14 .32 )-( 14 .33) and finally the fourth-order Runge-Kutta ... its length m The inital velocity is and . 14 .6. PHYSICS EXAMPLES 2 71 with an angular velocity and acceleration given by (14 .65 ) and (14 .66 ) For small angles, we...

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

20 389 0
w