{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 } {CSTYLE "2D Input" 2 19 "" 0 1 255 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "2 D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 256 " " 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 1 16 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 } {CSTYLE "" -1 259 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 263 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 271 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 273 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 275 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 278 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 279 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 280 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 281 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 283 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 284 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 285 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 286 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 287 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 288 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 289 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 290 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 291 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 292 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 293 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 294 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 295 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 296 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 297 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 298 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 299 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 300 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 301 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" 2 302 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 303 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 304 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 305 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 306 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 307 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 308 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 309 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 310 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 311 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 312 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 313 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 314 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" 2 315 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 316 " " 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 317 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 318 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 319 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 320 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 1" 0 3 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 }1 0 0 0 8 4 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT 257 34 "Nuclear decays: Poisson s tatistics" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 388 "We consider sequences of nuclear decay data as collected by a Gei ger counter. The data are shown to obey Poisson statistics, which repr esents evidence for the fact that the events originate from independen t spontaneous decays. In the limit of a large sample of data (very hig h count rate) the distribution of events becomes Gaussian. A detailed \+ discussion can be found in John R. Taylor: " }{TEXT 256 33 "An Introdu ction to Error Analysis" }{TEXT -1 185 ", chapters 10,11,12 (Universit y Science Books 1982, 2nd ed. 1998). For readers who are less familiar with probability distributions we begin with a section on the binomia l distribution." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{SECT 1 {PARA 3 "" 0 "" {TEXT 309 21 "Binomial distribution" }}{EXCHG {PARA 0 "" 0 "" {TEXT 258 0 "" } {TEXT -1 1666 "We begin with the binomial distribution to illustrate t he statistics of a familiar experiment with several possible discrete \+ outcomes (throwing of dice, etc.), before proceeding with nuclear deca ys. Any physics measurement (or indeed any measurement at all) results in a distribution of results (think of the repeated measurement of th e period of a pendulum, or the measurement of the height of the female population in Canada, etc.). Usually we think that in physics measure ments there should be a unique answer: in the case of the pendulum pe riod we might think that only one outcome is possible in the measureme nt, while for the second example we clearly expect a distribution with a well-defined average and spread of data around the average. Neverth eless, both situations have something in common: due to measurement er rors we will obtain a distribution of results in the first case as wel l (we can think of using an electronic stopwatch, and realize that at \+ some level of precision we will start to investigate our reaction time ). Even in automated experiments (using a light beam/photodetector set -up) the numerical result will differ from measurement to measurement \+ when the readout is precise. The remaining question then is to underst and the properties of the distribution of values in the limit that the number of measurements taken is very large, or infinite. This is the \+ main reason for studying the mathematical properties of distributions. The distribution will be able to tell us the probability for achievin g a certain result in a particular measurement. The knowledge of the l atter is crucial while assessing the validity of a certain measurement ." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 342 "The distribution used most often in any experimental science is the Gauss ian (or normal) distribution. It applies when a measurement is subject to many sources of small and random errors. It is a limiting distribu tion in the sense that an experiment can be said to follow this distri bution only in the limit of a large number of measurements." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 390 "The binomial d istribution (and also the Poisson distribution discussed later) applie s to finite sequences of data points. Let us illustrate it by the exam ple of throwing three dice, and asking the question of finding the num ber of times that a particular face appears (e.g., let us call an ace \+ when the face with six points appears). Obviously, we can have for a g iven throw the outcome of " }{TEXT 19 10 "nu=0,1,2,3" }{TEXT -1 6 " ac es." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 331 "T his situation is easily figured out: the probability to throw one ace \+ with a single die is 1/6 as there are six faces, and we assume that th e die is perfect. Whether a die is perfect or not, can later be assess ed by observing whether its runs are following the correct distributio n. The same is true for a random number generator." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 82 "Clearly the probability t o throw three aces with three independet dies is given by" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "P3w3:=(1/6)^3; evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%P3w3G#\"\"\"\"$;#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+I'H'HY!#7" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 373 " This corresponds to a probability of about than 0.5 percent. Now we fi gure out the probability for two aces. It should be the product of thr owing two aces, and not throwing a third one, i.e., we use the fact th at not throwing an ace with one of the dies is given by the complement of 1/6 with 1. Naively, we might think that the answer would be given by the simple product:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 " wrongP2w3:=(1/6)^2*(1-1/6); evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*wrongP2w3G#\"\"&\"$;#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+:[ \"[J#!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 379 "This answer correspo nds to the case that the first two dice show aces, and the third one d oes not. The problem is that unlike in the first case, where the three dice were undistinguishable (they showed the same result), we can hav e three different configurations now: suppose we mark by A a face call ed ace, and by B the 'no ace' face; then we can have with equal probab ilities (" }{TEXT 267 1 "A" }{TEXT -1 2 ", " }{TEXT 266 1 "A" }{TEXT -1 2 ", " }{TEXT 265 1 "B" }{TEXT -1 4 "), (" }{TEXT 264 1 "A" }{TEXT -1 2 ", " }{TEXT 263 1 "B" }{TEXT -1 2 ", " }{TEXT 262 1 "A" }{TEXT -1 8 "), and (" }{TEXT 261 1 "B" }{TEXT -1 2 ", " }{TEXT 260 1 "A" } {TEXT -1 2 ", " }{TEXT 259 1 "A" }{TEXT -1 215 "). The probability for any one of these three configurations is the one calculated above by \+ the simple product. Why is this so? One should think about the likelih ood of the outcome of any such configuration, e.g., (" }{TEXT 270 1 "A " }{TEXT -1 2 ", " }{TEXT 269 1 "B" }{TEXT -1 2 ", " }{TEXT 268 1 "A" }{TEXT -1 511 "). Its probability is given by the product of probabili ties to dial an ace with the first die (1/6) times the probability of \+ not getting an ace with the second die (5/6), times the probability of getting one with the third die (1/6). This example shows that we need to be careful in probability theory, and that it is always easiest to find the probabilities for the most detailed (exclusive) information. The more inclusive information is obtained by summing over the possib ilities that make up the statement. " }}{PARA 0 "" 0 "" {TEXT -1 47 "T herefore, we have a probability of almost 7 %:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 34 "P2w3:=3*(1/6)^2*(1-1/6); evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%P2w3G#\"\"&\"#s" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+WWWWp!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 113 "N ow we can figure out the remaining probabilities: we consider all the \+ explicit possibilites for dialing one ace:" }}{PARA 0 "" 0 "" {TEXT -1 1 "(" }{TEXT 279 1 "A" }{TEXT -1 2 ", " }{TEXT 278 1 "B" }{TEXT -1 2 ", " }{TEXT 277 1 "B" }{TEXT -1 4 "), (" }{TEXT 276 1 "B" }{TEXT -1 2 ", " }{TEXT 275 1 "A" }{TEXT -1 2 ", " }{TEXT 274 1 "B" }{TEXT -1 4 "), (" }{TEXT 273 1 "B" }{TEXT -1 2 ", " }{TEXT 272 1 "B" }{TEXT -1 2 ", " }{TEXT 271 1 "A" }{TEXT -1 57 ") - this gives us the same multipl icity as with two aces:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 " P1w3:=3*(1-1/6)^2*(1/6); evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%%P1w3G#\"#D\"#s" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+AAAsM!#5" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 34 "Finally we have an easy one again: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "P0w3:=(5/6)^3; evalf(%) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%P0w3G#\"$D\"\"$;#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Pq.(y&!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 93 "We can check that the probabilities sum to unity (otherwise the se wouldn't be probabilities!)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "P0w3+P1w3+P2w3+P3w3;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\" " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 77 "A graphical representation of these results can be obtained by vertical bars:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 120 "plot([[0,P0w3*t,t=0..1],[1,P1w3*t,t=0..1],[2 ,P2w3*t,t=0..1],[3,P3w3*t,t=0..1]],color=[red,blue,green,black],axes=b oxed);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6)-%'CUR VESG6$7S7$\"\"!F(7$F($\"1i]k%Q49E\"!#<7$F($\"1,a]X\\&*eBF,7$F($\"1#p*R <(fKf$F,7$F($\"1!eQDEYd$[F,7$F($\"1gY,'[FB2'F,7$F($\"1@=ftlz=sF,7$F($ \"1fnp2*)*eS)F,7$F($\"1BbzR=gL'*F,7$F($\"1$3-8vOd3\"!#;7$F($\"1%\\&y#4 ;;@\"FE7$F($\"1O;G&o!\\A8FE7$F($\"1#[ce@4tW\"FE7$F($\"1&=5LESEd\"FE7$F ($\"1&=53&*>Mp\"FE7$F($\"1M5'zZ*4.=FE7$F($\"1$QuT/>N$>FE7$F($\"1`(>R_, S/#FE7$F($\"1&os) *)>tLFE7$F($\"15')z:fV%\\$FE7$F($\"1%>Q!y)Qu$FE 7$F($\"17!H%*o9!fQFE7$F($\"1'H@W7qR)RFE7$F($\"1FOb?$)>(4%FE7$F($\"1P?T Ym)4A%FE7$F($\"1%)\\*R\"=YPVFE7$F($\"1%>8\\#>@fWFE7$F($\"1Q*HA\"\\CyXF E7$F($\"1aua=9'Gq%FE7$F($\"1@V:D>)G#[FE7$F($\"15d5G-iX\\FE7$F($\"1H%Rf >Ut1&FE7$F($\"1vS#G/#>z^FE7$F($\"1CEg*e%Q2`FE7$F($\"1.ejJ:/AaFE7$F($\" 1]P4Y(*GWbFE7$F($\"1&p]e)oHhcFE7$F($\"1Pq.Pq.(y&FE-%'COLOURG6&%$RGBG$ \"*++++\"!\")F(F(-F$6$7S7$$\"\"\"F(F(7$F`u$\"1q.(yIc%ov!#=7$F`u$\"1TKI nHP:9F,7$F`u$\"1:)R/$e&f:#F,7$F`u$\"1[J_dxW,HF,7$F`u$\"1'z3;\\'RVOF,7$ F`u$\"1#4bT%zFJVF,7$F`u$\"1b!=YMRN/&F,7$F`u$\"19t(Q5h,y&F,7$F`u$\"1*\\ 7y]?W^'F,7$F`u$\"1kHrclppsF,7$F`u$\"1:)*o6T%\\$zF,7$F`u$\"1)))Q^HbQo)F ,7$F`u$\"166')z:%eV*F,7$F`u$\"16h[q>0;5FE7$F`u$\"1?mx'of=3\"FE7$F`u$\" 1IY]E96g6FE7$F`u$\"1_=N94SE7FE7$F`u$\"1]ilZ8^.8FE7$F`u$\"1T2*p9n^'z^\"FE7$F`u$\"1ja_*=qBf\"FE7$F`u$ \"1d')oqspg;FE7$F`u$\"1/7ResRMFE7$F`u$\"166O#R>R-#FE7$F`u$\"1m\"z%\\:m'4#FE 7$F`u$\"1;H#oKVq;#FE7$F`u$\"16')z!H!>XAFE7$F`u$\"12ul8)3aJ#FE7$F`u$\"1 yFlu?Q!R#FE7$F`u$\"1w@L#**=$eCFE7$F`u$\"1As%y)>fKDFE7$F`u$\"1\"*pR)3xC g#FE7$F`u$\"1;z%\\:Fbn#FE7$F`u$\"1jzLZp%pu#FE7$F`u$\"1s%G6&or@GFE7$F`u $\"1$f#4b\"HP*GFE7$F`u$\"1EM'o8st'HFE7$F`u$\"1eOc<`SSIFE7$F`u$\"1WWpD_ ^2JFE7$F`u$\"1u:w`2V%=$FE7$F`u$\"1\"[\")*=\\A`KFE7$F`u$\"1]ilZQdELFE7$ F`u$\"1J%Feu7$ Fi^l$\"1'HY]^&*G!eFeu7$Fi^l$\"1#fP;#F,7$Fi^l$\"1f#4I&GA?BF, 7$Fi^l$\"1/PqG=!GX#F,7$Fi^l$\"1+DJ&pAqg#F,7$Fi^l$\"1\"[\")RHMNu#F,7$Fi ^l$\"1Mek`\"4L*GF,7$Fi^l$\"1y_rQ-$f.$F,7$Fi^l$\"1E40z.u%=$F,7$Fi^l$\"1 :tPTXR@LF,7$Fi^l$\"12Cy;XzoMF,7$Fi^l$\"12*>_m,>i$F,7$Fi^l$\"1L3-))==bP F,7$Fi^l$\"1_o(*z\"G\"**QF,7$Fi^l$\"1AAs%yQy/%F,7$Fi^l$\"1L$e*)4BL>%F, 7$Fi^l$\"1Lek`m3MVF,7$Fi^l$\"1Asf\"e!Q!\\%F,7$Fi^l$\"19[JFw\"3j%F,7$Fi ^l$\"1cbI\\Tw!y%F,7$Fi^l$\"1_Vm%)zj;\\F,7$Fi^l$\"1XWpvR=l]F,7$Fi^l$\"1 \")RzwT&\\?&F,7$Fi^l$\"1Le*)4V0^`F,7$Fi^l$\"1Efn%*Q*Q\\&F,7$Fi^l$\"1Wp D-PVVcF,7$Fi^l$\"1&=&=5$euy&F,7$Fi^l$\"1_ostUuMfF,7$Fi^l$\"1:t7N1\"33' F,7$Fi^l$\"1*)))Q^/.:iF,7$Fi^l$\"1[J_2:')ojF,7$Fi^l$\"1jH'z$)\\k]'F,7$ Fi^l$\"1+DJ&pZJl'F,7$Fi^l$\"1L3-jib$z'F,7$Fi^l$\"1XWWWWWWpF,-Fft6&FhtF (FitF(-F$6$7S7$$\"\"$F(F(7$Fahl$\"1\\gr2v745!#>7$Fahl$\"1@VScR;()=Ffhl 7$Fahl$\"1`(>Rx2Y(GFfhl7$Fahl$\"1k3.5qfoQFfhl7$Fahl$\"1F<\")))>'y&[Ffh l7$Fahl$\"1cM()es.vdFfhl7$Fahl$\"11u:E\">Zs'Ffhl7$Fahl$\"1>k$=Z\")oq(F fhl7$Fahl$\"1lmT5S*eo)Ffhl7$Fahl$\"1]RGU(GHp*Ffhl7$Fahl$\"14`A[D*z0\"F eu7$Fahl$\"1&=&ost%y:\"Feu7$Fahl$\"1[\"[1@7\"e7Feu7$Fahl$\"1[\"[1'fta8 Feu7$Fahl$\"1F)oBezCW\"Feu7$Fahl$\"11&R`B:oa\"Feu7$Fahl$\"1-e8>7?N;Feu 7$Fahl$\"1L$3-8:!QFeu7$Fah l$\"1_oZ#\\`R-#Feu7$Fahl$\"1<1q_-;B@Feu7$Fahl$\"1w[eFIE9AFeu7$Fahl$\"1 Q\\&ynHDJ#Feu7$Fahl$\"1/m956g9CFeu7$Fahl$\"1))Q,#faM]#Feu7$Fahl$\"1M7l 'y=%*f#Feu7$Fahl$\"1:[\")*=f&)p#Feu7$Fahl$\"1)))QEt[bz#Feu7$Fahl$\"1b0 V-6R*)GFeu7$Fahl$\"1[\")Raqe$*HFeu7$Fahl$\"14Ka^<@(3$Feu7$Fahl$\"1Pq`* 4wr=$Feu7$Fahl$\"1,HWc'exF$Feu7$Fahl$\"1H'HrJ*ywLFeu7$Fahl$\"1()f>^%p* pMFeu7$Fahl$\"1b0$*R&ptc$Feu7$Fahl$\"1]RyHffiOFeu7$Fahl$\"1jz$[8*GiPFe u7$Fahl$\"1dM7SbIeQFeu7$Fahl$\"1oX[#='\\cRFeu7$Fahl$\"1V:vcP(Q0%Feu7$F ahl$\"1f#fUj`L9%Feu7$Fahl$\"1)4#orw!fC%Feu7$Fahl$\"1U'3`ALwL%Feu7$Fahl $\"1+](ozJaV%Feu7$Fahl$\"1c0o3v.HXFeu7$Fahl$\"1H'H'H'H'HYFeu-Fft6&FhtF (F(F(-%+AXESLABELSG6$%!GFiam-%*AXESSTYLEG6#%$BOXG-%%VIEWG6$%(DEFAULTGF abm" 1 2 0 1 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT 280 12 "Exercise B1:" }}{PARA 0 "" 0 "" {TEXT -1 90 "Repeat the above steps with four dice, i.e., find the probabili tes for 0, 1, 2, 3, 4 aces." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 88 "The general form of the binomial d istribution is specified by the following parameters: " }}{PARA 0 "" 0 "" {TEXT -1 36 "- the number of independent trials: " }{TEXT 282 1 " n" }{TEXT -1 11 " (throwing " }{TEXT 281 1 "n" }{TEXT -1 15 " dice, to ssing " }{TEXT 283 1 "n" }{TEXT -1 19 " coins, looking at " }{TEXT 284 1 "n" }{TEXT -1 37 " electrons in an atomic shell, etc.)." }} {PARA 0 "" 0 "" {TEXT -1 27 "- the number of successes: " }{TEXT 19 2 "nu" }{TEXT -1 77 " (getting an ace, getting a head on a coin; having \+ an electron ionized, etc.)" }}{PARA 0 "" 0 "" {TEXT -1 30 "- the proba bility of success: " }{TEXT 287 1 "p" }{TEXT -1 18 ", and of failure: \+ " }{TEXT 286 1 "q" }{TEXT -1 3 "=1-" }{TEXT 285 1 "p" }{TEXT -1 19 " i n any one trial (" }{TEXT 288 1 "p" }{TEXT -1 34 "=1/6 for getting an \+ ace on a die; " }{TEXT 289 1 "p" }{TEXT -1 130 "=1/2 for getting heads on a coin; probability of ionization of an electron in a given atomic shell by an X ray of a given energy)." }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 20 "The probability for " }{TEXT 19 2 "nu " }{TEXT -1 14 " successes in " }{TEXT 290 1 "n" }{TEXT -1 53 " trials has to be proportional to the nu-th power of " }{TEXT 292 1 "p" } {TEXT -1 29 ", and the (n-nu)-th power of " }{TEXT 291 1 "q" }{TEXT -1 136 ". The combinatorial factor that counts the number of possibili tes can be figured out by using a mathematical identity. We know that \+ 1 = " }{TEXT 294 1 "p" }{TEXT -1 1 "+" }{TEXT 293 1 "q" }{TEXT -1 39 " , and can raise this expression to the " }{TEXT 295 1 "n" }{TEXT -1 35 "th power. Examples are given below:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "expand(1^4=(p+q)^3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/\"\"\",**$)%\"pG\"\"$\"\"\"F$*&)F(\"\"#F*%\"qGF$F)*&F(F$)F.F-F*F)* $)F.F)F*F$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "expand(1^4=(p +q)^4);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/\"\"\",,*$)%\"pG\"\"%\"\" \"F$*&)F(\"\"$F*%\"qGF$F)*&)F(\"\"#F*)F.F1F*\"\"'*&F(F$)F.F-F*F)*$)F.F )F*F$" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 60 "We recognize the followi ng from our example discussed above:" }}{PARA 0 "" 0 "" {TEXT -1 34 "- the expansion of the expression " }{TEXT 19 7 "(p+q)^n" }{TEXT -1 116 " involves the required powers of single-event probabilities for t he desired success rates for all possible outcomes;" }}{PARA 0 "" 0 " " {TEXT -1 119 "- the sum of the terms on the RHS equals one, i.e., we can think of the identity as summing all possible probabilities;" }} {PARA 0 "" 0 "" {TEXT -1 6 "- for " }{TEXT 296 1 "n" }{TEXT -1 8 "=3 ( and " }{TEXT 297 1 "n" }{TEXT -1 130 "=4 in the exercise) the binomial coefficients count the number of configurations for a given scenario \+ as counted explicitly above." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 29 "Maple has a built-in function" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "binomial(4,2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 87 "but it is straightforward to define the function known from the mathematical ex pansion:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "bc:=(n,nu)->n!/ (nu!*(n-nu)!);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#bcGR6$%\"nG%#nuG6 \"6$%)operatorG%&arrowGF)*&-%*factorialG6#9$\"\"\"*&-F/6#9%\"\"\"-F/6# ,&F1\"\"\"F6!\"\"\"\"\"!\"\"F)F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "bc(4,2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"'" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 68 "Armed with this coefficient we can define the binomial distribution:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "bc(n,nu)*p^nu*(1-p)^(n-nu);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&*(-%*factorialG6#%\"nG\"\"\")%\"pG%#nuGF)),&F)F)F+!\" \",&F(F)F,F/F)\"\"\"*&-F&6#F,\"\"\"-F&6#F0\"\"\"!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "BD:=unapply(%,p,n,nu);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#BDGR6%%\"pG%\"nG%#nuG6\"6$%)operatorG%&arrow GF**&*(-%*factorialG6#9%\"\"\")9$9&F4),&F4F4F6!\"\",&F3F4F7F:F4\"\"\"* &-F16#F7\"\"\"-F16#F;\"\"\"!\"\"F*F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 55 "For a given case (given single-event probability value " }{TEXT 299 1 "p" }{TEXT -1 28 " and fixed number of trials " }{TEXT 298 1 "n" }{TEXT -1 37 ") we can graph the probabilities for " }{TEXT 19 2 "nu" }{TEXT -1 11 " successes:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "ps:=0.35;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#psG$\" #N!\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 176 "plot([[0,BD(ps,4 ,0)*t,t=0..1],[1,BD(ps,4,1)*t,t=0..1],[2,BD(ps,4,2)*t,t=0..1],[3,BD(ps ,4,3)*t,t=0..1],[4,BD(ps,4,4)*t,t=0..1]],color=[red,blue,green,brown,m agenta],axes=boxed);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6*-%'CURVESG6$7S7$\"\"!F(7$F($\"1UN'4DG4*Q!#=7$F($\"1:mW< TSwsF,7$F($\"1[u^!es$36!#<7$F($\"1NrC*[G;\\\"F37$F($\"1*4V&\\G1t=F37$F ($\"1XU`G8qEAF37$F($\"1`%*)**[tGf#F37$F($\"1?6=j#p:(HF37$F($\"1\"RDUb] !\\LF37$F($\"1@ITgkLPPF37$F($\"16:$z*)Q$zSF37$F($\"17._RENkWF37$F($\"1 PM_%GZ4&[F37$F($\"1(oH#fG]B_F37$F($\"1q[_7!>=c&F37$F($\"1r$GQ!yF37$F($\"1:m%pHVj=)F37$F($\"1u\\ma9hP&)F37$F($\"15S W<;];*)F37$F($\"1=6V6?15$*F37$F($\"1s'Q8\\dEl*F37$F($\"1q[`**pE-5!#;7$ F($\"1iSx6H\\S5Ffp7$F($\"1TTE4&*)y2\"Ffp7$F($\"13IBSE296Ffp7$F($\"1%)* =\"QyCa6Ffp7$F($\"1YOGSqM!>\"Ffp7$F($\"1W[\"Qk!*)G7Ffp7$F($\"1()*>t(o \"QE\"Ffp7$F($\"1\"GodZ+?I\"Ffp7$F($\"1;2>m#GzL\"Ffp7$F($\"19WIfJ[v8Ff p7$F($\"1_DtJ**>79Ffp7$F($\"1X9t*)*Q1X\"Ffp7$F($\"1#zu6Pgw[\"Ffp7$F($ \"18s\"33?b_\"Ffp7$F($\"1T)[HHmIc\"Ffp7$F($\"1i!*\\itc(f\"Ffp7$F($\"1= tk\"[4rj\"Ffp7$F($\"1U&)RCkZs;Ffp7$F($\"1$QK>'\\=5A]Dsc\"F37$Fau$\"1Ugd>SE( Q#F37$Fau$\"1#HKXGQF@$F37$Fau$\"1@0S9#*GMSF37$Fau$\"1/wA#4sfz%F37$Fau$ \"1PM#Qf]Ye&F37$Fau$\"1/,RGII+kF37$Fau$\"1(oa_lRL@(F37$Fau$\"1&e*)3wZ' \\!)F37$Fau$\"1$zWbRoiy)F37$Fau$\"1]P/i(Gbh*F37$Fau$\"1D\"G^T>[/\"Ffp7 $Fau$\"1DcHa@1D6Ffp7$Fau$\"15:&)[-$z>\"Ffp7$Fau$\"1e9Bhud%G\"Ffp7$Fau$ \"1L3l:!zzN\"Ffp7$Fau$\"1%fBJZiLW\"Ffp7$Fau$\"1$3d1IT*=:Ffp7$Fau$\"1%4 hIDj=g\"Ffp7$Fau$\"1\"G$3hX#3o\"Ffp7$Fau$\"1r#>CS7Kw\"Ffp7$Fau$\"1-jYf ,()Q=Ffp7$Fau$\"1z%\\9Fx/#>Ffp7$Fau$\"15!Rb7W_+#Ffp7$Fau$\"119@@R.z?Ff p7$Fau$\"1/^2A*G(e@Ffp7$Fau$\"1](Ghlh5C#Ffp7$Fau$\"1(o*zt\"3;K#Ffp7$Fa u$\"1%4r72T&*R#Ffp7$Fau$\"17yrmA2'[#Ffp7$Fau$\"1$3dWRCQc#Ffp7$Fau$\"1D \"ymQTok#Ffp7$Fau$\"16:Q?r1AFFfp7$Fau$\"1v$Rqr3V!GFfp7$Fau$\"1FpD];p\" )GFfp7$Fau$\"1ozM7\"zD'HFfp7$Fau$\"1eRU\"Rh;/$Ffp7$Fau$\"1WBl3KXCJFfp7 $Fau$\"1nToo:>/KFfp7$Fau$\"1HKoebt&G$Ffp7$Fau$\"1FWFBVgmLFfp7$Fau$\"1] (G#>V\"4W$Ffp7$Fau$\"1Y6Zg>3ENFfp7$Fau$\"1m\"*39pD-OFfp7$Fau$\"1%fBEIv Mo$Ffp7$Fau$\"12Roj;@hPFfp7$Fau$\"1+++++vWQFfp-Fgt6&FitF(F(Fjt-F$6$7S7 $$\"\"#F(F(7$Fi^l$\"1]7.iP$)onF,7$Fi^l$\"1>F37$Fi^l$\"18`/**R!\\f#F37$Fi^l$\"1VtZlVYeKF37$Fi^l$\"1>#*zm(pO(Q F37$Fi^l$\"1p/yD$z1^%F37$Fi^l$\"1p/$H_v%p^F37$Fi^l$\"1%4@B&*eh#eF37$Fi ^l$\"1EJ5w&Q;]'F37$Fi^l$\"1j:S'R,m4(F37$Fi^l$\"1vo\\:')QmxF37$Fi^l$\"1 DcEXg#*Q%)F37$Fi^l$\"1DJ+x//(3*F37$Fi^l$\"1pH\\r/fv'*F37$Fi^l$\"1j:2ML aP5Ffp7$Fi^l$\"1]([\\?Ho4\"Ffp7$Fi^l$\"1sOvmFzl6Ffp7$Fi^l$\"1voP/s$oA \"Ffp7$Fi^l$\"1ACZFE\"QH\"Ffp7$Fi^l$\"1msco!*ed8Ffp7$Fi^l$\"1>Z:^:Ffp7$Fi^l$\"1ZICb-i>;Ffp 7$Fi^l$\"1GQf0(>#z;Ffp7$Fi^l$\"1pHvW()eV2E)35=Ffp7$Fi ^l$\"1$f`M1X^(=Ffp7$Fi^l$\"1AuDl34Q>Ffp7$Fi^l$\"1c^>Y9)z+#Ffp7$Fi^l$\" 1vonT7yq?Ffp7$Fi^l$\"1j:qFM$y8#Ffp7$Fi^l$\"1(HaTn#f)>#Ffp7$Fi^l$\"1)=Z \"z&=]E#Ffp7$Fi^l$\"1sO0z,_FBFfp7$Fi^l$\"1fL(*fB&GR#Ffp7$Fi^l$\"17yX3l scCFfp7$Fi^l$\"1ZIPAofBDFfp7$Fi^l$\"1]()R#)3+)e#Ffp7$Fi^l$\"1J&feLjQl# Ffp7$Fi^l$\"1AC&\\T!=>FFfp7$Fi^l$\"1v$pq-+#zFFfp7$Fi^l$\"1c^wv#*)z%GFf p7$Fi^l$\"1]7%)*>:&4HFfp7$Fi^l$\"1sO]nY6vHFfp7$Fi^l$\"1y+OG#=#>;F37$Fahl$ \"1/w2i\"4d&=F37$Fahl$\"1(=#>v;W\"4#F37$Fahl$\"1%eR,a@RL#F37$Fahl$\"1# zpr$z\\ZDF37$Fahl$\"1]PzuV$zy#F37$Fahl$\"1]7`n!e$HIF37$Fahl$\"1]i?IX,i KF37$Fahl$\"10E?()))GtMF37$Fahl$\"1$ekD`9Xs$F37$Fahl$\"1L$3vgNt$RF37$F ahl$\"1PMU6r*[=%F37$Fahl$\"1L32!*4./WF37$Fahl$\"1Q%)z&ebWk%F37$Fahl$\" 18`RcmRt[F37$Fahl$\"14x1YEF7^F37$Fahl$\"1@0gt`jJ`F37$Fahl$\"1#z>,rZ#ob F37$Fahl$\"1/w2H9-9eF37$Fahl$\"1il#Rl.fiF37$Fahl$ \"1+vyIFv(\\'F37$Fahl$\"1v=Pr**GJnF37$Fahl$\"1Q%)*)f#\\s&pF37$Fahl$\"1 DJ![_Q\"3sF37$Fahl$\"1L32)otNV(F37$Fahl$\"1]7.VQFuwF37$Fahl$\"1/E]eRQ# *yF37$Fahl$\"1^P9Df$38)F37$Fahl$\"1qnvU1?b$)F37$Fahl$\"1(=))F37$Fahl$\"1P4@;b1f!*F37$Fahl$\"1n;%ydf-H*F37$Fahl$ \"1#Hdz*))oE&*F37$Fahl$\"1r<$RYg6w*F37$Fahl$\"1+vy\"[:m(**F37$Fahl$\"1 YOjC^NA5Ffp7$Fahl$\"1n\"*y`8WW5Ffp7$Fahl$\"1WBk!))*)z1\"Ffp7$Fahl$\"1c ^cD)G04\"Ffp7$Fahl$\"1+++++v96Ffp-Fgt6&Fit$\")#)eqkF\\u$\"))eqk\"F\\uF gam-F$6$7S7$$\"\"%F(F(7$F]bm$\"17$F]bm$\"1Y64!*)ep6'Fbbm7$ F]bm$\"1\"z%*z!*4wJ*Fbbm7$F]bm$\"1aj%=pZRD\"F,7$F]bm$\"1*[Vd4.Yd\"F,7$ F]bm$\"1[\\!>q\"*=(=F,7$F]bm$\"1Jq&R$orz@F,7$F]bm$\"1)pte-q!)\\#F,7$F] bm$\"11k1bAS:GF,7$F]bm$\"13-)yX<=9$F,7$F]bm$\"1/,>YTKHMF,7$F]bm$\"1DJX \"e))Hv$F,7$F]bm$\"1vV)R,#)z2%F,7$F]bm$\"1vo/1M<\"R%F,7$F]bm$\"1)>6V> \"evYF,7$F]bm$\"13x1j&pP,&F,7$F]bm$\"1Lekr%f-I&F,7$F]bm$\"1\"y+BM:Nj&F ,7$F]bm$\"1$ekeD.&GfF,7$F]bm$\"1\"G8,v^@D'F,7$F]bm$\"1W)RGlT.c'F,7$F]b m$\"1Y64&[0>)oF,7$F]bm$\"1SAYk9?xrF,7$F]bm$\"1/^Jr!=d\\(F,7$F]bm$\"1)p t$3tcEyF,7$F]bm$\"1=U#=asX6)F,7$F]bm$\"1y>6%3EcU)F,7$F]bm$\"1\\i!*H[(p u)F,7$F]bm$\"1h!RX)=Nh!*F,7$F]bm$\"1\"Gj)\\y_l$*F,7$F]bm$\"1OMx)RjKq*F ,7$F]bm$\"1ek3\"Qt1+\"F37$F]bm$\"1PM-D`2L5F37$F]bm$\"1?h5-jVi5F37$F]bm $\"18yI0c`%4\"F37$F]bm$\"16M1[&QZ7\"F37$F]bm$\"1\"RDeK4j:\"F37$F]bm$\" 1@0qJb<(=\"F37$F]bm$\"1.2s)e*[>7F37$F]bm$\"1!o9F37$F]bm $\"1++++]i+:F3-Fgt6&FitFjtF(Fjt-%+AXESLABELSG6$%!GFe[n-%*AXESSTYLEG6#% $BOXG-%%VIEWG6$%(DEFAULTGF]\\n" 1 2 0 1 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT 300 12 "Exercis e B2:" }}{PARA 0 "" 0 "" {TEXT -1 64 "Explore the possibilites of the \+ binomial distribution for fixed " }{TEXT 301 1 "n" }{TEXT -1 44 "=4 by changing the single-event probability " }{TEXT 19 2 "ps" }{TEXT -1 40 " in the above example. What happens for " }{TEXT 19 6 "ps=0.5" } {TEXT -1 19 "? Explore also for " }{TEXT 302 1 "n" }{TEXT -1 55 "=5. G raph the probabilities for all possible successes " }{TEXT 19 7 "nu=0. .n" }{TEXT -1 47 " as a function of the single-event probability." }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 77 "T here are further properties that one should explore. The general expre ssion " }{TEXT 19 13 "BD(ps, n, nu)" }{TEXT -1 27 " gives the probabil ity for " }{TEXT 19 2 "nu" }{TEXT -1 14 " successes in " }{TEXT 19 1 " n" }{TEXT -1 51 " trials for given single-trial success probability " }{TEXT 19 2 "ps" }{TEXT -1 84 ". One can ask what the average number f or success is. This is obtained by averaging " }{TEXT 19 2 "nu" } {TEXT -1 37 " probabilistically, i.e., by summing " }{TEXT 19 2 "nu" } {TEXT -1 46 " with the probability as a weight coefficient." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 78 "We try a few ex amples (to anticipate the result of a mathematical derivation)." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "N:=5; sum(nu*BD(p,N,nu),nu=0 ..N); simplify(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\"\"&" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#,,*&%\"pG\"\"\"),&F&F&F%!\"\"\"\"%\"\" \"\"\"&*&)F%\"\"#F+)F(\"\"$F+\"#?*&)F%F1F+)F(F/F+\"#I*&)F%F*F+F(F&F2*$ )F%F,F+F," }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"pG\"\"&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 51 "The average number of successes is simply given by " }{TEXT 303 1 "n" }{TEXT -1 346 " times the single-event pr obability. Sometimes this quantity is all that one is interested in. O n the other hand one can argue that this is all that is required from \+ a measurement, and that the individual success probabilities can be re constructed via the deduced single-trial probability (assuming that th e experiment obeys binomial statistics)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 51 "The standard deviation for the num ber of successes " }{TEXT 19 2 "nu" }{TEXT -1 123 " can also be calcul ated. It is defined as follows: sigma-squared is given as the average \+ of the square of the deviation of " }{TEXT 19 2 "nu" }{TEXT -1 18 " fr om the average." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "nubar:=N *p;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&nubarG,$%\"pG\"\"&" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sigsq:=simplify(sum((nu-nuba r)^2*BD(p,N,nu),nu=0..N));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&sigsq G,&%\"pG\"\"&*$)F&\"\"#\"\"\"!\"&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 52 "We can generalize this finding (try other values of " }{TEXT 304 1 "N" }{TEXT -1 89 "). It is also possible to show that this result is identical to the following difference:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "simplify(sum(nu^2*BD(p,N,nu),nu=0..N))-nubar^2;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#,&%\"pG\"\"&*$)F$\"\"#\"\"\"!\"&" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 95 "The quantity sigma (or its square) measures the width of the distribution of probabilities for " } {XPPEDIT 18 0 "nu;" "6#%#nuG" }{TEXT -1 275 " successes. The two resul ts for the average and the width of the binomial distribution can be u sed to observe which parameters to choose in a Gaussian approximation \+ to the binomial distribution. The latter becomes a valid approximation to the binomial distribution for large " }{TEXT 305 1 "n" }{TEXT -1 2 ")." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "N:=7; ps:=0.4;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\"\"(" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#psG$\"\"%!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "P1:=plot([seq([nu,BD(ps,N,nu)*t,t=0..1],nu= 0..N)],color=blue):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "nuba r:=ps*N; sigma:=sqrt(N*ps*(1-ps));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%&nubarG$\"#G!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&sigmaG$\"+S \"[hH\"!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "P2:=plot(exp (-(nu-nubar)^2/(2*sigma^2))/(sqrt(2*Pi)*sigma),nu=-0.1..8,color=red): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "display(P1,P2,axes=boxe d);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6.-%'CURVES G6$7S7$\"\"!F(7$F($\"1****z[\"3=5'!#>7$F($\"1+!)Gze4T6!#=7$F($\"1+!))e Rl\"QRBF07$F($\"1+?**zLOPHF07$F($\"1+g,$>W>\\$F07$F ($\"1**zK-0=mSF07$F($\"1**>([^c+m%F07$F($\"1+?^fz-__F07$F($\"1,S'GNV4' eF07$F($\"1+!oFOwsR'F07$F($\"1**>z,.1,qF07$F($\"1++oNMK2wF07$F($\"1**z +..d\">)F07$F($\"1+?RIC7A()F07$F($\"1+?&fA+IN*F07$F($\"1+!o6ZOu))*F07$ F($\"1+%QvG640\"!#<7$F($\"1+G@$=Sf5\"F[o7$F($\"1+%yLg:j;\"F[o7$F($\"1+ #>4G2QA\"F[o7$F($\"1++;:Qz$G\"F[o7$F($\"1+COz,))Q8F[o7$F($\"1+S1l#)H)R \"F[o7$F($\"1+#fp&p,g9F[o7$F($\"1+sofKu8:F[o7$F($\"1+O^Q\"pF[o7$F($\"1++1=`$>)>F[o7$F( $\"1+[kr`\"=/#F[o7$F($\"1+k]Ny:)4#F[o7$F($\"1+[eW>0d@F[o7$F($\"1+?2,M#R#F[o7 $F($\"1+_dSCA^CF[o7$F($\"1+Oz+vK0DF[o7$F($\"1+ob0!Qtc#F[o7$F($\"1+'\\- %4!Gi#F[o7$F($\"1+%em*f$>o#F[o7$F($\"1+[5\"zN&QFF[o7$F($\"1+++++O*z#F[ o-%'COLOURG6&%$RGBGF(F($\"*++++\"!\")-F$6$7S7$$\"\"\"F(F(7$Fau$\"1++Wp /^ZGF07$Fau$\"1+SM.T6D`F07$Fau$\"1+S99&Q96)F07$Fau$\"1+3=\\Ki\"4\"F[o7 $Fau$\"1+'Htdp2P\"F[o7$Fau$\"1+3ucRdH;F[o7$Fau$\"1+k3,4b(*=F[o7$Fau$\" 1+OFSIpu@F[o7$Fau$\"1+c5\"QY4X#F[o7$Fau$\"1+K+)*o5NFF[o7$Fau$\"1+%efj& R&)HF[o7$Fau$\"1+'pTZhrE$F[o7$Fau$\"1+S)*pU3]NF[o7$Fau$\"1+/23GtAQF[o7 $Fau$\"1+'\\3!QKqSF[o7$Fau$\"1*fx(QMtkVF[o7$Fau$\"1+%y)>q89YF[o7$Fau$ \"1*>z^n_U!\\F[o7$Fau$\"1*RE$)=a5;&F[o7$Fau$\"1+#*4\\h!GW&F[o7$Fau$\"1 +'*GW156dF[o7$Fau$\"1++3/y.\"*fF[o7$Fau$\"1+7p.v5[iF[o7$Fau$\"1+?j.>RD lF[o7$Fau$\"1)f4eY7M\"oF[o7$Fau$\"1+O(=@NT1(F[o7$Fau$\"1*z'RYE#\\L(F[o 7$Fau$\"1*zGGcsYh(F[o7$Fau$\"1)f,'yKN))yF[o7$Fau$\"1)*>^VM:`\")F[o7$Fa u$\"1*z_%z)orW)F[o7$Fau$\"1*R5b#[N6()F[o7$Fau$\"1**zG7(HM**)F[o7$Fau$ \"1***zU[J!\\#*F[o7$Fau$\"1+CMMjBB.9z*F[o7$Fau$\"1S i!3CCm+\"!#;7$Fau$\"1+O$Qz%\\L5Fb\\l7$Fau$\"1++isdih5Fb\\l7$Fau$\"1Suj [#>()3\"Fb\\l7$Fau$\"1?F1TiU;6Fb\\l7$Fau$\"1g\"Fb\\l7$Fau$\"1![;@xtRA\"Fb\\l7$Fau$\"1 ?RxJ,d^7Fb\\l7$Fau$\"1SA=pL)zF\"Fb\\l7$Fau$\"1++++!ojI\"Fb\\lFft-F$6$7 S7$$\"\"#F(F(7$Fh^l$\"1***z)Q4-&p&F07$Fh^l$\"1+)o1#G-l5F[o7$Fh^l$\"1+) GGq(GA;F[o7$Fh^l$\"1+;O)\\YK=#F[o7$Fh^l$\"1+#fY:R:u#F[o7$Fh^l$\"1+;[8z 9fKF[o7$Fh^l$\"1+G<-=5&z$F[o7$Fh^l$\"1+sa!3'Q\\VF[o7$Fh^l$\"1*>6Aw#*=! \\F[o7$Fh^l$\"1,k+'z8-Z&F[o7$Fh^l$\"1+o\">F\"zqfF[o7$Fh^l$\"1+#R$[HKMl F[o7$Fh^l$\"1**z'*R&o,5(F[o7$Fh^l$\"1*zShhlak(F[o7$Fh^l$\"1+#*p,wkS\") F[o7$Fh^l$\"1*>bv(oYH()F[o7$Fh^l$\"1*zc(RSFG#*F[o7$Fh^l$\"1*Re.N0&3)*F [o7$Fh^l$\"1!Glw$3@K5Fb\\l7$Fh^l$\"1S)>)H7c)3\"Fb\\l7$Fh^l$\"1?z&)G,AU 6Fb\\l7$Fh^l$\"1+g\"3c2#)>\"Fb\\l7$Fh^l$\"1S#Q2]@'\\7Fb\\l7$Fh^l$\"1+k s!Qy]I\"Fb\\l7$Fh^l$\"1?>;$\\#oi8Fb\\l7$Fh^l$\"1?ZPUq#GT\"Fb\\l7$Fh^l$ \"1g$z#HX)pY\"Fb\\l7$Fh^l$\"1gdc7X$H_\"Fb\\l7$Fh^l$\"1?.sb1nx:Fb\\l7$F h^l$\"1+Cq)oI1j\"Fb\\l7$Fh^l$\"1g0*exL%*o\"Fb\\l7$Fh^l$\"1!3-^'4FUFb\\l7$Fh^l$\"1SEZY1Ge>Fb\\l7$Fh^l$\"1![7;[[K,#Fb\\l7$Fh^l$\"1+sm(e* )p1#Fb\\l7$Fh^l$\"1++CX:DB@Fb\\l7$Fh^l$\"1!)[F(\\Qu<#Fb\\l7$Fh^l$\"1Sa 7#[_GB#Fb\\l7$Fh^l$\"1?N?6w!yG#Fb\\l7$Fh^l$\"1gtSncIQBFb\\l7$Fh^l$\"1! o>&Q@='R#Fb\\l7$Fh^l$\"1gHBWv%zW#Fb\\l7$Fh^l$\"1Syaj-9.DFb\\l7$Fh^l$\" 1![k$Qn'fb#Fb\\l7$Fh^l$\"1++++gt7EFb\\lFft-F$6$7S7$$\"\"$F(F(7$F^hl$\" 1++?V5!yK'F07$F^hl$\"1+?j*oeL=\"F[o7$F^hl$\"1+?.#*=a-=F[o7$F^hl$\"1+Si U%HeU#F[o7$F^hl$\"1+!)G;Y:YIF[o7$F^hl$\"1+SUfaF@OF[o7$F^hl$\"1+?>!yzn@ %F[o7$F^hl$\"1+!31?^E$[F[o7$F^hl$\"1**zc8vaYaF[o7$F^hl$\"1,g*)G`,ygF[o 7$F^hl$\"1+?NCD@MmF[o7$F^hl$\"1**z[J)e.E(F[o7$F^hl$\"1++_bh2*)yF[o7$F^ hl$\"1+?r&zh\\\\)F[o7$F^hl$\"1**z)=+k^/*F[o7$F^hl$\"1+!G#yj`-\"Fb\\l7$F^hl$\"1+w\"y#R$)*3\"Fb\\l7$F^hl$\"1+#R2$4!p9\" Fb\\l7$F^hl$\"1+wdxC^47Fb\\l7$F^hl$\"1+)3lZL\"p7Fb\\l7$F^hl$\"1++CMFb\\l7$F^hl$\"1+S1O*R&)*>Fb\\l 7$F^hl$\"1++%)H.Mb?Fb\\l7$F^hl$\"1+s=j\"Qu6#Fb\\l7$F^hl$\"1+'p\\Qne<#F b\\l7$F^hl$\"1+sMdF%pB#Fb\\l7$F^hl$\"1+!3u4bmH#Fb\\l7$F^hl$\"1++g$Qo\" fBFb\\l7$F^hl$\"1+K33hP>CFb\\l7$F^hl$\"1+;O-s%4[#Fb\\l7$F^hl$\"1+G*oX3 ?a#Fb\\l7$F^hl$\"1+/B3u6)f#Fb\\l7$F^hl$\"1+_N)fCCm#Fb\\l7$F^hl$\"1+W\" erT*>FFb\\l7$F^hl$\"1+w\\fpE\"y#Fb\\l7$F^hl$\"1+siUI'*RGFb\\l7$F^hl$\" 1+++++/.HFb\\lFft-F$6$7S7$$\"\"%F(F(7$Fdam$\"1++!)GS`=UF07$Fdam$\"1,+) 3$z0*)yF07$Fdam$\"1+!)o%f%p,7F[o7$Fdam$\"1+gTG'>sh\"F[o7$Fdam$\"1+?>W( p2.#F[o7$Fdam$\"1+ghRO=9CF[o7$Fdam$\"1+!G,_'=6GF[o7$Fdam$\"1+?2nuw@KF[ o7$Fdam$\"1+?rv;.JOF[o7$Fdam$\"1+SE>-,_SF[o7$Fdam$\"1,!o&\\$3GU%F[o7$F dam$\"1+?*4AR-%[F[o7$Fdam$\"1++o.TQf_F[o7$Fdam$\"1+!3Q'yIjcF[o7$Fdam$ \"1+?fM$4,.'F[o7$Fdam$\"1+?:hFb\\lFft-F$6$7S7$$\"\"&F(F(7$Fjj m$\"1++_6OT(o\"F07$Fjjm$\"1+?NsJibJF07$Fjjm$\"1+?vy$yn![F07$Fjjm$\"1** Rm8&y)okF07$Fjjm$\"1**zww*yI7)F07$Fjjm$\"1**RYeXtc'*F07$Fjjm$\"1+703YZ C6F[o7$Fjjm$\"1+)Go)pq)G\"F[o7$Fjjm$\"1+[GqET_9F[o7$Fjjm$\"1+cq(3/3i\" F[o7$Fjjm$\"1+s#)RL7pF[o7$Fjjm$\"1+?ZTOv.@F[o7$ Fjjm$\"1+K_XJKlAF[o7$Fjjm$\"1+o$QtV?T#F[o7$Fjjm$\"1+3Y/(3le#F[o7$Fjjm$ \"1+s=>MIMFF[o7$Fjjm$\"1+O^2QA1HF[o7$Fjjm$\"1*>r>[-%eIF[o7$Fjjm$\"1+O( ogm`A$F[o7$Fjjm$\"1+oNPfN%Q$F[o7$Fjjm$\"1++kCYC]NF[o7$Fjjm$\"1+'4WA#e- PF[o7$Fjjm$\"1+g&e!*)*o'QF[o7$Fjjm$\"1*z;vvxv.%F[o7$Fjjm$\"1*zG'*>ah=% F[o7$Fjjm$\"1*R%\\F3iYVF[o7$Fjjm$\"1*R]db)R7XF[o7$Fjjm$\"1*z_1szXn%F[o 7$Fjjm$\"1**f*oL)\\J[F[o7$Fjjm$\"1*RUjqHd+&F[o7$Fjjm$\"1*>V5M%Gi^F[o7$ Fjjm$\"1**R]H)R%H`F[o7$Fjjm$\"1***Ria24[&F[o7$Fjjm$\"1*>*\\o<]YcF[o7$F jjm$\"1*f&eEIJ-eF[o7$Fjjm$\"1*>fio!=lfF[o7$Fjjm$\"1**z3$f8W7'F[o7$Fjjm $\"1****f*o:6H'F[o7$Fjjm$\"1*>b:ip;X'F[o7$Fjjm$\"1*f(H1#feh'F[o7$Fjjm$ \"1+3Q=#*oynF[o7$Fjjm$\"1*R9')38$GpF[o7$Fjjm$\"1+s%*G*)z*4(F[o7$Fjjm$ \"1*RQ)3z<`sF[o7$Fjjm$\"1*f$*>*=r;uF[o7$Fjjm$\"1+#R.yMKd(F[o7$Fjjm$\"1 *********R9u(F[oFft-F$6$7S7$$\"\"'F(F(7$F`dn$\"1****fD!3)\\PF,7$F`dn$ \"1++cQf\\7qF,7$F`dn$\"1+g0tH&QJRF07$F`dn $\"1+S5j[V-VF07$F`dn$\"1++;#43]n%F07$F`dn$\"1**f\\M90M]F07$F`dn$\"1**R I'=(4g`F07$F`dn$\"1**R-@rzZdF07$F`dn$\"1**fT?)Hi2'F07$F`dn$\"1**z!y7v# ekF07$F`dn$\"1)*f$*f*\\kz'F07$F`dn$\"1+!3'[8[nrF07$F`dn$\"1**R!>(4z?vF 07$F`dn$\"1****>*\\K%*)yF07$F`dn$\"1)*zo)\\gzA)F07$F`dn$\"1)**z'z'3Jf) F07$F`dn$\"1(*Rq;]Rs*)F07$F`dn$\"1)*R1))[c-$*F07$F`dn$\"1**>V%Rd\"f'*F 07$F`dn$\"1+7hM_v-5F[o7$F`dn$\"1+%yM\\&zQ5F[o7$F`dn$\"1+!)3(HmO2\"F[o7 $F`dn$\"1+s=zVQ76F[o7$F`dn$\"1+'4!)Hur9\"F[o7$F`dn$\"1+?6i*>V=\"F[o7$F `dn$\"1++sa%zz@\"F[o7$F`dn$\"1+wxq\"yZD\"F[o7$F`dn$\"1+oz%*GS*G\"F[o7$ F`dn$\"1+w03dfD8F[o7$F`dn$\"1+SY)z!)4O\"F[o7$F`dn$\"1++!)3d-)R\"F[o7$F `dn$\"1+cMgVqL9F[o7$F`dn$\"1+GtM4>q9F[o7$F`dn$\"1+C'=Qvj]\"F[o7$F`dn$ \"1+K!38D'R:F[o7$F`dn$\"1+;K&4Lxd\"F[o7$F`dn$\"1+_28t\"=h\"F[o7$F`dn$ \"1+3A)>e\"[;F[o7$F`dn$\"1+wHi5%Ho\"F[o7$F`dn$\"1+++++K?Aus#F,7$Ff]o$\"1***z#\\#*)Q2$F,7$F f]o$\"1LL\\,KFIMF,7$Ff]o$\"1nme!QwTu$F,7$Ff]o$\"1++[7qb(4%F,7$Ff]o$\"1 ++?f')Q_WF,7$Ff]o$\"1++_*pMVz%F,7$Ff]o$\"1mm9jW&[5&F,7$Ff]o$\"1LL@'\\# 4uaF,7$Ff]o$\"1LLDVb)oy&F,7$Ff]o$\"1***fp@Q2:'F,7$Ff]o$\"1KLl&[4GZ'F,7 $Ff]o$\"1++'Hrsh#oF,7$Ff]o$\"1++[o!eE;(F,7$Ff]o$\"1MLLP_u8vF,7$Ff]o$\" 1LL*e&G:OyF,7$Ff]o$\"1mmEGI\"R=)F,7$Ff]o$\"1lm96\"Q^a)F,7$Ff]o$\"1***z wg&ef))F,7$Ff]o$\"1lm]=v>*>*F,7$Ff]o$\"1***RM&)\\+b*F,7$Ff]o$\"1)**zST !H$*)*F,7$Ff]o$\"1++cI$RD-\"F07$Ff]o$\"1+SY%p8%f5F07$Ff]o$\"1L`[]pa#4 \"F07$Ff]o$\"1++WtP#z7\"F07$Ff]o$\"1nm1L/)*f6F07$Ff]o$\"1+?JHo-&>\"F07 $Ff]o$\"1L$\\pv-!G7F07$Ff]o$\"1+?\"H5sCE\"F07$Ff]o$\"1LL,*=shH\"F07$Ff ]o$\"1+++cIXJ8F07$Ff]o$\"1n'QTsKaO\"F07$Ff]o$\"1nEgU==+9F07$Ff]o$\"1L8 #3AVYV\"F07$Ff]o$\"1+SQ['4jY\"F07$Ff]o$\"1L`#z%Hg-:F07$Ff]o$\"1n1p2T1N :F07$Ff]o$\"1+g\\uWnp:F07$Ff]o$\"1+?r65!Gg\"F07$Ff]o$\"1+++++SQ;F0Fft- F$6$7Y7$$!1+++++++5Fb\\l$\"1f*=sN\"**=DF[o7$$\"1'****\\([pbwF[o$\"16Q4 cN*\\Q$F[o7$$\"1+]7e?y,BFb\\l$\"1NE&z*Q$=J%F[o7$$\"1++DhPTHSFb\\l$\"16 /h'[*RmbF[o7$$\"1++vGt\\odFb\\l$\"1WY\"Q)3;qqF[o7$$\"1**\\iN`J*\\(Fb\\ l$\"1pg$Q\"[06))F[o7$$\"1)*\\7)=0S5*Fb\\l$\"1I4ho/^j5Fb\\l7$$\"1+DJeob w5!#:$\"1ULd>/cr7Fb\\l7$$\"1+DJ#)fR[7Fbio$\"1Wf!o2:M]\"Fb\\l7$$\"1+D1; So>9Fbio$\"1=L`'3&yXFb\\l7$$ \"1+]PiP1^Fbio$\"12XJ*>5,@Fbio$\"15w,g;chEFb\\l7$$\"1++vOgCqAFbio$\"1`!>'y>FJGFb \\l7$$\"1+DJxAwBCFbio$\"1//:M'H4&HFb\\l7$$\"1]7`1^.:DFbio$\"1zm914T/IF b\\l7$$\"1++vNzI1EFbio$\"1%3#y53tVIFb\\l7$$\"1+]P%*zi$o#Fbio$\"1+147g_ lIFb\\l7$$\"1*****H0[4w#Fbio$\"1)fY]\")4l2$Fb\\l7$$\"1DJqd)=f!GFbio$\" 1b!eIUuy2$Fb\\l7$$\"1\\iSi'*)3&GFbio$\"1XB`(4Nb2$Fb\\l7$$\"1u$4rYge*GF bio$\"147B!G+&pIFb\\l7$$\"1*\\7=FJ3%HFbio$\"1cE\")G+YRFbio$\"1hO]ts5#3#Fb\\l7$$\"1+DJ-heCTFbio$\"1;GZ _>*e#=Fb\\l7$$\"1*\\(=aT/!G%Fbio$\"1K(f$3#F[o7$$\"1*\\7=be5(fFbio$\"1LiV'RiNa\"F[ o7$$\"1+D1Z(p99'Fbio$\"1ETYo3P46F[o7$$\"1***\\(Qy23jFbio$\"1,@qjl,**yF 07$$\"1*\\7.2,D['Fbio$\"1Do?KBGQaF07$$\"1+++t8\\]mFbio$\"1^Yc/OfJPF07$ $\"1+](3s&GAoFbio$\"1E&)oyJ/&\\#F07$$\"1+D\"[_dE*pFbio$\"1$z7yDK]k\"F0 7$$\"1++]=:@\\rFbio$\"1^'3\"e5\"\\5\"F07$$\"1**\\Pi+kGtFbio$\"1:0E;u`x oF,7$$\"1+++\"4B\"*[(Fbio$\"1[:3'[0#GWF,7$$\"1+D\"=ZJ-m(Fbio$\"1g)*)yH ')Gs#F,7$$\"1+v=DS+CyFbio$\"1)fq8d()=o\"F,7$$\"\")F($\"1Mh1b;2X)*F[^o- Fgt6&FitFjtF(F(-%+AXESLABELSG6$%!GFehp-%*AXESSTYLEG6#%$BOXG-%%VIEWG6$; $!\"\"F_ipF\\hp%(DEFAULTG" 1 2 0 1 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 153 "The comparison sh ows that the discrete binomial distribution can be approximated reason ably by a continuous Gaussian distribution evaluated for discrete " } {TEXT 19 2 "nu" }{TEXT -1 31 " values for moderate values of " }{TEXT 306 1 "n" }{TEXT -1 9 " already." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 742 "We can also understand now how the two s eemingly different tasks can be put on a common footing: examples of m easurements where we obviously expect a broad distribution (height of \+ some species population with average and deviation the meaningful quan tities), and others where we are seeking a 'single number' only. In th e latter case the 'single number' is given by the average of the distr ibution, and the deviation becomes a measure for the uncertainty in th e measured result. While we would like in an ideal world the uncertain ty to be negligible compared to the 'single number', it is essential t o determine the uncertainty from the distribution of repeated measurem ents in order to assess the quality (accuracy) of the measured quantit y." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 324 "To really understand the concept of the probability distribution it is e ssential to analyze sequences of real data. We encourage the reader to perform 'real' experiments with coins, dies, etc.. It is possible to \+ analyze sequences produced with Maple's random number generator (whose quality falls short of state-of-the-art)." }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 16 "die:=rand(1..6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%$dieGR6\"6#%\"tGF&F&C%>%&_seedG-%%iremG6$,$F+\"-\"3p'>uU\"-*)****** ****>8$F+,&-F-6$F3\"\"'\"\"\"F8F8F&6#F+F&" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 6 "die();" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"%" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "We generate a sequence of tosses w ith three dies:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "M:=200; \+ n:=3;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MG\"$+#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "S1:=[seq([die(),die(),die()],i=1..M)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "S1[1],S1[2];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$ 7%\"\"$\"\"%\"\"'7%\"\"&F$F&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "W e can analyze the presence of particular faces:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 5 "F:=6;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"FG \"\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "DieAnalysis:=proc( S) local nu,C,i,toss,isucc,j; global F,M,n,p_s;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "for nu from 0 to n do: C[nu]:=0; od: p_s:=0;" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "for i from 1 to M do: toss:=S[i]: i succ:=0:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 22 " for j from 1 to n do: " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " if toss[j]=F then isucc:=isuc c+1; p_s:=p_s+1; fi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 26 " od: C[isuc c]:=C[isucc]+1;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "od: p_s:=evalf(p _s/(M*n)); print(\"Single-event probability: \",p_s);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "[seq(evalf(C[nu]/M),nu=0..n)]; end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "DieAnalysis(S1);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6$Q;Single-event~probability:~6\"$\"++++] " 0 "" {MPLTEXT 1 0 15 "ps:=evalf(1/6);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#psG$\"+nmmm;!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "[seq(BD(ps,n,nu),nu=0..n)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"+Oq.(y&!#5$\"+BAAsMF&$\"+ZWWWp!#6$\"+K'H'HY!# 7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "[seq(BD(p_s,n,nu),nu=0 ..n)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"+]i::c!#5$\"+]7GtNF&$\" ++vozv!#6$\"++]Pf`!#7" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 78 "We appea r to have results that are consistent with the binomial distribution. \+ " }}{PARA 0 "" 0 "" {TEXT -1 323 "We note that in addition to measurin g directly from the sample the nu-fold event probabilities we have the option to measure the single-event probability, and to generate nu-fo ld event probabilities from the binomial distribution. This allows one to 'beat' the statistics, i.e., to generate predictions for unlikely \+ events." }}{PARA 0 "" 0 "" {TEXT -1 211 "We leave it as an exercise to consider longer sequences of tosses to observe whether the consistenc y improves with the sample size. We can also look at the distribution \+ of frequencies of success for other faces:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 5 "F:=1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"FG\"\" \"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "DieAnalysis(S1);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6$Q;Single-event~probability:~6\"$\"+nmm m9!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"++++]i!#5$\"+++++JF&$\"+ ++++l!#6\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "[seq(BD(p_ s,n,nu),nu=0..n)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"+.Py8i!#5$ \"+BA)R?$F&$\"+ZW%o]&!#6$\"+lH'\\:$!#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "F:=2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"FG\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "DieAnalysis(S1);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6$Q;Single-event~probability:~6\"$\"++++ +:!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"++++]g!#5$\"++++]MF&$\"+ ++++X!#6$\"+++++]!#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "[se q(BD(p_s,n,nu),nu=0..n)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&$\"+++D Th!#5$\"+++D^KF&$\"+++]Pd!#6$\"++++vL!#7" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 337 "The scatter in the answers for the different faces demon strates that we have to be very careful when drawing conclusions from \+ measurements with small sample sizes (low number of repetitions)! Even the single-event probabilities vary substantially. The deviations wil l shrink only slowly when the sample size is increased (typically with " }{TEXT 19 9 "1/sqrt(M)" }{TEXT -1 3 " )." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 307 12 "Exercise B3:" }}{PARA 0 "" 0 " " {TEXT -1 256 "Determine the single-event probabilities for two diffe rent faces. Determine whether the measurement from samples approaches \+ the expected value of 1/6 when the sample size is increased by differe nt factors (e.g., 2, 10, 50). Is the approach consistent with " } {TEXT 19 9 "1/sqrt(M)" }{TEXT -1 11 " behaviour?" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 308 12 "Exercise B4:" }}{PARA 0 "" 0 "" {TEXT -1 82 "Analyze the case where four dies (or some other n umber) are tossed simultaneously." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 " " {TEXT 320 42 "Poisson distribution with chi-squared test" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 1401 "The Poisson distribution is discrete (a s the binomial distribution) and applies to the following case: suppos e a random process occurs, such as the decay of radioactive nuclei in \+ a sample. We cannot predict at what time any single nucleus will under go a transition, but we do know that nuclei in the sample decay with s ome rate, i.e., a certain number of events is expected per unit time. \+ If we pick some time interval, we expect on average a certain number o f decayed atoms (a fractional number as it is an average). In reality, when we record time sequences, we obtain different numbers of decays \+ in such intervals. This is due to the fact that the number of decays i s somehow distributed around the average. The Poisson distribution is \+ used in this situation as a simplified version of the binomial distrib ution. It depends on some average rate, and a chosen time interval, an d a decay is considered to be a success. The reason why we can simplif y from the binomial distribution is that the number of trials is extre mely large (a fraction of Avogadro's number), while the probability fo r success (decay of an individual nucleus) is extremely small (of the \+ order of the inverse of Avogadro's number per unit time). Thus, we can ignore the factor for the probability of no success raised to a huge \+ power (it is practically one), and use a mathematical property for the binomial coefficient for large " }{TEXT 310 1 "n" }{TEXT -1 1 "." }} {PARA 0 "" 0 "" {TEXT -1 37 "The Poisson distribution is given as:" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "PD:=(mu,nu)->exp(-mu)*mu^nu /nu!;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#PDGR6$%#muG%#nuG6\"6$%)ope ratorG%&arrowGF)*&*&-%$expG6#,$9$!\"\"\"\"\")F39%F5\"\"\"-%*factorialG 6#F7!\"\"F)F)F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 5 "Here " }{TEXT 19 2 "mu" }{TEXT -1 71 " is the expected mean numberof counts in the t ime interval chosen, and " }{TEXT 19 2 "nu" }{TEXT -1 58 " labels the \+ possibilities of observing 0, 1, 2,... counts." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 525 "Note that in an actual e xperiment with a Geiger-Mueller counter one does not observe all decay s, but a fraction thereof (depending on the efficiency of the tube for particular forms of radiation; the count rate, which leads to a certa in dead time for the tube, the supply voltage on the tube, etc.). This does not alter, however, the fact that one can determine the statisti cal nature of the radiation, as one simply looks at a subsample of the actual radiation sequence (one is also limited by the direction of ob servation)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 156 "Why is mu the mean (or average) count number for the chosen ti me interval? Let us calculate the probabilistic average number of obse rved decays (successes):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 " mu:='mu':" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "simplify(sum(n u*PD(mu,nu),nu=0..infinity));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%#muG " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "simplify(sum(PD(mu,nu), nu=0..infinity));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\"" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 84 "Note that the sum started at one, \+ as no contribution is made towards the average by " }{XPPEDIT 18 0 "nu = 0;" "6#/%#nuG\"\"!" }{TEXT -1 108 ". The mathematical steps require d in the calculation (done for us by Maple) are shown in J.R. Taylor's book." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 54 "The procedure when analyzing a sequence is as follows:" }}{PARA 0 "" 0 "" {TEXT -1 425 "First one can determine the average count rate, by \+ taking the total number of decays and divide by the total time interva l (which is subdivided into a certain number of sub-intervals). Then o ne arrives at an average number of decays per subinterval. Finally one forms a histogram by recording the frequency of intervals with 0, 1 , 2, ... counts. This histogram is to be compared with the predictions \+ of the PD for the given " }{TEXT 19 2 "mu" }{TEXT -1 112 " value. It i s possible to re-analyze the same decay sequence by using a different \+ choice of sub-interval length." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 65 "Let us look at the Poisson distribution f or two cases, one where " }{XPPEDIT 18 0 "mu < 1;" "6#2%#muG\"\"\"" } {TEXT -1 148 ", for which the number of zero-count intervals dominates , and another one, where one can see how the Gaussian distribution eme rges as a limit again." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "mu :=0.6;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "plot([seq([nu,PD(mu,nu)*t ,t=0..1],nu=0..10)],color=blue,axes=boxed);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#muG$\"\"'!\"\"" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "61-%'CURVESG6#7S7$\"\"!F(7$F($\"1FK0%=`i>\"!#<7$F ($\"1Ap;$o1rB#F,7$F($\"1m%*Q-`l2MF,7$F($\"1r5pxI'fe%F,7$F($\"1s*>=cq'e dF,7$F($\"1IXh;0#f%oF,7$F($\"1C\"H)o[qrzF,7$F($\"1:1McE*f8*F,7$F($\"1g ^\"ela'H5!#;7$F($\"1rrj*)=.\\6FE7$F($\"1#)f;A%zTD\"FE7$F($\"1x!o'R1bs8 FE7$F($\"13@w+zS\"\\\"FE7$F($\"1@6a&*)[fg\"FE7$F($\"1rQ-#4j*4FE7$F($\"1q0qA2Ig?FE7$F($\"1n+>/V=o@FE7$ F($\"1')=r\"p\\lG#FE7$F($\"10&3;`h#*R#FE7$F($\"1p@H7X'o^#FE7$F($\"1FrB g2'[i#FE7$F($\"1!4KMI\\8u#FE7$F($\"1=IIh\"[B'GFE7$F($\"1LwxA!yw'HFE7$F ($\"1kgxerV\"3$FE7$F($\"17>gx9'*)>$FE7$F($\"1`ZJ&4ORJ$FE7$F($\"1Bm$)\\ *z^U$FE7$F($\"1z^s$G(p[NFE7$F($\"1.%z=5$ofOFE7$F($\"1%f2A:%=yPFE7$F($ \"1rS_;Pc&)QFE7$F($\"1a;)H%z&H+%FE7$F($\"1SALhnT8TFE7$F($\"1+/x`!y)GUF E7$F($\"1J:`rDwTVFE7$F($\"1?[.+A%*fWFE7$F($\"18bEGKwtXFE7$F($\"1in;p;; !p%FE7$F($\"1V!4AF'f0[FE7$F($\"1.3$Qoo;\"\\FE7$F($\"1`[\"yjRK.&FE7$F($ \"1\\gJ[T(>9&FE7$F($\"1Lj]:y!zD&FE7$F($\"1Y+\"y7r)o`FE7$F($\"1,++hj6)[ &FE-F$6#7S7$$\"\"\"F(F(7$Fit$\"1X7>0\">v<(!#=7$Fit$\"1/K15SEU8F,7$Fit$ \"1Y?o\"=$fW?F,7$Fit$\"1*))[p%yd^FF,7$Fit$\"1,<^PB?bMF,7$Fit$\"1\"on/J _v5%F,7$Fit$\"1*[y=#H-$y%F,7$Fit$\"1V-Z%f&f\"[&F,7$Fit$\"1=9kNz#z<'F,7 $Fit$\"1'\\g'Q8>%*oF,7$Fit$\"1(**4R`w]_(F,7$Fit$\"1S)3!RQIN#)F,7$Fit$ \"1a'fcSZ%[*)F,7$Fit$\"1>sTuLpN'*F,7$Fit$\"1_*Q`&y(f-\"FE7$Fit$\"1Y)*R >w=+6FE7$Fit$\"1*>qO-`I;\"FE7$Fit$\"120xL/=O7FE7$Fit$\"1n?n#e54I\"FE7$ Fit$\"1'y$>:)H>P\"FE7$Fit$\"1t*R\">pbR9FE7$Fit$\"1U(evq=,^\"FE7$Fit$\" 1)eLjX;\\d\"FE7$Fit$\"1c!fAe4[k\"FE7$Fit$\"1K/*p*)3ur\"FE7$Fit$\"1yG)Q \"og!y\"FE7$Fit$\"1G-\\&Hi)[=FE7$Fit$\"1-Vz')oP>>FE7$Fit$\"1(QIulh$))> FE7$Fit$\"1<;:qz5b?FE7$Fit$\"1`P\\q$=#H@FE7$Fit$\"1xVRh)4e>#FE7$Fit$\" 1G**f\"\\5pE#FE7$Fit$\"1Tw>I#Q8L#FE7$Fit$\"1Z23mZx,CFE7$Fit$\"1\\\"*4d +0oCFE7$Fit$\"1g/dKoKPDFE7$Fit$\"1m`BVv00EFE7$Fit$\"1`fM?`'fn#FE7$Fit$ \"1lEHPzDWFFE7$Fit$\"1)*=%=+(49GFE7$Fit$\"1!ovOwdL)GFE7$Fit$\"1nkl57+Z HFE7$Fit$\"1dd0$yV*>IFE7$Fit$\"1+WO*[%=&3$FE7$Fit$\"1@qo*oWZ:$FE7$Fit$ \"1Ot2xEK@KFE7$Fit$\"1+++<)pGH$FE-F$6#7S7$$\"\"#F(F(7$F`^l$\"1-%R8tbK: #F^u7$F`^l$\"1&)>yH?zESF^u7$F`^l$\"1A_UW&zP8'F^u7$F`^l$\"1]5,SNta#)F^u 7$F`^l$\"1\"e[6qgl.\"F,7$F`^l$\"1kb,$plAB\"F,7$F`^l$\"1$H=k(o!\\V\"F,7 $F`^l$\"1/Yh(F,7$F`^l$\"1y\\\"*GE<:yF,7$F`^l$\"11_Agf*y-)F,7$F`^l$\"1-Y/6QxK#)F ,7$F`^l$\"1#4rY+\"HU%)F,7$F`^l$\"1/9:!Ht+l)F,7$F`^l$\"1RW2JO+T))F,7$F` ^l$\"1d,D[8$)f!*F,7$F`^l$\"1ti:nMbb#*F,7$F`^l$\"15I5oSBk%*F,7$F`^l$\"1 RPDI!oRm*F,7$F`^l$\"1+++]%4'y)*F,-F$6#7S7$$\"\"$F(F(7$Ffgl$\"10)yEY6lI %!#>7$Ffgl$\"1qRcfSe`!)F[hl7$Ffgl$\"1W]))3fvE7F^u7$Ffgl$\"15A+3n%4l\"F ^u7$Ffgl$\"1irH-97t?F^u7$Ffgl$\"1F6.'QJXY#F^u7$Ffgl$\"1'eOGv8)pGF^u7$F fgl$\"14#\\jNd*)G$F^u7$Ffgl$\"1@'45wcnq$F^u7$Ffgl$\"1kvx-[^OTF^u7$Ffgl $\"1Y*)))>f/:XF^u7$Ffgl$\"19^!HI#=T\\F^u7$Ffgl$\"1*G_GWo!p`F^u7$Ffgl$ \"1'3lS-;9y&F^u7$Ffgl$\"1k0TJr'e:'F^u7$Ffgl$\"1^3t:d7,mF^u7$Ffgl$\"1!z 989=$ypF^u7$Ffgl$\"1=A(=g#38E;:\"F,7$Ffgl$\"1kuL%*p,$>\"F ,7$Ffgl$\"1\\h'>ykIB\"F,7$Ffgl$\"1Hp;A5`x7F,7$Ffgl$\"1fKq;f[<8F,7$Ffgl $\"1rA#[HY,O\"F,7$Ffgl$\"1')p(z$H!))R\"F,7$Ffgl$\"1rD]f[1T9F,7$Ffgl$\" 1(e4T.I3[\"F,7$Ffgl$\"1m\")Q*4'RA:F,7$Ffgl$\"1'*HyDX.j:F,7$Ffgl$\"1T]/ #>zbg\"F,7$Ffgl$\"1@*3Awalk\"F,7$Ffgl$\"1=U$4?e%)o\"F,7$Ffgl$\"1\"GI!e Y,I\"=F,7$Ffgl$\"1b7V$p56&=F, 7$Ffgl$\"1-1i8o%G*=F,7$Ffgl$\"1[201OzK>F,7$Ffgl$\"1+++!*=sv>F,-F$6#7S7 $$\"\"%F(F(7$F]am$\"11#=S>n(fk!#?7$F]am$\"1&fM*3w.37F[hl7$F]am$\"1mvKj Q8S=F[hl7$F]am$\"1:L+i+UwCF[hl7$F]am$\"1UdW.@o4JF[hl7$F]am$\"1\"pY!zqz 'p$F[hl7$F]am$\"1y[DH1s/VF[hl7$F]am$\"17Q_MgVL\\F[hl7$F]am$\"1JW^T^8gb F[hl7$F]am$\"1Xj;/Ax/iF[hl7$F]am$\"1=M$)z)oDx'F[hl7$F]am$\"1rwNaMx6uF[ hl7$F]am$\"1L%yUm-O0)F[hl7$F]am$\"1Gw4OS7s')F[hl7$F]am$\"1Xe6(p+QB*F[h l7$F]am$\"1tift&)o,**F[hl7$F]am$\"1=s>@xuY5F^u7$F]am$\"1K3G!RiD6\"F^u7 $F]am$\"1SjGC&>3<\"F^u7$F]am$\"1<%\\N$otM7F^u7$F]am$\"1B[4F7g&H\"F^u7$ F]am$\"1(Glm$o5f8F^u7$F]am$\"1Xnb5[U<9F^u7$F]am$\"1)H$)QiG.[\"F^u7$F]a m$\"1B\\82!oca\"F^u7$F]am$\"1mBLKha-;F^u7$F]am$\"1jFxlg(Rm\"F^u7$F]am$ \"10+%z>Rus\"F^u7$F]am$\"1(>1:\\D&*y\"F^u7$F]am$\"1B#\\H<(f\\=F^u7$F]a m$\"1$R]K`'H;>F^u7$F]am$\"1))[0v)Gi(>F^u7$F]am$\"11MBU%>-/#F^u7$F]am$ \"1ya'pS/#)4#F^u7$F]am$\"1cQD*G(fh@F^u7$F]am$\"1!Qk60X7A#F^u7$F]am$\"1 \\A3\\Tf$G#F^u7$F]am$\"1$\\u')y^XM#F^u7$F]am$\"1iv1)yo$3CF^u7$F]am$\"1 \"Q8L9K)pCF^u7$F]am$\"1F8S,toKDF^u7$F]am$\"1@a/()>-&f#F^u7$F]am$\"1JBK *3,Bl#F^u7$F]am$\"1Z]Z/%\\zr#F^u7$F]am$\"1#)o9SgmwFF^u7$F]am$\"1.4V?-F RGF^u7$F]am$\"1@h24/>**GF^u7$F]am$\"1+++NGejHF^u-F$6#7S7$$\"\"&F(F(7$F djm$\"1Z=#Gj?T8(>BVm*Fbam7$Fdjm$\"1:c3;![&=F[hl7$Fdjm$\"1S))ze`0B>F[hl7$Fdjm$\"1;t#*y7 x'*>F[hl7$Fdjm$\"11!Gv.FH2#F[hl7$Fdjm$\"1Ou!)*eIu9#F[hl7$Fdjm$\"1n!Rvg ;&>AF[hl7$Fdjm$\"1s/!*Reb*H#F[hl7$Fdjm$\"1le1]YZrBF[hl7$Fdjm$\"1(3!oIL E[CF[hl7$Fdjm$\"1u&e$)GXy^#F[hl7$Fdjm$\"1FY5Zn\"Rf#F[hl7$Fdjm$\"1csRhS \\lEF[hl7$Fdjm$\"1)p)*)yHJSFF[hl7$Fdjm$\"1#R4k9iM\"GF[hl7$Fdjm$\"1u5oX D/!*GF[hl7$Fdjm$\"1cg(>d)zjHF[hl7$Fdjm$\"1#f\"ohZARIF[hl7$Fdjm$\"10XX% QES6$F[hl7$Fdjm$\"1(z'=28w#=$F[hl7$Fdjm$\"1c+P&GR:E$F[hl7$Fdjm$\"1ei<[ #**>L$F[hl7$Fdjm$\"1$3BVm*Fijm7$F[dn$\"1#4RS)[lS5Fb am7$F[dn$\"1EBM$3c!36Fbam7$F[dn$\"1Tu\\GE?)=\"Fbam7$F[dn$\"1fM3ls4c7Fb am7$F[dn$\"1)ehz'[2N8Fbam7$F[dn$\"1R&[(GM)\\S\"Fbam7$F[dn$\"1jE%)*>%o \"[\"Fbam7$F[dn$\"19m2s9sa:Fbam7$F[dn$\"1U(RN?G4j\"Fbam7$F[dn$\"18)*=s (44q\"Fbam7$F[dn$\"1_/;[VRwFbam7$F[dn$\"1TeOy7x'*>Fbam7$F[dn$\"1=^%p.FH2#Fbam7$F[dn$\"1(f.#*e Iu9#Fbam7$F[dn$\"1f\\\"pg;&>AFbam7$F[dn$\"1dQDReb*H#Fbam7$F[dn$\"1G!*R \\YZrBFbam7$F[dn$\"1d;**HLE[CFbam7$F[dn$\"1y0l(GXy^#Fbam7$F[dn$\"1T_PY n\"Rf#Fbam7$F[dn$\"1VxkgS\\lEFbam7$F[dn$\"1Z\"G\"yHJSFFbam7$F[dn$\"1r# =c9iM\"GFbam7$F[dn$\"1?%o[aU+*GFbam7$F[dn$\"1jE9r&)zjHFbam7$F[dn$\"1*) p#3wC#RIFbam7$F[dn$\"1o)yNQES6$Fbam7$F[dn$\"1L=H18w#=$Fbam7$F[dn$\"1SH X%GR:E$Fbam7$F[dn$\"1H$RsC**>L$Fbam7$F[dn$\"1I!fNECrS$Fbam7$F[dn$\"1tI \"**[G!zMFbam7$F[dn$\"1+++,%*HcNFbam-F$6#7S7$$\"\"(F(F(7$Fb]o$\"1&4ZD( >LWm!#B7$Fb]o$\"1K5[pHbU7F`dn7$Fb]o$\"13!4q64F*=F`dn7$Fb]o$\"1lX^@\\$F`dn7$Fb]o$\"1XZi`)>C!QF`dn7$Fb]o$\"1=6F[ErF WF`dn7$Fb]o$\"1%)\\)4N\"Ru]F`dn7$Fb]o$\"14y*o9'**=dF`dn7$Fb]o$\"1]^\\6 90#Q'F`dn7$Fb]o$\"1\"*4Y\\32mpF`dn7$Fb]o$\"1bG)>TQNi(F`dn7$Fb]o$\"1>Y[ *f0PG)F`dn7$Fb]o$\"1Dxk5!**)>*)F`dn7$Fb]o$\"1Q=SwNi(\\*F`dn7$Fb]o$\"1$ )*zjRf%=5Fijm7$Fb]o$\"1**z>*zam2\"Fijm7$Fb]o$\"1]TT())\\V9\"Fijm7$Fb]o $\"18n(4^rU?\"Fijm7$Fb]o$\"1$G\"4H],q7Fijm7$Fb]o$\"1&3C`S=EL\"Fijm7$Fb ]o$\"1$e(*3YQzR\"Fijm7$Fb]o$\"1n&GEmAzX\"Fijm7$Fb]o$\"17'zxsBE_\"Fijm7 $Fb]o$\"1f%o0&*H)*e\"Fijm7$Fb]o$\"1n42lJL[;Fijm7$Fb]o$\"14;[R#=:r\"Fij m7$Fb]o$\"1$y6Sg%zwFij m7$Fb]o$\"1,PBjy/r>Fijm7$Fb]o$\"1EcndDpK?Fijm7$Fb]o$\"1a*)yE9^)4#Fijm7 $Fb]o$\"1gI\"[5`\"e@Fijm7$Fb]o$\"1d\\K7sNBAFijm7$Fb]o$\"1H(=g>4ZG#Fijm 7$Fb]o$\"1z#\\npR)[BFijm7$Fb]o$\"13$*F^go7$Fi fo$\"1F:gP=`>9Fg]o7$Fifo$\"1>q#4>\"Q5>Fg]o7$Fifo$\"1Wz&Hi(*))R#Fg]o7$F ifo$\"11n!**)[\"=&GFg]o7$Fifo$\"1.-%e[%y?LFg]o7$Fifo$\"1mD#G^$z0QFg]o7 $Fifo$\"1>Vq4rC*G%Fg]o7$Fifo$\"1Wz4e&Qly%Fg]o7$Fifo$\"1FW_OJbC_Fg]o7$F ifo$\"1.>O3QlzF@'Fg]o7$Fifo$\"1PUDdU#**o'Fg]o7$F ifo$\"1SC_\"o5$*4\"\\2)Fg ]o7$Fifo$\"1)fnY:CEe)Fg]o7$Fifo$\"1ywLJj.K!*Fg]o7$Fifo$\"1KI9W$4\" F`dn7$Fifo$\"1K)4dzn>9\"F`dn7$Fifo$\"1cfziCP#>\"F`dn7$Fifo$\"1R!oO()\\ iB\"F`dn7$Fifo$\"1Q3Zz'QOG\"F`dn7$Fifo$\"1:J'G&ffK8F`dn7$Fifo$\"1j,i$4 \"\\!Q\"F`dn7$Fifo$\"1g[-0@$oU\"F`dn7$Fifo$\"1AOE(*eGy9F`dn7$Fifo$\"15 +4=%>X_\"F`dn7$Fifo$\"13'>*pN)Qd\"F`dn7$Fifo$\"1'zK%G[h=;F`dn7$Fifo$\" 1r814z^n;F`dn7$Fifo$\"1pm#o*=`8F`dn7$Fifo$\"1 #[RqM(y`>F`dn7$Fifo$\"1kBa/T(=+#F`dn7$Fifo$\"1u+fS31Y?F`dn7$Fifo$\"1b8 dEQq'4#F`dn7$Fifo$\"1JC&p^**>9#F`dn7$Fifo$\"1E/KqTH!>#F`dn7$Fifo$\"19D &eJ=lB#F`dn7$Fifo$\"1+++tC>'G#F`dn-F$6#7S7$$\"\"*F(F(7$F``p$\"1*o$=&)f ;AL!#D7$F``p$\"1BqOX[w7iFe`p7$F``p$\"1j/%>eXNY*Fe`p7$F``p$\"1u%R.Y(et7 F^go7$F``p$\"1'>A;3l#*f\"F^go7$F``p$\"1o')=E*47!>F^go7$F``p$\"1*G4MKcQ @#F^go7$F``p$\"1[,muc>PDF^go7$F``p$\"1G3^s!)\\fGF^go7$F``p$\"1R2q/d-\" >$F^go7$F``p$\"1kyeBa.$[$F^go7$F``p$\"1_4u/#p<\"QF^go7$F``p$\"1[NQ)z_= 9%F^go7$F``p$\"1[2'Q]\\*fWF^go7$F``p$\"1TIk'y6)[ZF^go7$F``p$\"1_$H-)pH #4&F^go7$F``p$\"1wRA%*RF$Q&F^go7$F``p$\"1&p$>N%\\-Bp *)pF^go7$F``p$\"1J9v5Lh*G(F^go7$F``p$\"1E0SO'=Jh(F^go7$F``p$\"1GXB](\\ \"\\zF^go7$F``p$\"1.6lAemT#)F^go7$F``p$\"1w1g%>\"fd&)F^go7$F``p$\"1tW9 W<\"Fg]o7$F``p$\"1%H07\\pd?\"Fg]o7$F``p$\"1`'=Bm*eQ7Fg]o7$F``p$ \"1#oqyQ*>q7Fg]o7$F``p$\"1a9u(*[_-8Fg]o7$F``p$\"1(ppgt#eM8Fg]o7$F``p$ \"1h]4g0/k8Fg]o7$F``p$\"1l=u]D!yR\"Fg]o7$F``p$\"1yflx'**zU\"Fg]o7$F``p $\"1+wA8h>g9Fg]o7$F``p$\"1_AC5A,\"\\\"Fg]o7$F``p$\"1+++:$GT_\"Fg]o-F$6 #7S7$$\"#5F(F(7$Fgip$\"1&=G8f*H$*>!#E7$Fgip$\"1Tyi2*ews$F\\jp7$Fgip$\" 1$>&y\\t7ycF\\jp7$Fgip$\"1gC(GwC:k(F\\jp7$Fgip$\"1vCy!\\!f&f*F\\jp7$Fg ip$\"1Uz$e&fsS6Fe`p7$Fgip$\"1F3>%z8$G8Fe`p7$Fgip$\"1eD'\\SFe`p7$Fgip$\"1W7Qa7#)*3#Fe`p7$ Fgip$\"1mY4B:1(G#Fe`p7$Fgip$\"1\")=Iz;6&[#Fe`p7$Fgip$\"1r!4Eqpfn#Fe`p7 $Fgip$\"1+uHsqG\\GFe`p7$Fgip$\"1C%Fe`p7$Fgip$\"1R\"Hp)zwtVF e`p7$Fgip$\"1A)RB=ryc%Fe`p7$Fgip$\"1qAm])*[pZFe`p7$Fgip$\"1398%\\**\\% \\Fe`p7$Fgip$\"1z=K_ &Fe`p7$Fgip$\"1K`Z>%Gtq&Fe`p7$Fgip$\"1tyS)eVJ\"fFe`p7$Fgip$\"1-Kprw2)4 'Fe`p7$Fgip$\"1,+**yU`&H'Fe`p7$Fgip$\"1*=BIJfWZ'Fe`p7$Fgip$\"1(4;bjr+n 'Fe`p7$Fgip$\"1hrb'eFT&oFe`p7$Fgip$\"1Knq)3>l/(Fe`p7$Fgip$\"1#)G-[phMs Fe`p7$Fgip$\"1qXsuz`JuFe`p7$Fgip$\"1([d!Gj>@wFe`p7$Fgip$\"1CLI(Q\\^\"y Fe`p7$Fgip$\"1>QH " 0 "" {MPLTEXT 1 0 8 "mu:=3.4;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "plot([seq([nu,PD(mu,nu)*t,t=0..1],nu=0..10)],color=blue,axes=boxed); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#muG$\"#M!\"\"" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "61-%'CURVESG6#7S7$\"\"!F(7$F($\" 1b!RO5BWF(!#>7$F($\"1V')*\\(fQg8!#=7$F($\"1$z\"RBt>s?F07$F($\"16dX>qs) y#F07$F($\"1%Q]\\<^=]$F07$F($\"1`*RLb3I;%F07$F($\"16C$Hs)fZ[F07$F($\"1 *GTc&GgbbF07$F($\"1!Q<,TO8E'F07$F($\"1old;,F()pF07$F($\"1u(f.#HnEwF07$ F($\"1suFN!*[Y$)F07$F($\"1O:K%yg#p!*F07$F($\"12Q(3H&yl(*F07$F($\"1^j5F 'H)R5!#<7$F($\"1(3G/IT]6\"FX7$F($\"1HIE`avy6FX7$F($\"1dhrf,(GD\"FX7$F( $\"1$*3=LUZ=8FX7$F($\"11!e%*H_/R\"FX7$F($\"1gP))RC**e9FX7$F($\"1U^3')o ]I:FX7$F($\"1+?4s%zhf\"FX7$F($\"1(fe\"Gi,n;FX7$F($\"1)=/%4dfSFX7$F($\"1z*y[h1_,# FX7$F($\"1Q!zC2aG3#FX7$F($\"1<\"o0/lz:#FX7$F($\"1<94'ebaA#FX7$F($\"1Ug +\\h^(H#FX7$F($\"1IAJHP\"GO#FX7$F($\"1-!Q>B,UV#FX7$F($\"1\"o]&38P,DFX7 $F($\"1i()y#R$erDFX7$F($\"1j?tQ&G-k#FX7$F($\"1zO$[k$47FFX7$F($\"1,+Q$G 38y#FX7$F($\"1l`=[-4_GFX7$F($\"1)4&e-hGAHFX7$F($\"1P!>D!))y')HFX7$F($ \"1&Gpd<;21$FX7$F($\"1@*)p*pPo7$FX7$F($\"1*e>`.Pt>$FX7$F($\"1+6-&*Q\"[ E$FX7$F($\"1+++'*pKPLFX-F$6#7S7$$\"\"\"F(F(7$Fjt$\"1w>-cQItCF07$Fjt$\" 1VGY;BJDYF07$Fjt$\"1#Rn \"FX7$Fjt$\"1^Ye3HU:9FX7$Fjt$\"1M*>jc$=[;FX7$Fjt$\"1Et^r\\!*))=FX7$Fjt $\"1C`6!Qa)G@FX7$Fjt$\"1b(*QS=nvBFX7$Fjt$\"1=]u$zoIf#FX7$Fjt$\"1$o%*GF 1y$GFX7$Fjt$\"1Iwkn'[N3$FX7$Fjt$\"1P0&***pO?LFX7$Fjt$\"1mKGL2UNNFX7$Fj t$\"1,$eEUS6z$FX7$Fjt$\"1LeOU&ox+%FX7$Fjt$\"18kQW&e(fUFX7$Fjt$\"1#GPUR 7G[%FX7$Fjt$\"16rl>y`FZFX7$Fjt$\"15'yrHu0'\\FX7$Fjt$\"1^/M9Ms._FX7$Fjt $\"14Y$o?5qU&FX7$Fjt$\"1au`x^&ym&FX7$Fjt$\"1K=&QTD!=fFX7$Fjt$\"1RS)*R$ *zNhFX7$Fjt$\"1$[=j9,5P'FX7$Fjt$\"1v_xEw)Rh'FX7$Fjt$\"1\\Bw#\\-<&oFX7$ Fjt$\"1>bn[Qq\"3(FX7$Fjt$\"15%f-9\"3PtFX7$Fjt$\"1M96&**[lc(FX7$Fjt$\"1 #*))44\\b6yFX7$Fjt$\"1m.@#omN.)FX7$Fjt$\"1(*\\@\">%Gw#)FX7$Fjt$\"1k0x^ Cm/&)FX7$Fjt$\"1wdDQNQV()FX7$Fjt$\"1[q`Mqxw*)FX7$Fjt$\"1l+O&R=6A*FX7$F jt$\"1@A\\m\"[kX*FX7$Fjt$\"18o!p%o5(p*FX7$Fjt$\"1F;%>vsd$**FX7$Fjt$\"1 g'yJ>3b,\"!#;7$Fjt$\"1><4+NkS5F]]l7$Fjt$\"1G\\B=[7j5F]]l7$Fjt$\"1fNN#f %4(3\"F]]l7$Fjt$\"1_$fYsO+6\"F]]l7$Fjt$\"1+++z6pM6F]]l-F$6#7S7$$\"\"#F (F(7$Fa^l$\"1XMealh/UF07$Fa^l$\"1^RwO4.jyF07$Fa^l$\"1#=9J0Ix>\"FX7$Fa^ l$\"19h]U%*=iS#FX7$Fa^l$\" 1E\"3B17>!GFX7$Fa^l$\"1Z+3^%Q6@$FX7$Fa^l$\"11KjXC0>OFX7$Fa^l$\"1#[M!G@ kQSFX7$Fa^l$\"1^4o[p@3WFX7$Fa^l$\"1w1?V115F]]l7$Fa^l$\"1j]j()e3V5F]]l7$ Fa^l$\"1+dq%>qI3\"F]]l7$Fa^l$\"1IpP'*yPC6F]]l7$Fa^l$\"1[%)eB%*yk6F]]l7 $Fa^l$\"11w3a'*)Q?\"F]]l7$Fa^l$\"1:,l$z.tC\"F]]l7$Fa^l$\"1#*)o*GLJ'G\" F]]l7$Fa^l$\"1$GSVLkzK\"F]]l7$Fa^l$\"1kLwNjql8F]]l7$Fa^l$\"1MxG7$opS\" F]]l7$Fa^l$\"1Uhd;EzX9F]]l7$Fa^l$\"1k+k@=\"F]]l7$Fggl$\"1t *[WbzuA\"F]]l7$Fggl$\"1)\\$p\\\\Hu7F]]l7$Fggl$\"1y*\\0o%4?8F]]l7$Fggl$ \"11>8v#3WO\"F]]l7$Fggl$\"1*3am'4h89F]]l7$Fggl$\"1#)[Ot<#yX\"F]]l7$Fgg l$\"1)ROhCE]]\"F]]l7$Fggl$\"1U:V90!ya\"F]]l7$Fggl$\"1SH%y3kXf\"F]]l7$F ggl$\"1rD0m\\cQ;F]]l7$Fggl$\"1+%e?ceXo\"F]]l7$Fggl$\"1a!=!4d_H#=F]]l7$Fggl$\"1jr*GD4$o=F]]l7$Fgg l$\"1zoFF]]l7$Fggl$\"1,f$fyXl&>F]]l7$Fggl$\"1=>$R9t\\+#F]]l7$Fggl $\"1\"GI-:(G[?F]]l7$Fggl$\"1&*=/:*oW4#F]]l7$Fggl$\"1Bqic2nQ@F]]l7$Fggl $\"1+++s;<'=#F]]l-F$6#7S7$$\"\"%F(F(7$F]am$\"1L>yzsW]SF07$F]am$\"1yiH+ )>Zd(F07$F]am$\"1>r<%G8Q:\"FX7$F]am$\"1NIw&z\"y_:FX7$F]am$\"1;mlzS&)\\ >FX7$F]am$\"15jh:4*zJ#FX7$F]am$\"1T[$fGv\"*p#FX7$F]am$\"1U\"3qq'R$4$FX 7$F]am$\"1*Q%)\\&QN'[$FX7$F]am$\"1v6D#ee0*QFX7$F]am$\"1!)4/EBeYUFX7$F] am$\"1CRbP2QZYFX7$F]am$\"1gC&e'[#)\\]FX7$F]am$\"1w;baVlPaFX7$F]am$\"1# e>]oS)*y&FX7$F]am$\"16r8'HC'3iFX7$F]am$\"1j.pY*)QjlFX7$F]am$\"1tI:tX1w pFX7$F]am$\"1#[mxU_8M(FX7$F]am$\"1uyG*4J@u(FX7$F]am$\"1+hZqrwB\")FX7$F ]am$\"1=\"3XPk>_)FX7$F]am$\"1?N_6Pj()))FX7$F]am$\"1khq\"\\d?G*FX7$F]am $\"1\")>`[Hv\"p*FX7$F]am$\"1cs![SR[+\"F]]l7$F]am$\"1dyg?wNV5F]]l7$F]am $\"12+a12:$3\"F]]l7$F]am$\"1vGux/3A6F]]l7$F]am$\"1!>8J.Z(f6F]]l7$F]am$ \"1Q+)3Kp:?\"F]]l7$F]am$\"1W*fl][\"R7F]]l7$F]am$\"1A)*Q3BFz7F]]l7$F]am $\"1;sTO/j:8F]]l7$F]am$\"1O7zt%z`N\"F]]l7$F]am$\"1q_C?-y#R\"F]]l7$F]am $\"1y\\#oxu=V\"F]]l7$F]am$\"1,gr^o4q9F]]l7$F]am$\"1i$*)=#=65:F]]l7$F]a m$\"1F$eS*4l[:F]]l7$F]am$\"1lq$R'G1)e\"F]]l7$F]am$\"1%eMOQ[ri\"F]]l7$F ]am$\"1f:(p\"R1j;F]]l7$F]am$\"1\"fT79qy\"=F]]l7$F]am$\"1+++?fCe=F]]l- F$6#7S7$$\"\"&F(F(7$Fcjm$\"1,O7^TIaFF07$Fcjm$\"1y>$eY43:&F07$Fcjm$\"1t !))[LIf%yF07$Fcjm$\"16JX@;*e0\"FX7$Fcjm$\"1=ieu2!fK\"FX7$Fcjm$\"1a!=J# QBw:FX7$Fcjm$\"12n,&>Ra$=FX7$Fcjm$\"169Vh(4N5#FX7$Fcjm$\"1Y)R@i?2P#FX7 $Fcjm$\"1o#on$)zbk#FX7$Fcjm$\"1s>i#)fn()GFX7$Fcjm$\"1[q`-*=-;$FX7$Fcjm $\"1#om)44)QV$FX7$Fcjm$\"1KcEi\\g(p$FX7$Fcjm$\"1='fqm\"4PRFX7$Fcjm$\"1 $o4F_k=A%FX7$Fcjm$\"1gA;&[/JY%FX7$Fcjm$\"1NdC6RsVZFX7$Fcjm$\"1$[hCl>@* \\FX7$Fcjm$\"1+B=\\\"\\YE&FX7$Fcjm$\"1WCnx;;CbFX7$Fcjm$\"1I**\\wd$\\z& FX7$Fcjm$\"1<\"\\d#4fVgFX7$Fcjm$\"1<#ej4*z6jFX7$Fcjm$\"1jz32?R!f'FX7$F cjm$\"1CB&[&z!H$oFX7$Fcjm$\"1H$zB#=$[4(FX7$Fcjm$\"1/O!o![UltFX7$Fcjm$ \"11p1^s9IwFX7$Fcjm$\"1HhmF)zi)yFX7$Fcjm$\"1>2d%Qr1<)FX7$Fcjm$\"1r\\FZ )4iU)FX7$Fcjm$\"14X!)*p^!*p)FX7$Fcjm$\"1]!p/(pGY*)FX7$Fcjm$\"1[zpC/e;# *FX7$Fcjm$\"1))em+b!4Z*FX7$Fcjm$\"1e?\\&[[nt*FX7$Fcjm$\"1(GL[fem***FX7 $Fcjm$\"1G)4#Rg(o-\"F]]l7$Fcjm$\"1?IHwE3`5F]]l7$Fcjm$\"1Y!>yu#))z5F]]l 7$Fcjm$\"1r$*)F07$Ficn$\"1'z:.@#3S5FX7$Ficn$\"1v\"*ek)))> >\"FX7$Ficn$\"1fdH_$3MM\"FX7$Ficn$\"1iGA2>;*\\\"FX7$Ficn$\"1#3\"=B(\\j j\"FX7$Ficn$\"1,U562z!z\"FX7$Ficn$\"1<@c[e'e%>FX7$Ficn$\"1FqAy%4`4#FX7 $Ficn$\"1]$=4h=5B#FX7$Ficn$\"1rmv&*)*Q#R#FX7$Ficn$\"1557TD4HDFX7$Ficn$ \"1N?(e@5\")o#FX7$Ficn$\"1\\Z`-y')GGFX7$Ficn$\"1@6)R=,L)HFX7$Ficn$\"1' [JMGe.8$FX7$Ficn$\"1d^SRpz$G$FX7$Ficn$\"1fM&R_,ZU$FX7$Ficn$\"1\")p$R:' owNFX7$Ficn$\"1ux)*pabMPFX7$Ficn$\"1A`pt6)>(QFX7$Ficn$\"1`'*f)p//-%FX7 $Ficn$\"1!>yIsSP<%FX7$Ficn$\"1'zK\"3,vBVFX7$Ficn$\"1&*>^,>*)oWFX7$Ficn $\"1a71hFX7$Ficn$\"1dev0^%*piFX7$Ficn$\"13dcF'R$3kFX7$Ficn$\"1 uNPn`&pc'FX7$Ficn$\"1!pHf*)\\f@F07$F_]o$\"1KN)[)f@1HF07$F_]o$\"1 a(zW.%Q\\OF07$F_]o$\"1.we[fRQVF07$F_]o$\"1j'=l(y#=0&F07$F_]o$\"1'))fp> g'*y&F07$F_]o$\"1s>zzi7DlF07$F_]o$\"1*ediaV;G(F07$F_]o$\"1N9MmV)z%zF07 $F_]o$\"1[=O2j7)p)F07$F_]o$\"1ZUE.%[8X*F07$F_]o$\"1fE8Yx#pln!)=FX7$F_]o$\"1)\\Gc*zw_>FX7$F_] o$\"1V0?$\\Xs-#FX7$F_]o$\"1b;[%>2,5#FX7$F_]o$\"13+&pj/1<#FX7$F_]o$\"15 iKF+))[AFX7$F_]o$\"1dy,ZR@>BFX7$F_]o$\"1dKw=/J%R#FX7$F_]o$\"1+\\p$)*eB Y#FX7$F_]o$\"1`l9RSvODFX7$F_]o$\"1;Ve))Rv1EFX7$F_]o$\"1*R:`5C*zEFX7$F_ ]o$\"1B(*po7Y^FFX7$F_]o$\"1kQ/\\SNEGFX7$F_]o$\"1dLt(p%[)*GFX7$F_]o$\"1 \"zySq[A(HFX7$F_]o$\"1T!4s!>SXIFX7$F_]o$\"1oz[.@i7JFX7$F_]o$\"1x)=U.k' *=$FX7$F_]o$\"1E\")[b7deKFX7$F_]o$\"1U=\\D2/KLFX7$F_]o$\"1P&yFQgBS$FX7 $F_]o$\"1+++e%GzZ$FX-F$6#7S7$$\"\")F(F(7$Fefo$\"1*e)*)o+)=A$F,7$Fefo$ \"1P\"[&*o?_-'F,7$Fefo$\"1_9r4/(y<*F,7$Fefo$\"1nz'QzT^B\"F07$Fefo$\"1] 6-:#))4b\"F07$Fefo$\"1Cjey#=Q%=F07$Fefo$\"1@)G![o-Z@F07$Fefo$\"1#f!H%e 01Y#F07$Fefo$\"1T#o&p'yJx#F07$Fefo$\"1gA*y])p%4$F07$Fefo$\"1_\\\\cL*yP $F07$Fefo$\"1=T]\"o.np$F07$Fefo$\"1P%QB2Bo,%F07$Fefo$\"1\")z3AuJDVF07$ Fefo$\"1Wu%f%)fag%F07$Fefo$\"1)prQGw&Q\\F07$Fefo$\"1S'[Y'F07$Fefo$\"1y&Q$4^pynF07$Fefo$\"1?OCICcpqF07$Fefo$\"1r)\\18-LQ (F07$Fefo$\"1JbV/&*=4xF07$Fefo$\"1r1K@a(G*zF07$Fefo$\"1yiQ$)RE*H)F07$F efo$\"1>C9[Lz:')F07$Fefo$\"1n/mybXD*)F07$Fefo$\"18>A4(p]A*F07$Fefo$\"1 U&*R=,ud&*F07$Fefo$\"1D)frFfm&)*F07$Fefo$\"1N`@G>e<5FX7$Fefo$\"1J\"=$o D]Y5FX7$Fefo$\"1qD*oY?\"y5FX7$Fefo$\"1k1TX/(y5\"FX7$Fefo$\"1P(G]un*Q6F X7$Fefo$\"1vUZR5Pp6FX7$Fefo$\"1r!G6Z+7?\"FX7$Fefo$\"1l$Go*f&=B\"FX7$Fe fo$\"1(fK&*p0KE\"FX7$Fefo$\"1:.(3\"eH%H\"FX7$Fefo$\"1CcH%RkGK\"FX7$Fef o$\"19S'[@2cN\"FX7$Fefo$\"1z.\"Ry#*[Q\"FX7$Fefo$\"1^$p'3t6;9FX7$Fefo$ \"1uK-jJ+Y9FX7$Fefo$\"1+++&f>\"y9FX-F$6#7S7$$\"\"*F(F(7$F[`p$\"1&=bJpa r@\"F,7$F[`p$\"1K(yN[%>wAF,7$F[`p$\"1PSU*[&>nMF,7$F[`p$\"1mI/B74mYF,7$ F[`p$\"11t^9))GfeF,7$F[`p$\"12tt@d`lpF,7$F[`p$\"1F07$F[`p$ \"1WBh%G*H'4#F07$F[`p$\"1)yb(fy11AF07$F[`p$\"1592%Q,lK#F07$F[`p$\"1yuE zD=TCF07$F[`p$\"1zrgv.%3c#F07$F[`p$\"10%odeB2n#F07$F[`p$\"1aLjsuC*y#F0 7$F[`p$\"1\"[a0[gB\"HF07$F[`p$\"1=F'erI&>IF07$F[`p$\"1^,-1vFNJF07$F[`p $\"1YV3m_&[D$F07$F[`p$\"1\"RNaxQ=P$F07$F[`p$\"12D\"fLE][$F07$F[`p$\"1q :(Q#=q5OF07$F[`p$\"1DfWRoiBPF07$F[`p$\"11N$pR)>WQF07$F[`p$\"1kUZ[TX`RF 07$F[`p$\"1rR8V&**G2%F07$F[`p$\"1\"[64N)G&=%F07$F[`p$\"1_>SQqw-VF07$F[ `p$\"1fHdGRi[%oy`%F07$F[`p$\"1w)**>@yOl%F07$F[`p$\"1 T1EA/6sZF07$F[`p$\"18`-a>c*)[F07$F[`p$\"1QGn8x[(*\\F07$F[`p$\"1HYk8G=@ ^F07$F[`p$\"1#RRc=<=B&F07$F[`p$\"12,W,lx\\`F07$F[`p$\"1U4$zhyEY&F07$F[ `p$\"1+++ht+%e&F0-F$6#7S7$$\"#5F(F(7$Faip$\"1RY6bfKQT!#?7$Faip$\"1Z7:T 71RxFfip7$Faip$\"1&pWfYY)y6F,7$Faip$\"1'G;_,rke\"F,7$Faip$\"1/%fh>e@*> F,7$Faip$\"1.wYW@GoBF,7$Faip$\"119*pqOxv#F,7$Faip$\"1MSns6^gJF,7$Faip$ \"1sq`NR*>c$F,7$Faip$\"1+ik0$o\\(RF,7$Faip$\"1b\\ibjrQVF,7$Faip$\"18(4 k]5#[ZF,7$Faip$\"1#)4_zgQf^F,7$Faip$\"1ZY>L*Hcb&F,7$Faip$\"1C*pJJda\"f F,7$Faip$\"1Y)47(oKVjF,7$Faip$\"1&)>h5&)y0nF,7$Faip$\"1GR!\\c,K=Dk9\"F07$Faip$\"17%[P&*3\\=\"F07$Faip$\"1Ryj>'QwA\"F07 $Faip$\"1b\"=\\7LgE\"F07$Faip$\"1azWau-28F07$Faip$\"1Lt&*4W]\"F07$Faip$\"1Ukyr_(Ga\"F07$Faip$\"1YKY\"f]Ae\"F07$Faip$\"17 i$Ha#e9*p\"F07$Faip$\"130' f:-7u\"F07$Faip$\"1mSUUy\")y " 0 "" {MPLTEXT 1 0 57 "sigsq:=simplify(sum((nu-mu)^2*PD(mu,nu),nu=0..infinity));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%&sigsqG$\"+,+++M!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 40 "Again, this can also be calculated from: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "simplify(sum(nu^2*PD(mu ,nu),nu=0..infinity)-mu^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"++++ +M!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 65 "This result may appear \+ to be puzzling: the deviation is given as " }{TEXT 19 8 "sqrt(mu)" } {TEXT -1 77 "! Indeed, if one lists the average with deviation as an i nterval one obtains " }{XPPEDIT 18 0 "[mu-sqrt(mu), mu+sqrt(mu)];" "6# 7$,&%#muG\"\"\"-%%sqrtG6#F%!\"\",&F%F&-F(6#F%F&" }{TEXT -1 128 " . Whi le using a larger count interval we do also obtain a larger uncertaint y, which seems improper. The fractional uncertainty " }{XPPEDIT 18 0 " sqrt(mu)/mu;" "6#*&-%%sqrtG6#%#muG\"\"\"F'!\"\"" }{TEXT -1 53 ", howev er, is reduced for an increased time interval." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "mu:=8.3;" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "P1:=plot([seq([nu,PD(mu,nu)*t,t=0.. 1],nu=0..20)],color=blue):" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#muG$ \"#$)!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "sigma:=sqrt(m u);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&sigmaG$\"+e?(4)G!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "P2:=plot(exp(-(nu-mu)^2/(2*s igma^2))/(sqrt(2*Pi)*sigma),nu=-0.1..20,color=red):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "display(P1,P2,axes=boxed);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6;-%'CURVESG6$7S7$\"\"!F(7$F($\"1W& 4LffpT&!#@7$F($\"1jA6rhbiYF07$F($\"16=4O(GJ?&F07$F($\"1&e'\\X`EzcF07$F($\" 1Qr.&R#G:iF07$F($\"1iD3#Q+Nv'F07$F($\"1T'o>9t@F(F07$F($\"19Oo(pxJu(F07 $F($\"1*pk4nZKI)F07$F($\"17$F($\"1_.KS3X'3\"Fao7$F($\"1U,wmYqR6Fao7$F($ \"1/tZ.$3')=\"Fao7$F($\"1/'GB[d8C\"Fao7$F($\"1n%oU0\\hH\"Fao7$F($\"1G= ry_%QM\"Fao7$F($\"1D9z2&e`R\"Fao7$F($\"1SP+9nd[9Fao7$F($\"1+RCY.k+:Fao 7$F($\"1Z$)\\(e95b\"Fao7$F($\"1#RQQfYpg\"Fao7$F($\"1+9@*)R?d;Fao7$F($ \"1S$)*f[k3r\"Fao7$F($\"1(H0;())[fFao7$F($\"1[c%R6og'>Fao7$F($\"11J>iIe>?F ao7$F($\"1!)yvi/h<#Fao7$F($\"1 GqF$4PTA#Fao7$F($\"1D\"))es(=zAFao7$F($\"1X)ezmD%GBFao7$F($\"11@(QeB4Q #Fao7$F($\"1HU)zyq6V#Fao7$F($\"1+++r#o^[#Fao-%'COLOURG6&%$RGBGF(F($\"* ++++\"!\")-F$6$7S7$$\"\"\"F(F(7$Fau$\"1-&*fik2'\\%F07$Fau$\"1?rJ_&)33% )F07$Fau$\"1twAi@*pQFao7$Fau$\"1gF5^of=VFao7$Fau$\"1M@#pB!z8ZFao7$Fau$\"1!*e0)Q%o e^Fao7$Fau$\"1a(=tJ0ag&Fao7$Fau$\"14)Q!3P!f.'Fao7$Fau$\"1.bH*[PoU'Fao7 $Fau$\"1zX8nbp\"*oFao7$Fau$\"1I()fM8\\&G(Fao7$Fau$\"18)f)G%pNu(Fao7$Fa u$\"1))oigh/\\\")Fao7$Fau$\"1F-)**R=Rf)Fao7$Fau$\"1S\\'\\(>aIWt[&f%*Fao7$Fau$\"1>k>>*[a')*Fao7$Fau$\"12xO5nKI5!#=7$Fau$\"1O $z]r.e2\"F[z7$Fau$\"1p'o8e\"R:6F[z7$Fau$\"1'=1:cZ\"e6F[z7$Fau$\"12Rms' =B?\"F[z7$Fau$\"11ZU([JbC\"F[z7$Fau$\"1;tm2@M(G\"F[z7$Fau$\"1G6$HnlPL \"F[z7$Fau$\"1T@46$zaP\"F[z7$Fau$\"17oUBv,?9F[z7$Fau$\"1c>[jdPg9F[z7$F au$\"1].()zz\\/:F[z7$Fau$\"1+r'pa8ga\"F[z7$Fau$\"1@,vw#4%*e\"F[z7$Fau$ \"1n-jMl$=j\"F[z7$Fau$\"1krnTRDw;F[z7$Fau$\"1ytz$=L!>F[z7$Fau$\"164oul;w>F[z7$Fau$\"1!R&4a<(y,#F[z7 $Fau$\"1+++l'*oi?F[zFft-F$6$7S7$$\"\"#F(F(7$Fh^l$\"1'GL?tre'=Fao7$Fh^l $\"17NJ\\nN*[$Fao7$Fh^l$\"1sDCYa8:`Fao7$Fh^l$\"116`*e;I:(Fao7$Fh^l$\"1 XGc/E;#)*)Fao7$Fh^l$\"1J&[2^,y1\"F[z7$Fh^l$\"1*4lY@(RV7F[z7$Fh^l$\"1jy %Ga)*\\U\"F[z7$Fh^l$\"10!oKZF[z7$Fh^l$\"1=d5@S&39#F[z7$Fh^l$\"1@_v1KCEBF[z7$Fh^l$\"1h#4*Q +!\\]#F[z7$Fh^l$\"1Gb8ev8nEF[z7$Fh^l$\"1yVWgO0gGF[z7$Fh^l$\"1un%R!*yM- $F[z7$Fh^l$\"1v1288e8KF[z7$Fh^l$\"1E*en:a=Q$F[z7$Fh^l$\"1`f4OhZmNF[z7$ Fh^l$\"1W(>(>*zAu$F[z7$Fh^l$\"1%=.\\s7d#RF[z7$Fh^l$\"1%3'e,8;%4%F[z7$F h^l$\"1f3:[e&eF%F[z7$Fh^l$\"1M'4xT&ekWF[z7$Fh^l$\"1l^\"G1v)GYF[z7$Fh^l $\"1R5R!Q7j![F[z7$Fh^l$\"1O4q\"\\A'*)\\F[z7$Fh^l$\"1!\\jHnb*o^F[z7$Fh^ l$\"1*QvpupCM&F[z7$Fh^l$\"1Ne#GaF^`&F[z7$Fh^l$\"1%4*>T\"R#3dF[z7$Fh^l$ \"1zLRAF2$*eF[z7$Fh^l$\"1=ri=%f01'F[z7$Fh^l$\"1*HWmhmOC'F[z7$Fh^l$\"1V 35?i&fT'F[z7$Fh^l$\"1Tcq)\\[gf'F[z7$Fh^l$\"1(Q8*=;7snF[z7$Fh^l$\"1'Rj \"eQXcpF[z7$Fh^l$\"1nu'Hr()R8(F[z7$Fh^l$\"1LuG@F[z7$F^hl$\"1eu:()\\1&[#F[z7$F^hl$\" 1H&zi%3DaHF[z7$F^hl$\"1`[4Fc1SMF[z7$F^hl$\"1'Gz$of\\URF[z7$F^hl$\"1*>( =4]JVWF[z7$F^hl$\"1kK9^*o%e\\F[z7$F^hl$\"1-BiBo@7aF[z7$F^hl$\"1$RU\"e% HI#fF[z7$F^hl$\"1$p<%Q&RfV'F[z7$F^hl$\"1n'*oSMAIpF[z7$F^hl$\"1=n>x03zt F[z7$F^hl$\"1N`*o7:G\"zF[z7$F^hl$\"1/ic+$e\\O)F[z7$F^hl$\"1E,7*H34*))F [z7$F^hl$\"1hzjmJYc$*F[z7$F^hl$\"1p@yU1Dn)*F[z7$F^hl$\"1X(y56k`.\"!#<7 $F^hl$\"1?S%Q&Q6'3\"F_\\m7$F^hl$\"1bw0gzrK6F_\\m7$F^hl$\"1@^g%y')H=\"F _\\m7$F^hl$\"1xT^l>?N7F_\\m7$F^hl$\"1a9*RVb1G\"F_\\m7$F^hl$\"1t`=DkuH8 F_\\m7$F^hl$\"1PckAAY!Q\"F_\\m7$F^hl$\"1#=Ehtx+V\"F_\\m7$F^hl$\"1!y+mH $3y9F_\\m7$F^hl$\"1?Qy'G&QJ:F_\\m7$F^hl$\"1m]m&\\z#z:F_\\m7$F^hl$\"1#* *fky;/j\"F_\\m7$F^hl$\"1OgTsZvw;F_\\m7$F^hl$\"1%)>`IWTF3A\"3vF_\\m7$F^hl$\"1EV]+mtt>F_\\m7$F^hl$\"1jS!GamR-#F_\\m7$F^h l$\"125#eb!yt?F_\\m7$F^hl$\"1A.IgWb>@F_\\m7$F^hl$\"1)3;/_;?<#F_\\m7$F^ hl$\"1(Q)Q7#R*=AF_\\m7$F^hl$\"1F**z)eo*oAF_\\m7$F^hl$\"1g.j&3`oJ#F_\\m 7$F^hl$\"1+++<:JoBF_\\mFft-F$6$7S7$$\"\"%F(F(7$Feam$\"1!)f1je;r5F[z7$F eam$\"1Wp-`:=.?F[z7$Feam$\"1Eh'GnI80$F[z7$Feam$\"1Q+X0wU1TF[z7$Feam$\" 1?#*[Y)4l:&F[z7$Feam$\"1ue\"Q]q+8'F[z7$Feam$\"1iNkth8QrF[z7$Feam$\"1+A VJ\"z1=)F[z7$Feam$\"1mKtS'y)>#*F[z7$Feam$\"1+1*pI#))G5F_\\m7$Feam$\"1$ *f\\!*\\.B6F_\\m7$Feam$\"1G177'G!H7F_\\m7$Feam$\"16$RPXdaL\"F_\\m7$Fea m$\"1+PUj8-Q9F_\\m7$Feam$\"1s_@>#f6`\"F_\\m7$Feam$\"1LxBQ\"4>k\"F_\\m7 $Feam$\"1^0uY)Gdt\"F_\\m7$Feam$\"1r5SrM'[%=F_\\m7$Feam$\"1HhP1hYT>F_\\ m7$Feam$\"1qxj2XXZ?F_\\m7$Feam$\"17CrH0Q[@F_\\m7$Feam$\"1BB\"fC'o`AF_ \\m7$Feam$\"1'oqiw*Q]BF_\\m7$Feam$\"1,!prd(paCF_\\m7$Feam$\"1fh^xS/jDF _\\m7$Feam$\"1VCdC+OdEF_\\m7$Feam$\"1LIEJQAfFF_\\m7$Feam$\"1p]&46fW'GF _\\m7$Feam$\"18+R,8TnHF_\\m7$Feam$\"1E)Q$RG-nIF_\\m7$Feam$\"1%oG*oWixJ F_\\m7$Feam$\"1u8K_\\+xKF_\\m7$Feam$\"1?DoI[6$Q$F_\\m7$Feam$\"1$e1::m# zMF_\\m7$Feam$\"11#f?%\\Q%e$F_\\m7$Feam$\"1:.(=`$H$o$F_\\m7$Feam$\"1Xw V2:o'y$F_\\m7$Feam$\"1Ih:^@w()QF_\\m7$Feam$\"1#)GY%F_\\m7$Feam$\"1MVx^Y4.VF_\\m7$Feam$ \"11Y`.b2)R%F_\\m7$Feam$\"1,bByU$p]%F_\\m7$Feam$\"1cG/k))H/YF_\\m7$Fea m$\"1.$f+#)4\"3ZF_\\m7$Feam$\"1kl%4:qu![F_\\m7$Feam$\"1+++'RYU\"\\F_\\ mFft-F$6$7S7$$\"\"&F(F(7$F[[n$\"1t\\EL`8yZIC'yP[u()R$F_\\m7$F[[n$\"1VnZ!o6jc$F_\\m7$F[[n$\"1n^Up\">6u$F_\\m7$F [[n$\"1,?F8qk,RF_\\m7$F[[n$\"1R#>%zzzuSF_\\m7$F[[n$\"1EP0qJlaUF_\\m7$F [[n$\"1Dy@Uw@6WF_\\m7$F[[n$\"1F*y8c6.e%F_\\m7$F[[n$\"1TUsD@+bZF_\\m7$F [[n$\"1c:!>w-f#\\F_\\m7$F[[n$\"1(4]\\\"zD\"4&F_\\m7$F[[n$\"1!oG,ic[F&F _\\m7$F[[n$\"1LRhA#G)RaF_\\m7$F[[n$\"1R/v?1(fh&F_\\m7$F[[n$\"1P+P8=evd F_\\m7$F[[n$\"1jQu0!z+&fF_\\m7$F[[n$\"1RM)[EnU6'F_\\m7$F[[n$\"1Z2Q/8*e G'F_\\m7$F[[n$\"1Nx+t^o`kF_\\m7$F[[n$\"1#zLJ!)\\$HmF_\\m7$F[[n$\"1<<=? d`)z'F_\\m7$F[[n$\"1fmUcFbrpF_\\m7$F[[n$\"1wqDCr8VrF_\\m7$F[[n$\"1`,BQ `!3I(F_\\m7$F[[n$\"1>@H/4^\"[(F_\\m7$F[[n$\"1KYy;h8VwF_\\m7$F[[n$\"1/x #QIia\"yF_\\m7$F[[n$\"1MR:`/S!)zF_\\m7$F[[n$\"1+++?!\\w:)F_\\mFft-F$6$ 7S7$$\"\"'F(F(7$Fadn$\"1S_YxQvfCF[z7$Fadn$\"1L_,uD(**f%F[z7$Fadn$\"1!e rnLso+(F[z7$Fadn$\"1h6`$fE(H%*F[z7$Fadn$\"1[PpWl5%=\"F_\\m7$Fadn$\"1)R !pB&owS\"F_\\m7$Fadn$\"1Iis.S:R;F_\\m7$Fadn$\"1*oA$=mby=F_\\m7$Fadn$\" 1Qzpl9><@F_\\m7$Fadn$\"1dY0;NF_\\m7$Fadn$\"1zfK&oq.x$F_\\m7$Fadn$\"1[_L4?\"e)RF_\\m7$ Fadn$\"1;H#e[@kB%F_\\m7$Fadn$\"12hlOLDeWF_\\m7$Fadn$\"1>rZN!Q;q%F_\\m7 $Fadn$\"1@BACyRL\\F_\\m7$Fadn$\"1X&z0&[@v^F_\\m7$Fadn$\"1>8)yOysR&F_\\ m7$Fadn$\"1S/.rQ!oj&F_\\m7$Fadn$\"1+hQ^Ng&)eF_\\m7$Fadn$\"1Xs_/W=-hF_ \\m7$Fadn$\"1TgMfw4OjF_\\m7$Fadn$\"15B+SHvxlF_\\m7$Fadn$\"18Lp'[lT\"oF _\\m7$Fadn$\"1DNs\"y1H/(F_\\m7$Fadn$\"1\"RJ0*\\)oH(F_\\m7$Fadn$\"1NCG2 d4DvF_\\m7$Fadn$\"1Uj,e$f(oxF_\\m7$Fadn$\"1b07T[b*)zF_\\m7$Fadn$\"1cH) QiU4B)F_\\m7$Fadn$\"1xDn:(p!e%)F_\\m7$Fadn$\"1\"[*GqH[&p)F_\\m7$Fadn$ \"1z&p=#yfF*)F_\\m7$Fadn$\"1\"4*=g0gq\"*F_\\m7$Fadn$\"188vy5k/%*F_\\m7 $Fadn$\"1uA?78)Rk*F_\\m7$Fadn$\"1&)*zwoR8))*F_\\m7$Fadn$\"1vcOQZ%*45!# ;7$Fadn$\"1)po;MU\\.\"Fd\\o7$Fadn$\"1::1'\\+t0\"Fd\\o7$Fadn$\"1%pL>&*Q 63\"Fd\\o7$Fadn$\"1+=#fRbR5\"Fd\\o7$Fadn$\"1+++\"yu%G6Fd\\oFft-F$6$7S7 $$\"\"(F(F(7$Fh]o$\"19'p?Jll\"HF[z7$Fh]o$\"1/9J&>`UX&F[z7$Fh]o$\"1=s]G '[\"3$)F[z7$Fh]o$\"1UlJ`h4=6F_\\m7$Fh]o$\"1*3cf/7SS\"F_\\m7$Fh]o$\"1%y 9#yE4p;F_\\m7$Fh]o$\"18#\\)=$oN%>F_\\m7$Fh]o$\"1+[.>9VFAF_\\m7$Fh]o$\" 1G!em;%Q5DF_\\m7$Fh]o$\"1w@!p%eV,GF_\\m7$Fh]o$\"1_RGDVzdIF_\\m7$Fh]o$ \"1n8g`KRYLF_\\m7$Fh]o$\"1pk%o>xhj$F_\\m7$Fh]o$\"1`?#=sPa\"RF_\\m7$Fh] o$\"1'HO@.L!pTF_\\m7$Fh]o$\"1)[jVQ#eqWF_\\m7$Fh]o$\"1qN**RU.EZF_\\m7$F h]o$\"1%eAkZ&=B]F_\\m7$Fh]o$\"1)f,5n9iG&F_\\m7$Fh]o$\"12&R`Q*zubF_\\m7 $Fh]o$\"1b$HxF+'\\eF_\\m7$Fh]o$\"1#\\J/=Ej8'F_\\m7$Fh]o$\"1Y=Dl+j*R'F_ \\m7$Fh]o$\"1'))\\h()QOo'F_\\m7$Fh]o$\"1)y]9@W'ypF_\\m7$Fh]o$\"1B;0BsW NsF_\\m7$Fh]o$\"1C*G4^,G^(F_\\m7$Fh]o$\"1sFWHcL*z(F_\\m7$Fh]o$\"1TtRjL lz!)F_\\m7$Fh]o$\"1=xm%=v3N)F_\\m7$Fh]o$\"1PQ8v2-_')F_\\m7$Fh]o$\"1q)e @[8E#*)F_\\m7$Fh]o$\"1^\\w`*G:@*F_\\m7$Fh]o$\"1pBK7$HLZ*F_\\m7$Fh]o$\" 1O?waiaf(*F_\\m7$Fh]o$\"1')*=#4a)G+\"Fd\\o7$Fh]o$\"13G)z$e.J5Fd\\o7$Fh ]o$\"1i$Hv-e&e5Fd\\o7$Fh]o$\"1H3O47P(3\"Fd\\o7$Fh]o$\"1CC*zs@^6\"Fd\\o 7$Fh]o$\"1W%*QT1]V6Fd\\o7$Fh]o$\"1L`0jXkr6Fd\\o7$Fh]o$\"1a;:/f](>\"Fd \\o7$Fh]o$\"1#)))\\>f9F7Fd\\o7$Fh]o$\"1Kmc-jl`7Fd\\o7$Fh]o$\"1ddC+>#>G \"Fd\\o7$Fh]o$\"1wO$Qov*38Fd\\o7$Fh]o$\"1+++p'[!Q8Fd\\oFft-F$6$7S7$$\" \")F(F(7$F^go$\"1*\\\"36l$f-$F[z7$F^go$\"12*y(*o(yecF[z7$F^go$\"1PH`,U q>')F[z7$F^go$\"1'z#**eZ-g6F_\\m7$F^go$\"1Gwe(\\imX\"F_\\m7$F^go$\"1!p QNl$oJOaQy&F_\\m7 $F^go$\"1y;w(Gg*ogF_\\m7$F^go$\"1&R'z6%QkO'F_\\m7$F^go$\"1%H!G#>;'RmF_ \\m7$F^go$\"1Q%oNLvU$pF_\\m7$F^go$\"1KQ#*oLMSsF_\\m7$F^go$\"11HW$\\wn] (F_\\m7$F^go$\"1nu4n:`%z(F_\\m7$F^go$\"1Zo.:2\"=4)F_\\m7$F^go$\"1!Q'*R OSEQ)F_\\m7$F^go$\"1ima.I.k')F_\\m7$F^go$\"1mbj`:Zw*)F_\\m7$F^go$\"1Vj :([6sD*F_\\m7$F^go$\"1g(G*)Ghpb*F_\\m7$F^go$\"1g\\\"4;z&G)*F_\\m7$F^go $\"1Sx'Q#Hb75Fd\\o7$F^go$\"1:\"**>h$\\S5Fd\\o7$F^go$\"1S(R=o*pp5Fd\\o7 $F^go$\"1!\\#\\`RD)4\"Fd\\o7$F^go$\"1`2kHw9G6Fd\\o7$F^go$\"1=#>FzQp:\" Fd\\o7$F^go$\"1rnO:>Q'=\"Fd\\o7$F^go$\"1!H$4M7e:7Fd\\o7$F^go$\"1J!zaP7 CC\"Fd\\o7$F^go$\"1U)*R\"*Q;t7Fd\\o7$F^go$\"1**3$QTo1I\"Fd\\o7$F^go$\" 1V6UrR**H8Fd\\o7$F^go$\"1S=*=Fi!e8Fd\\o7$F^go$\"1+++%\\D#)Q\"Fd\\oFft- F$6$7S7$$\"\"*F(F(7$Fd`p$\"1+X]$*ee!z#F[z7$Fd`p$\"1w5Ee(f'=_F[z7$Fd`p$ \"1^z\\jFG\\zF[z7$Fd`p$\"1lD061!)p5F_\\m7$Fd`p$\"1<(Hb_mLM\"F_\\m7$Fd` p$\"1FY7Zr*pf\"F_\\m7$Fd`p$\"13s#\\'phf=F_\\m7$Fd`p$\"1n,=c(=78#F_\\m7 $Fd`p$\"1S#F_\\m7$Fd`p$\"15Ij\"QH/o#F_\\m7$Fd`p$\"1XF%zf9d#HF_ \\m7$Fd`p$\"1I?'elZ=?$F_\\m7$Fd`p$\"1]ukUX6zMF_\\m7$Fd`p$\"1q*pVc7ju$F _\\m7$Fd`p$\"1o>7mR&*))RF_\\m7$Fd`p$\"1IReh!yuF%F_\\m7$Fd`p$\"1?\"3i\" e*=_%F_\\m7$Fd`p$\"1[2/q<@1[F_\\m7$Fd`p$\"1?qmk&zy0&F_\\m7$Fd`p$\"1[\\5Lm&f*F_\\m7$Fd`p$\"1e_zw\" 3]')*F_\\m7$Fd`p$\"10m`,U$G,\"Fd\\o7$Fd`p$\"1pR)R\"GSS5Fd\\o7$Fd`p$\"1 #4yVmap1\"Fd\\o7$Fd`p$\"14[Pmx5%4\"Fd\\o7$Fd`p$\"1(oxp-O57\"Fd\\o7$Fd` p$\"1u2\"=I!yX6Fd\\o7$Fd`p$\"1v%fJ.ST<\"Fd\\o7$Fd`p$\"1kl&\\J0&*>\"Fd \\o7$Fd`p$\"1v[%e**\\lA\"Fd\\o7$Fd`p$\"19&yi?NCD\"Fd\\o7$Fd`p$\"1++++C D!G\"Fd\\oFft-F$6$7S7$$\"#5F(F(7$Fjip$\"1N(3;H'=;BF[z7$Fjip$\"1#pc$*f( [JVF[z7$Fjip$\"1)H.Pp/zf'F[z7$Fjip$\"1*GO<2X$z))F[z7$Fjip$\"1l'*=;U*\\ 6\"F_\\m7$Fjip$\"1TM6Kw]D8F_\\m7$Fjip$\"1#e*)33#[V:F_\\m7$Fjip$\"1R&Hw c6*oF_\\m7$Fjip$\"1)RbnQcZA#F_\\m7$Fjip$\"1yC H;JMGCF_\\m7$Fjip$\"1&[lVbLvl#F_\\m7$Fjip$\"1$Q(Rq]m()GF_\\m7$Fjip$\"1 vqUG%R%4JF_\\m7$Fjip$\"1L7)=zJ3J$F_\\m7$Fjip$\"1iY6\"z1.b$F_\\m7$Fjip$ \"1RFXFN<`PF_\\m7$Fjip$\"1?Q8pd:*)RF_\\m7$Fjip$\"1EOnQ+/)>%F_\\m7$Fjip $\"1^Xt)=>sU%F_\\m7$Fjip$\"1*\\Sd$>XXYF_\\m7$Fjip$\"1wE78\\:t[F_\\m7$F jip$\"15h`FtD#3&F_\\m7$Fjip$\"1wfBID!yI&F_\\m7$Fjip$\"1An^213UbF_\\m7$ Fjip$\"1#G#\\i.-YdF_\\m7$Fjip$\"1X!y\\)3GmfF_\\m7$Fjip$\"1sM@d:$Q>'F_ \\m7$Fjip$\"15aG'\\XkT'F_\\m7$Fjip$\"1$*H$*Gf$=j'F_\\m7$Fjip$\"1rYK\"G !*4(oF_\\m7$Fjip$\"1@\"oP5\")e3(F_\\m7$Fjip$\"1eTDCHK:tF_\\m7$Fjip$\"1 5UP:8BBvF_\\m7$Fjip$\"1o3n-,`]xF_\\m7$Fjip$\"1:$3xa+W'zF_\\m7$Fjip$\"1 i+u'ycz=)F_\\m7$Fjip$\"1@QvsQ_1%)F_\\m7$Fjip$\"1Rp1cLMN')F_\\m7$Fjip$ \"1k\"QVrBd&))F_\\m7$Fjip$\"154\"4Y%4\"3*F_\\m7$Fjip$\"1+[\"R-*f/$*F_ \\m7$Fjip$\"1B%H]]w*4&*F_\\m7$Fjip$\"1POAvAOX(*F_\\m7$Fjip$\"1\"[RT6%* e&**F_\\m7$Fjip$\"1Y7b'\\O!=5Fd\\o7$Fjip$\"1m6@@7_R5Fd\\o7$Fjip$\"1+++ #\\4E1\"Fd\\oFft-F$6$7S7$$\"#6F(F(7$F`cq$\"1sm)*QymZ)z'y;F_\\m7$F`cq$\" 1kaGh_HK=F_\\m7$F`cq$\"1._jk)H_+#F_\\m7$F`cq$\"16:3aX()y@F_\\m7$F`cq$ \"1Z:%zQ8iM#F_\\m7$F`cq$\"1WKV!3t\")\\#F_\\m7$F`cq$\"11?m_y')yEF_\\m7$ F`cq$\"1h?^E+%>$GF_\\m7$F`cq$\"10ZSa$***4IF_\\m7$F`cq$\"13mM&F_\\m7$F`cq$\"1O*oMmV(>bF_\\m7$F`cq$\"1* z>]**>mn&F_\\m7$F`cq$\"1DS(=EF\"[eF_\\m7$F`cq$\"1/+6MA]4gF_\\m7$F`cq$ \"1x7=yb=yhF_\\m7$F`cq$\"1H\">&zV5VjF_\\m7$F`cq$\"1I3k!3fd^'F_\\m7$F`c q$\"17!>$yg/#o'F_\\m7$F`cq$\"1*=33b)4_oF_\\m7$F`cq$\"1:v'H!Hu?qF_\\m7$ F`cq$\"1vq(ya4d<(F_\\m7$F`cq$\"1S.#4\")=LN(F_\\m7$F`cq$\"1G;3`Y<7vF_\\ m7$F`cq$\"1czt8\"[:o(F_\\m7$F`cq$\"1tE:F,mVyF_\\m7$F`cq$\"1+++D;(y,)F_ \\mFft-F$6$7S7$$\"#7F(F(7$Ff\\r$\"1IVi(e.)37F[z7$Ff\\r$\"1*fyd7x0E#F[z 7$Ff\\r$\"1)\\$*4Z1MW$F[z7$Ff\\r$\"1Y5\"pPwSj%F[z7$Ff\\r$\"1%)pgOs3>eF [z7$Ff\\r$\"1&)Q\\&*Ru#QEVR`0)F[z7$Ff\\r$\"1z*z5EU=B* F[z7$Ff\\r$\"10I+.sXS5F_\\m7$Ff\\r$\"1/2Uuo3h6F_\\m7$Ff\\r$\"16bmavLn7 F_\\m7$Ff\\r$\"1#*f@b1&pQ\"F_\\m7$Ff\\r$\"1C&\\')*[02:F_\\m7$Ff\\r$\"1 cG(>f(zA;F_\\m7$Ff\\r$\"1B'ee/.zs\"F_\\m7$Ff\\r$\"1`fa'f$)G&=F_\\m7$Ff \\r$\"1VRq;&e(e>F_\\m7$Ff\\r$\"1VOR=i\">3#F_\\m7$Ff\\r$\"1([eElJ4>#F_ \\m7$Ff\\r$\"1yQ1i)Q0J#F_\\m7$Ff\\r$\"1^*e+@LWU#F_\\m7$Ff\\r$\"11,*Q,q Ka#F_\\m7$Ff\\r$\"1iBk'H*R_EF_\\m7$Ff\\r$\"1#oW4(*4,x#F_\\m7$Ff\\r$\"1 lCcz#yB*GF_\\m7$Ff\\r$\"1&*Rk,K\"))*HF_\\m7$Ff\\r$\"1??PLew8JF_\\m7$Ff \\r$\"1z,2]J_KKF_\\m7$Ff\\r$\"1<4l:Sq[LF_\\m7$Ff\\r$\"1)34q(\\6hMF_\\m 7$Ff\\r$\"1(3`S=Gfe$F_\\m7$Ff\\r$\"1kP(=Sy!)p$F_\\m7$Ff\\r$\"1D?Q(pAy \"QF_\\m7$Ff\\r$\"18UhE)Gj#RF_\\m7$Ff\\r$\"1b%*3'oa\\/%F_\\m7$Ff\\r$\" 1`i7vBdcTF_\\m7$Ff\\r$\"1:\\M@]CtUF_\\m7$Ff\\r$\"1)zT)eQJ(Q%F_\\m7$Ff \\r$\"1!Q^.OLn]%F_\\m7$Ff\\r$\"1O&\\O`[j\\F_\\m7$Ff\\r$\"11\\ k1b/'3&F_\\m7$Ff\\r$\"1(\\\"Q(z?f>&F_\\m7$Ff\\r$\"1cxBu228`F_\\m7$Ff\\ r$\"1B%*3]#)>DaF_\\m7$Ff\\r$\"1+++OXpXbF_\\mFft-F$6$7S7$$\"#8F(F(7$F\\ fr$\"1Mr1%)fu#F[z7$F \\fr$\"1eAaL&z'eHF[z7$F\\fr$\"1y&o-Rj_r$F[z7$F\\fr$\"1x)G\\Z8nT%F[z7$F \\fr$\"1W+-3W-V^F[z7$F\\fr$\"1e9-9i<%*eF[z7$F\\fr$\"1=Ylm!>Hk'F[z7$F\\ fr$\"1#4$yAJ48uF[z7$F\\fr$\"1Gbj]GY\"4)F[z7$F\\fr$\"1i*F[z7$F\\fr$\"1c(*QZQ4O5F_\\m7$F\\fr$\"1v<>P'*>.6F_\\ m7$F\\fr$\"1Et(=\"\\*H=\"F_\\m7$F\\fr$\"1+%*Hw>f]7F_\\m7$F\\fr$\"15$4? VB#H8F_\\m7$F\\fr$\"11+!RfD))R\"F_\\m7$F\\fr$\"1EL%o>!>v9F_\\m7$F\\fr$ \"1p3]dt!za\"F_\\m7$F\\fr$\"1Df6'3!yB;F_\\m7$F\\fr$\"1hq)e$\\X$p\"F_\\ m7$F\\fr$\"1'[5Tn3'oF_\\m7$F\\fr$\"1E5yC(>!))>F_\\m7$F\\fr$\"1^)=RqTQ1#F_\\m7$F\\fr$\"1[@ 6=(=!Q@F_\\m7$F\\fr$\"1r^X$z)y4AF_\\m7$F\\fr$\"1H2mDsZ*G#F_\\m7$F\\fr$ \"1FOGE33hBF_\\m7$F\\fr$\"106V`H`PCF_\\m7$F\\fr$\"1TqkV*4o]#F_\\m7$F\\ fr$\"1kSS#*za#e#F_\\m7$F\\fr$\"1k-)Qf6Ql#F_\\m7$F\\fr$\"1(H'4\"f-$GFF_ \\m7$F\\fr$\"1L'*ec68,GF_\\m7$F\\fr$\"1\"*4E:fPxGF_\\m7$F\\fr$\"1Hgwz' 33&HF_\\m7$F\\fr$\"17;g(G/f-$F_\\m7$F\\fr$\"1i$*)e1y.5$F_\\m7$F\\fr$\" 1+%)eL<\")oJF_\\m7$F\\fr$\"1Qp]RWCZKF_\\m7$F\\fr$\"1_No'e&R#R$F_\\m7$F\\fr$\"1bK!>U!yjMF_\\m7$F\\fr$\"1+++mErSNF_ \\mFft-F$6$7S7$$\"#9F(F(7$Fb_s$\"1cPvL3_vXFao7$Fb_s$\"1Yf-Yqlc&)Fao7$F b_s$\"1Zk24hQ.8F[z7$Fb_s$\"1yq,&HuSv\"F[z7$Fb_s$\"1g7oQ/i-AF[z7$Fb_s$ \"1)*4wJ,[=EF[z7$Fb_s$\"1o:Tc(y!\\IF[z7$Fb_s$\"1')G+qWS%\\$F[z7$Fb_s$ \"1>HQ=?IQRF[z7$Fb_s$\"1Ch$)3_!\\R%F[z7$Fb_s$\"1.;_Z:3(z%F[z7$Fb_s$\"1 3Zn\\o$)\\_F[z7$Fb_s$\"1)G)f37X/dF[z7$Fb_s$\"1#R89QcD9'F[z7$Fb_s$\"1R^ b\\yRSlF[z7$Fb_s$\"1AjAN)pM,(F[z7$Fb_s$\"1LV/gQA9uF[z7$Fb_s$\"1\")zyhg R!)yF[z7$Fb_s$\"1f3<]u.$H)F[z7$Fb_s$\"1>(oLXqdu)F[z7$Fb_s$\"1Q@WM$zo<* F[z7$Fb_s$\"1&y]bA'pE'*F[z7$Fb_s$\"1([A\\&o(R+\"F_\\m7$Fb_s$\"1'*RW&GK &[5F_\\m7$Fb_s$\"1A-obG\"[4\"F_\\m7$Fb_s$\"1U47G-5N6F_\\m7$Fb_s$\"1ab1 %p6'y6F_\\m7$Fb_s$\"1Wj)fHjNA\"F_\\m7$Fb_s$\"1[**oo'RvE\"F_\\m7$Fb_s$ \"1SO]q!*358F_\\m7$Fb_s$\"1HLuEHLd8F_\\m7$Fb_s$\"1#3IGj$y*R\"F_\\m7$Fb _s$\"1![*\\A)3^W\"F_\\m7$Fb_s$\"1;x)4=!='[\"F_\\m7$Fb_s$\"1UF=C?3J:F_ \\m7$Fb_s$\"1%zY$f6Lt:F_\\m7$Fb_s$\"1x\"G\"zO\\<;F_\\m7$Fb_s$\"1t\")Q9 1ng;F_\\m7$Fb_s$\"1!HvT'G(eq\"F_\\m7$Fb_s$\"12^8.!3%\\F_\\m7$Fb_s$\"1P`Rp;um>F_\\m7$Fb_s$\"1)=\"pm[36?F_\\m7$Fb _s$\"1,yckp_`?F_\\m7$Fb_s$\"1+++&zO\"*4#F_\\mFft-F$6$7S7$$\"#:F(F(7$Fh hs$\"1]]*[>)yJDFao7$Fhhs$\"1=`(ec$oMZFao7$Fhhs$\"1bsF[z7$Fhhs$\"1JGp$y$>z@F[z7$Fhhs$\"1 *o)4\\v%=V#F[z7$Fhhs$\"1![!>B^QaEF[z7$Fhhs$\"1nntB(4\\!HF[z7$Fhhs$\"10 GO-IYcJF[z7$Fhhs$\"1L!3!)>\")))R$F[z7$Fhhs$\"1'3KP9?!>OF[z7$Fhhs$\"1xE \\Umy!)QF[z7$Fhhs$\"1o5\"H?PD5%F[z7$Fhhs$\"1F9B@e[gVF[z7$Fhhs$\"1yi:cS \"))e%F[z7$Fhhs$\"1Gt%zJE$R[F[z7$Fhhs$\"1C=;KK(y2&F[z7$Fhhs$\"1Q+$>\"> xE`F[z7$Fhhs$\"1Ugr(fQ`b&F[z7$Fhhs$\"1_HiY'y=!eF[z7$Fhhs$\"1&y='oW'z0' F[z7$Fhhs$\"1FmZHz)3G'F[z7$Fhhs$\"1t3Ru8l@lF[z7$Fhhs$\"1EtqrNQqnF[z7$F hhs$\"1=#)o+&>P,(F[z7$Fhhs$\"1$exS_f\"\\sF[z7$Fhhs$\"1ojiGbd5vF[z7$Fhh s$\"1d*fcVpau(F[z7$Fhhs$\"1&))=&=)oi*zF[z7$Fhhs$\"1$p1bLINA)F[z7$Fhhs$ \"1Bs?,s)>Z)F[z7$Fhhs$\"12^8\\dw0()F[z7$Fhhs$\"1lk%>pJ,&*)F[z7$Fhhs$\" 1XB2nS/*=*F[z7$Fhhs$\"1%f#e#=j\"R%*F[z7$Fhhs$\"1))4==w0!o*F[z7$Fhhs$\" 11YtrpSE**F[z7$Fhhs$\"1ph#Gkrq,\"F_\\m7$Fhhs$\"1\\6eW6_R5F_\\m7$Fhhs$ \"1Gt,23Dl5F_\\m7$Fhhs$\"1!yDQsj#)3\"F_\\m7$Fhhs$\"1j#)*HH+G6\"F_\\m7$ Fhhs$\"15?\\?\\GO6F_\\m7$Fhhs$\"1+++gB_h6F_\\mFft-F$6$7S7$$\"#;F(F(7$F ^bt$\"1$z5i7lLJ\"Fao7$F^bt$\"1N;++r6cCFao7$F^bt$\"14uA67ETPFao7$F^bt$ \"1\"*)R=4C\\.&Fao7$F^bt$\"13.g=%QCK'Fao7$F^bt$\"1&48G,Hh^(Fao7$F^bt$ \"18G'=[E@v)Fao7$F^bt$\"1omC$oRI+\"F[z7$F^bt$\"1kWQvnXI6F[z7$F^bt$\"1H '*>;4_h7F[z7$F^bt$\"1D$pFGipP\"F[z7$F^bt$\"1&\\8$>>#p]\"F[z7$F^bt$\"1 \"e6E=:uj\"F[z7$F^bt$\"1zHh@'pJw\"F[z7$F^bt$\"12or*pmt(=F[z7$F^bt$\"1K ,'42eJ,#F[z7$F^bt$\"1=n=u6>G@F[z7$F^bt$\"14:&*>?+iAF[z7$F^bt$\"1.'G=HZ /Q#F[z7$F^bt$\"1&\\Vw_+/^#F[z7$F^bt$\"1LDSRn9MEF[z7$F^bt$\"13IaHHEjFF[ z7$F^bt$\"1QV!y*=$=)GF[z7$F^bt$\"15a=[Xs4IF[z7$F^bt$\"1snNo!pD9$F[z7$F ^bt$\"1ujo81@eKF[z7$F^bt$\"1.gu]m5$Q$F[z7$F^bt$\"1?Y6yk87NF[z7$F^bt$\" 1LRG\"*pOQOF[z7$F^bt$\"1,tU.:]gPF[z7$F^bt$\"1O#)zb56'*QF[z7$F^bt$\"1A? 3qA'z,%F[z7$F^bt$\"1')3&*\\W1[TF[z7$F^bt$\"1&4)*=OcfE%F[z7$F^bt$\"1S4* 4OV[R%F[z7$F^bt$\"1%fQF)f6;XF[z7$F^bt$\"10De:3)Gk%F[z7$F^bt$\"1R$*=5m \"ow%F[z7$F^bt$\"1[5P,fc'*[F[z7$F^bt$\"1&3[)*))H:-&F[z7$F^bt$\"1HY,iz$f^#R&F[z7$F^bt$\"1Z4'=c) )f_&F[z7$F^bt$\"1'p9y1o`k&F[z7$F^bt$\"1:vS?:lsdF[z7$F^bt$\"1S9aD!yW*eF [z7$F^bt$\"1+++VsRDgF[zFft-F$6$7S7$$\"#\"Fao7$Fd[u$\"1ZU9ndhE=Fao7$Fd[u$\"1m/\"RkF#eCFao7$Fd[ u$\"1g%z'zv$o3$Fao7$Fd[u$\"1e,\"yXR'pOFao7$Fd[u$\"1cATXq4tUFao7$Fd[u$ \"1Z`wXP>(*[Fao7$Fd[u$\"1k`3g*)G>bFao7$Fd[u$\"1@_B`.>fhFao7$Fd[u$\"10s i[e\"Gs'Fao7$Fd[u$\"1oI5iYKdtFao7$Fd[u$\"1;%H(H)QW*zFao7$Fd[u$\"1_=i%Q ;%3')Fao7$Fd[u$\"1J^cCo'f;*Fao7$Fd[u$\"1534D#[*G)*Fao7$Fd[u$\"1a\\ZV!e !R5F[z7$Fd[u$\"1o$)3@#*Q/6F[z7$Fd[u$\"1r#=KQ=A;\"F[z7$Fd[u$\"1;R]RhmD7 F[z7$Fd[u$\"1-_d2M3'G\"F[z7$Fd[u$\"1k*>')[A\"\\8F[z7$Fd[u$\"1#Q!pV?,29 F[z7$Fd[u$\"1s5W\\PXp9F[z7$Fd[u$\"1VVzKJJM:F[z7$Fd[u$\"1C#pmTt2f\"F[z7 $Fd[u$\"1:T.q?v^;F[z7$Fd[u$\"1tbs&)*[Zr\"F[z7$Fd[u$\"1BhUx\"zjx\"F[z7$ Fd[u$\"1?='ob4g$=F[z7$Fd[u$\"1&yBl#*=A!>F[z7$Fd[u$\"1Qd%**)4rh>F[z7$Fd [u$\"1b`(fYJ_-#F[z7$Fd[u$\"1RsDe/z#3#F[z7$Fd[u$\"1%\\V,kFao7$Fjdu$\"1xwy\\/:eAFao7$Fjdu$ \"1GVKla+XDFao7$Fjdu$\"1Frw262SGFao7$Fjdu$\"1B$Qi3l**4$Fao7$Fjdu$\"1/b F]Ta#R$Fao7$Fjdu$\"1T+1sXK'o$Fao7$Fjdu$\"1OyjWkVpRFao7$Fjdu$\"1\"eSv9H lA%Fao7$Fjdu$\"1b0$yBPA`%Fao7$Fjdu$\"1$Q-,477z%Fao7$Fjdu$\"1_W949Y#4&F ao7$Fjdu$\"1&RH75=\"f`Fao7$Fjdu$\"1N(e_v#o^cFao7$Fjdu$\"1x&39\\t-$fFao 7$Fjdu$\"1[N\\)p`4A'Fao7$Fjdu$\"1xv8p)))y['Fao7$Fjdu$\"1DhO7X\"ex'Fao7 $Fjdu$\"1_FlN$))[2(Fao7$Fjdu$\"1<_+yHBNtFao7$Fjdu$\"1@/7CBT;wFao7$Fjdu $\"1Ui2-\")*o!zFao7$Fjdu$\"1#y]Qw\"3\">)Fao7$Fjdu$\"1*Gp!pS/m%)Fao7$Fj du$\"1[ZAHGMr()Fao7$Fjdu$\"1Q`$RMnc/*Fao7$Fjdu$\"1bTy;tcQ$*Fao7$Fjdu$ \"1<3oD#yRg*Fao7$Fjdu$\"1$f7LXUT*)*Fao7$Fjdu$\"11([yV;n,\"F[z7$Fjdu$\" 1F')z&4b_/\"F[z7$Fjdu$\"1GkL]o:t5F[z7$Fjdu$\"1xu`,uO-6F[z7$Fjdu$\"1\\/ 5M1]I6F[z7$Fjdu$\"1,*ox#4Ff6F[z7$Fjdu$\"1q*z)))H!y=\"F[z7$Fjdu$\"12e\" )*)4-97F[z7$Fjdu$\"1l\")ex(pSC\"F[z7$Fjdu$\"1hferd%4F\"F[z7$Fjdu$\"1\" zsn=,'*H\"F[z7$Fjdu$\"1z^Xv!GqK\"F[z7$Fjdu$\"1+++8?]c8F[zFft-F$6$7S7$$ \"#>F(F(7$F`^v$\"1/0;['\\;H\"F07$F`^v$\"1?=h+r]:CF07$F`^v$\"1**\\ZYBSz OF07$F`^v$\"1>@X5bn^\\F07$F`^v$\"1$fm@y,z@'F07$F`^v$\"1i5HCc&=R(F07$F` ^v$\"1D1vHoT2')F07$F`^v$\"1z^Jn>bk)*F07$F`^v$\"13]lXbw66Fao7$F`^v$\"1- E7UEmS7Fao7$F`^v$\"169&[F&>a8Fao7$F`^v$\"1BJFao7$F`^v$\"1IJswJ+$4#Fao7$F`^v$\"1)G&\\j:gCAFao7$F`^v$\"1`rFM%)3T BFao7$F`^v$\"1,2q$*H*)oCFao7$F`^v$\"1wuTdJf!f#Fao7$F`^v$\"1&Qq5^uvr#Fa o7$F`^v$\"1-B,LI=MGFao7$F`^v$\"1!z!>x5'*fHFao7$F`^v$\"18p\"e6414$Fao7$ F`^v$\"1JW]<'QV?$Fao7$F`^v$\"1i%[]fprK$Fao7$F`^v$\"1Y7(Q,mSX$Fao7$F`^v $\"1.J*HS4#yNFao7$F`^v$\"1QEI:XK)p$Fao7$F`^v$\"1\\jXo=pJQFao7$F`^v$\"1 >w,n$G:&RFao7$F`^v$\"1Q-#ySz%zSFao7$F`^v$\"1)H^!o?U&>%Fao7$F`^v$\"1XN: ?!y@K%Fao7$F`^v$\"1?%fF\\X9W%Fao7$F`^v$\"16F*pO9hc%Fao7$F`^v$\"148&*y4 +)o%Fao7$F`^v$\"1y'3e'\\g:[Fao7$F`^v$\"13fCrF]Q\\Fao7$F`^v$\"1A)Gy_$=k ]Fao7$F`^v$\"1&Q#H%eB))=&Fao7$F`^v$\"1%*oc^[N.`Fao7$F`^v$\"1k.r31iMaFa o7$F`^v$\"1L*HnEE?b&Fao7$F`^v$\"19^t&=0sn&Fao7$F`^v$\"1LKp$QxDfFaoFft-F$6$7S7$$\"#?F(F(7$Ffgv$\"1MryQgMg`F,7$Ffgv$\"1I @bYaV-5F07$Ffgv$\"1r(zM(>&p_\"F07$Ffgv$\"1+SUO`%\\0#F07$Ffgv$\"1M(z!R# H/e#F07$Ffgv$\"1CLY$3?w1$F07$Ffgv$\"1.vi$)z2sNF07$Ffgv$\"1!y6c1*y$4%F0 7$Ffgv$\"1(eRN,FQh%F07$Ffgv$\"1JCxj*\\([^F07$Ffgv$\"1^@]*Q5*>cF07$Ffgv $\"1#\\I[VD.:'F07$Ffgv$\"1\"='H:%=Ho'F07$Ffgv$\"13V0R/<'>(F07$Ffgv$\"1 ,$\\qW_Am(F07$Ffgv$\"1=B(\\.pk@)F07$Ffgv$\"1NEl\"=jfo)F07$Ffgv$\"1frk' ['4K#*F07$Ffgv$\"1%f%H+q^:(*F07$Ffgv$\"1c&z@f!fC5Fao7$Ffgv$\"1_=6gh4v5 Fao7$Ffgv$\"1-s'=U$zF6Fao7$Ffgv$\"1dLYef=w6Fao7$Ffgv$\"1ApyYQQG7Fao7$F fgv$\"1]e!Gy-EG\"Fao7$Ffgv$\"1qm*f_+)H8Fao7$Ffgv$\"1#=n;Qv2Q\"Fao7$Ffg v$\"1P!oaRPMV\"Fao7$Ffgv$\"1QK%>!p&\\[\"Fao7$Ffgv$\"1h;atY![`\"Fao7$Ff gv$\"16$*3D@:!f\"Fao7$Ffgv$\"1FA)>F%))R;Fao7$Ffgv$\"1EY!R&R)Hp\"Fao7$F fgv$\"1J3*y:+6u\"Fao7$Ffgv$\"1vE+yQq$z\"Fao7$Ffgv$\"1]T7z)*>V=Fao7$Ffg v$\"1_0#>YP\\*=Fao7$Ffgv$\"1*=teS?b%>Fao7$Ffgv$\"1P$3/1w%)*>Fao7$Ffgv$ \"1AXl\\'y%\\?Fao7$Ffgv$\"1Lfijhj,@Fao7$Ffgv$\"1+z/(yhL:#Fao7$Ffgv$\"1 as&HE#*3?#Fao7$Ffgv$\"1H3;_vObAFao7$Ffgv$\"1V\"H-!44/BFao7$Ffgv$\"1Oeg ^,/cBFao7$Ffgv$\"1Z%[P@idS#Fao7$Ffgv$\"1+++g`>fCFaoFft-F$6$7en7$$!1+++ ++++5Fd\\o$\"1#)*Rfj[S(>F[z7$$\"1++]()zA\"Q$Fd\\o$\"1D)=d*R?SIF[z7$$\" 1,]7L5J$>(Fd\\o$\"1!3'>2GBWVF[z7$$\"1+]76(R![6!#:$\"1)yonwobN'F[z7$$\" 1+](y2!fz:F]bw$\"1,W:JqX:\"*F[z7$$\"1+D1h$*34?F]bw$\"1)oTN+)Qw7F_\\m7$ $\"1+DJw,H2CF]bw$\"1&>ccng'4GF]bw$\"1C$yccOxE#F_ \\m7$$\"1+D\")f'>gC$F]bw$\"1`o# )3QF_\\m7$$\"1,+vU>G3TF]bw$\"10u/D#=\\![F_\\m7$$\"1+]PZ/Q$\\%F]bw$\"1$ R!zeIb%y&F_\\m7$$\"1,+D`,\"p#\\F]bw$\"1r9f&=_v(pF_\\m7$$\"1++vt*>AO&F] bw$\"1AvH$)\\LL#)F_\\m7$$\"1,+vY;s\"y&F]bw$\"1XZ5cJa]%*F_\\m7$$\"1+D\" [NpE;'F]bw$\"1sta@]h^5Fd\\o7$$\"1++vDIl:mF]bw$\"1YZozp?n6Fd\\o7$$\"1,+ +8'*Q**pF]bw$\"1w,\\N1f]7Fd\\o7$$\"1,DJ*[mdW(F]bw$\"1%)4e.c>D8Fd\\o7$$ \"1]i!\\2ELk(F]bw$\"1'4R$[-C\\8Fd\\o7$$\"1++]gc)3%yF]bw$\"1(*H+klFn8Fd \\o7$$\"1CJ&*[FE\\zF]bw$\"1w!z_#\\_u8Fd\\o7$$\"1\\iSP)Rw0)F]bw$\"1_KQ\"Fd\\o7$$\"1,DJ9SRu#)F]b w$\"1]1xwTp%Q\"Fd\\o7$$\"1\\(o9j%fx$)F]bw$\"1fWA/nC%Q\"Fd\\o7$$\"1+]i[ _z![)F]bw$\"1u9zb[-#Q\"Fd\\o7$$\"1^7yle*Re)F]bw$\"1fMLtr.y8Fd\\o7$$\"1 -v$H['>(o)F]bw$\"1'zdW'*)Hs8Fd\\o7$$\"1_7.=Xb-*)F]bw$\"1@&*4G2za8Fd\\o 7$$\"1+]7`D\"z6*F]bw$\"1]C>RH0I8Fd\\o7$$\"1,DcgVW8&*F]bw$\"1T^Z.%4sE\" Fd\\o7$$\"1+]P\\)y+%**F]bw$\"1)=YKx+w<\"Fd\\o7$$\"1]7)z@B$Q5!#9$\"1\"H h_Zzh1\"Fd\\o7$$\"1](oc%)**o2\"Fd[x$\"1T=!o;+#\\F_\\m7$$\"1 +v`bdp*G\"Fd[x$\"1>iZfE6xQF_\\m7$$\"1++0aQMI8Fd[x$\"1k)4;tx[1$F_\\m7$$ \"1+]PvVut8Fd[x$\"1>wn=HrKBF_\\m7$$\"1]7G)prIT\"Fd[x$\"1\"p&4ZW<'y\"F_ \\m7$$\"1+]#*\\o1c9Fd[x$\"1lc$z*4\"fI\"F_\\m7$$\"1]7$p$>_'\\\"Fd[x$\"1 p$R()H!HI&*F[z7$$\"1]iXf\"4)Q:Fd[x$\"1\\)=YE)\\8nF[z7$$\"1+]()[E:!e\"F d[x$\"1z[ZXY6oYF[z7$$\"1]7y]aVB;Fd[x$\"1`v@/kR@JF[z7$$\"1++IL>7l;Fd[x$ \"1=9h^Elt?F[z7$$\"1+ve#Q_xq\"Fd[x$\"1Rrb0RmN8F[z7$$\"1]7BI)H+v\"Fd[x$ \"1%yH6EY)\\%)Fao7$$\"1++&QTy))y\"Fd[x$\"1XA'RE_NW&Fao7$$\"1+vtuKSL=Fd [x$\"1$G11]\"*e@$Fao7$$\"1++5ToAt=Fd[x$\"1IF#\\E$>o>Fao7$$\"1]7$*oqo:> Fd[x$\"1z?AT@wT6Fao7$$\"1](oZ#pKc>Fd[x$\"1*\\s&p'[Lk'F07$Ffgv$\"1&G0#z [(3j$F0-Fgt6&FitFjtF(F(-%+AXESLABELSG6$%!GF]cx-%*AXESSTYLEG6#%$BOXG-%% VIEWG6$;$!\"\"FgcxFfgv%(DEFAULTG" 1 2 0 1 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 299 "Note th at there is a systematic tendency for the Gaussian distribution to und erestimate the count probability to the left of the average. The Poiss on distribution is asymmetric, and this asymmetry becomes negligible o nly for sufficiently large average count probability for the chosen ti me interval." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 389 "To properly understand the Poisson distribution one should ana lyze data runs from a Geiger counter. These are usually computer-inter faced in undergraduate laboratories today, which has the drawback that the software performs all the data analysis, and the students usually no longer manipulate the raw data. For this reason we use Maple to ge nerate event sequences which are then analyzed." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 357 "A sequence of events, co rresponds to the sequence of clicks that one hears as the Geiger count er picks up the radiation (either from natural background, i.e., mostl y cosmic radiation, or from a small radioactive source, such as a cera mic plate or tile that has been coated with a radiactive layer contain ing, e.g., uranium oxide, to enhance the brightness)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 150 "We use a random numbe r generator to produce a time sequence of events; we assume that a tic k represents a second, and generate events at a given rate." }}{PARA 0 "" 0 "" {TEXT -1 224 "A random number generator for integers in the \+ interval [0,100] is used in conjunction with a given rate to decide in a time unit whether an observable decay occurs or not. A decay is rec orded as a one, a non-decay as a zero." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "rnd:=rand(0..100);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%$rndGR6\"6#%\"tGF&F&C%>%&_seedG-%%iremG6$,$F+\"-\"3p'>uU\"-*)****** ****>8$F+-F-6$F3\"$,\"F&6#F+F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "rate:=1/10;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%rateG#\"\"\" \"#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "tmax:=3000;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%tmaxG\"%+I" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 112 "TS:=[]: for i from 1 to tmax do: if rnd() <= 10 0*rate then TS:=[op(TS),[i,1]]: else TS:=[op(TS),[i,0]]: fi: od:" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 36 "We graph a subset of the count dat a:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "plot([seq(TS[i],i=1.. 200)],style=point);" }}{PARA 13 "" 1 "" {GLPLOT2D 691 299 299 {PLOTDATA 2 "6&-%'CURVESG6$7dw7$$\"\"\"\"\"!F*7$$\"\"#F*F*7$$\"\"$F*F* 7$$\"\"%F*F*7$$\"\"&F*F*7$$\"\"'F*F*7$$\"\"(F*F*7$$\"\")F*F*7$$\"\"*F* F*7$$\"#5F*F*7$$\"#6F*F*7$$\"#7F*F*7$$\"#8F*F*7$$\"#9F*F*7$$\"#:F*F(7$ $\"#;F*F*7$$\"#F*F*7$$\"#?F*F*7$$\"#@F*F*7$$\"# AF*F*7$$\"#BF*F*7$$\"#CF*F*7$$\"#DF*F(7$$\"#EF*F*7$$\"#FF*F*7$$\"#GF*F *7$$\"#HF*F*7$$\"#IF*F*7$$\"#JF*F*7$$\"#KF*F*7$$\"#LF*F*7$$\"#MF*F*7$$ \"#NF*F*7$$\"#OF*F*7$$\"#PF*F*7$$\"#QF*F*7$$\"#RF*F*7$$\"#SF*F*7$$\"#T F*F*7$$\"#UF*F*7$$\"#VF*F*7$$\"#WF*F(7$$\"#XF*F*7$$\"#YF*F*7$$\"#ZF*F* 7$$\"#[F*F*7$$\"#\\F*F*7$$\"#]F*F*7$$\"#^F*F*7$$\"#_F*F(7$$\"#`F*F*7$$ \"#aF*F*7$$\"#bF*F*7$$\"#cF*F*7$$\"#dF*F*7$$\"#eF*F(7$$\"#fF*F*7$$\"#g F*F*7$$\"#hF*F*7$$\"#iF*F*7$$\"#jF*F*7$$\"#kF*F*7$$\"#lF*F*7$$\"#mF*F( 7$$\"#nF*F*7$$\"#oF*F*7$$\"#pF*F*7$$\"#qF*F*7$$\"#rF*F(7$$\"#sF*F*7$$ \"#tF*F*7$$\"#uF*F*7$$\"#vF*F(7$$\"#wF*F(7$$\"#xF*F*7$$\"#yF*F*7$$\"#z F*F*7$$\"#!)F*F*7$$\"#\")F*F*7$$\"##)F*F(7$$\"#$)F*F*7$$\"#%)F*F*7$$\" #&)F*F*7$$\"#')F*F*7$$\"#()F*F(7$$\"#))F*F*7$$\"#*)F*F*7$$\"#!*F*F*7$$ \"#\"*F*F*7$$\"##*F*F*7$$\"#$*F*F*7$$\"#%*F*F*7$$\"#&*F*F(7$$\"#'*F*F* 7$$\"#(*F*F*7$$\"#)*F*F*7$$\"#**F*F(7$$\"$+\"F*F*7$$\"$,\"F*F*7$$\"$- \"F*F(7$$\"$.\"F*F*7$$\"$/\"F*F*7$$\"$0\"F*F*7$$\"$1\"F*F*7$$\"$2\"F*F *7$$\"$3\"F*F*7$$\"$4\"F*F*7$$\"$5\"F*F*7$$\"$6\"F*F*7$$\"$7\"F*F*7$$ \"$8\"F*F*7$$\"$9\"F*F*7$$\"$:\"F*F*7$$\"$;\"F*F*7$$\"$<\"F*F*7$$\"$= \"F*F*7$$\"$>\"F*F*7$$\"$?\"F*F*7$$\"$@\"F*F*7$$\"$A\"F*F(7$$\"$B\"F*F *7$$\"$C\"F*F*7$$\"$D\"F*F*7$$\"$E\"F*F(7$$\"$F\"F*F*7$$\"$G\"F*F*7$$ \"$H\"F*F*7$$\"$I\"F*F*7$$\"$J\"F*F*7$$\"$K\"F*F*7$$\"$L\"F*F*7$$\"$M \"F*F*7$$\"$N\"F*F*7$$\"$O\"F*F*7$$\"$P\"F*F*7$$\"$Q\"F*F*7$$\"$R\"F*F *7$$\"$S\"F*F*7$$\"$T\"F*F(7$$\"$U\"F*F*7$$\"$V\"F*F*7$$\"$W\"F*F*7$$ \"$X\"F*F(7$$\"$Y\"F*F*7$$\"$Z\"F*F*7$$\"$[\"F*F*7$$\"$\\\"F*F*7$$\"$] \"F*F*7$$\"$^\"F*F*7$$\"$_\"F*F(7$$\"$`\"F*F*7$$\"$a\"F*F*7$$\"$b\"F*F *7$$\"$c\"F*F*7$$\"$d\"F*F(7$$\"$e\"F*F*7$$\"$f\"F*F(7$$\"$g\"F*F*7$$ \"$h\"F*F*7$$\"$i\"F*F*7$$\"$j\"F*F*7$$\"$k\"F*F(7$$\"$l\"F*F*7$$\"$m \"F*F*7$$\"$n\"F*F(7$$\"$o\"F*F*7$$\"$p\"F*F*7$$\"$q\"F*F*7$$\"$r\"F*F *7$$\"$s\"F*F*7$$\"$t\"F*F*7$$\"$u\"F*F*7$$\"$v\"F*F*7$$\"$w\"F*F*7$$ \"$x\"F*F*7$$\"$y\"F*F*7$$\"$z\"F*F*7$$\"$!=F*F*7$$\"$\"=F*F*7$$\"$#=F *F*7$$\"$$=F*F*7$$\"$%=F*F*7$$\"$&=F*F*7$$\"$'=F*F*7$$\"$(=F*F(7$$\"$) =F*F*7$$\"$*=F*F*7$$\"$!>F*F*7$$\"$\">F*F*7$$\"$#>F*F*7$$\"$$>F*F*7$$ \"$%>F*F*7$$\"$&>F*F*7$$\"$'>F*F*7$$\"$(>F*F*7$$\"$)>F*F*7$$\"$*>F*F*7 $$\"$+#F*F*-%'COLOURG6&%$RGBG$FE!\"\"F*F*-%+AXESLABELSG6$%!GFcam-%&STY LEG6#%&POINTG-%%VIEWG6$%(DEFAULTGF[bm" 1 5 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 302 "This se quence could have come from a Geiger counter. Note the irregularity: t here are stretches with no counts followed by segments where the count s appear to cluster. This irregularity is noticable when listening to \+ the clicks from a Geiger counter. Now we pick a time interval for the \+ decay analysis:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "tint:=25 ; # length of interval in seconds (time units)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%tintG\"#D" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 73 "We can determine the expected average number of events per time interval :" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "nu_bar_exp:=tint*rate; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%+nu_bar_expG#\"\"&\"\"#" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "Nt:=tmax/tint; # number of t ime intervals contained in sample" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %#NtG\"$?\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 123 "Now we need to ge nerate the histogram (frequency distribution): for each time interval \+ we find out how many events occured." }}{PARA 0 "" 0 "" {TEXT -1 157 " We begin by resetting the counters for the frequency. Assume that 10 e vents during the interval is the maximum (increase it, if the program \+ below complains):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "Emax:=1 5;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%EmaxG\"#:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "for icnt from 0 to Emax do: C[icnt]:=0; od: myrate:=0:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 240 "for i from 1 to Nt \+ do: icnt:=0: for j from 1 to tint do: if TS[(i-1)*tint+j][2]=1 then ic nt:=icnt+1; myrate:=myrate+1; fi; od: if icnt <= Emax then C[icnt]:=C[ icnt]+1; else print(\"Increase Emax, please\"); fi; od: myrate:=evalf( myrate/tmax);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'myrateG$\"++++!4\" !#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "mymu:=myrate*tint;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%mymuG$\"++++DF!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 138 "To compare the found results with the Po isson distribution we have to convert the probabilities calculated fro m the Poissonian to numbers." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 72 "P2:=plot([seq([nu-0.05,Nt*PD(mymu,nu)*t,t=0..1],nu=0..Emax)],col or=red):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 97 "P1:=plot([seq([icnt+0.0 5,C[icnt]*t,t=0..1],icnt=0..Emax)],color=blue): display(P1,P2,axes=box ed);" }{TEXT -1 0 "" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6E-%'CURVESG6$7S7$$\"1+++++++]!#<\"\"!7$F($\"1LLLLBxV5E$F/7$F($\"1LLLLAKn\\F/7$F($\"1LLLLc$\\o'F/7$F($\"1mmm; bQ%R)F/7$F($\"1LLL$Qk#z**F/7$F($\"1+++l9.i6!#:7$F($\"1LLL=\"\\qM`FB7$F($\"1++++.W2bFB7$F($\"1LL Lep'Rm&FB7$F($\"1+++S>4NeFB7$F($\"1mmm6s5'*fFB7$F($\"1+++lXTkhFB7$F($ \"1mmmmd'*GjFB7$F($\"1+++DcB,lFB7$F($\"1MLLt>:nmFB7$F($\"1LLL.a#o$oFB7 $F($\"1nmm^Q40qFB7$F($\"1+++!3:(frFB7$F($\"1nmmc%GpL(FB7$F($\"1LLL8-V& \\(FB7$F($\"1+++XhUkwFB7$F($\"1+++:o5FB7$Fcu$\"1nm\"z>)G_:FB7$Fcu$\"1nmT&QU!*3#FB7$Fcu$\"1L$eR ZXKi#FB7$Fcu$\"1n;z>,_=JFB7$Fcu$\"1+]7G$[8j$FB7$Fcu$\"1m;z%*frhTFB7$Fc u$\"1**\\ilFQ!p%FB7$Fcu$\"1LL$3_\"=M_FB7$Fcu$\"1nmTg(fJr&FB7$Fcu$\"1++ ]7eP_iFB7$Fcu$\"1++]Pf!Qz'FB7$Fcu$\"1++](=ubJ(FB7$Fcu$\"1m;zW(*Q*y(FB7 $Fcu$\"1LL$3F-GN)FB7$Fcu$\"1LLL$e'3I))FB7$Fcu$\"1**\\7.\"Fay7$Fcu$\"1n;/Egw[7Fay7$Fcu$\"1m\"z%*f%)Q I\"Fay7$Fcu$\"1+voza'=N\"Fay7$Fcu$\"1n;zWho.9Fay7$Fcu$\"1++]i>Ad9Fay7$ Fcu$\"1+]i:jf4:Fay7$Fcu$\"1+DJ&>r-c\"Fay7$Fcu$\"1+]P4q`;;Fay7$Fcu$\"1L L$eM%4n;Fay7$Fcu$\"1++v$4v5s\"Fay7$Fcu$\"1n\"zWn*)*pFay7$Fcu$\"1LL3-=!y(>Fay7$ Fcu$\"1+D\"G8O;.#Fay7$Fcu$\"1nmm\"*\\[$3#Fay7$Fcu$\"1n;aQz]O@Fay7$Fcu$ \"1LekG=4*=#Fay7$Fcu$\"1++]i4TPAFay7$Fcu$\"1L$3F9!z#H#Fay7$Fcu$\"1nmmT >KUBFay7$Fcu$\"1+DJqJ8&R#Fay7$Fcu$\"1+voa-oXCFay7$Fcu$\"#DF+Fht-F$6$7S 7$$\"1++++++]?FBF+7$Fj^lFfu7$Fj^lFiu7$Fj^lF\\v7$Fj^lF_v7$Fj^lFbv7$Fj^l Fev7$Fj^lFhv7$Fj^lF[w7$Fj^lF^w7$Fj^lFaw7$Fj^lFdw7$Fj^lFgw7$Fj^lFjw7$Fj ^lF]x7$Fj^lF`x7$Fj^lFcx7$Fj^lFfx7$Fj^lFix7$Fj^lF\\y7$Fj^lF_y7$Fj^lFcy7 $Fj^lFfy7$Fj^lFiy7$Fj^lF\\z7$Fj^lF_z7$Fj^lFbz7$Fj^lFez7$Fj^lFhz7$Fj^lF [[l7$Fj^lF^[l7$Fj^lFa[l7$Fj^lFd[l7$Fj^lFg[l7$Fj^lFj[l7$Fj^lF]\\l7$Fj^l F`\\l7$Fj^lFc\\l7$Fj^lFf\\l7$Fj^lFi\\l7$Fj^lF\\]l7$Fj^lF_]l7$Fj^lFb]l7 $Fj^lFe]l7$Fj^lFh]l7$Fj^lF[^l7$Fj^lF^^l7$Fj^lFa^l7$Fj^lFd^lFht-F$6$7S7 $$\"1++++++]IFBF+7$F`bl$\"1MLLeaM8]F/7$F`bl$\"1nmTN0Vv$*F/7$F`bl$\"1LL 3U^5G9FB7$F`bl$\"1LLe%**=>#>FB7$F`bl$\"1m;/OeQ8CFB7$F`bl$\"1L$3-^Q!pGF B7$F`bl$\"1+](=YS3M$FB7$F`bl$\"1L$3_ry(GQFB7$F`bl$\"1**\\PW@::VFB7$F`b l$\"1nm;**pW:[FB7$F`bl$\"1MLezp5c_FB7$F`bl$\"1++]Zd=_dFB7$F`bl$\"1++]i 9I]iFB7$F`bl$\"1++]_#G.t'FB7$F`bl$\"1L$3_cQi;(FB7$F`bl$\"1mm;*3yXo(FB7 $F`bl$\"1lmmc'zO7)FB7$F`bl$\"1**\\(o;fWj)FB7$F`bl$\"1mmm\"e&e'3*FB7$F` bl$\"1,](o\"*REe*FB7$F`bl$\"1+D1&4+b+\"Fay7$F`bl$\"1n;a8gya5Fay7$F`bl$ \"1nT5se/+6Fay7$F`bl$\"1L$eRuk)[6Fay7$F`bl$\"1L3_JQd*>\"Fay7$F`bl$\"1+ DJTirV7Fay7$F`bl$\"1L$3KD\"R\"H\"Fay7$F`bl$\"1++]0UkS8Fay7$F`bl$\"1+]P 5'G))Q\"Fay7$F`bl$\"1+vo*\\\\aV\"Fay7$F`bl$\"1+]i[S@([\"Fay7$F`bl$\"1n m;)zEP`\"Fay7$F`bl$\"1++D'3*Q$e\"Fay7$F`bl$\"1L3_+0RG;Fay7$F`bl$\"1++v K*)ex;Fay7$F`bl$\"1mTNB3)Qs\"Fay7$F`bl$\"1+vV(=pAx\"Fay7$F`bl$\"1nm\"z lx&>=Fay7$F`bl$\"1+v=U_5p=Fay7$F`bl$\"1LLL#>1o\">Fay7$F`bl$\"1L$eMI(el >Fay7$F`bl$\"1nTN#[kR,#Fay7$F`bl$\"1++]&3=%e?Fay7$F`bl$\"1n;HJpO4@Fay7 $F`bl$\"1LLL'=O\\:#Fay7$F`bl$\"1+vo;D_.AFay7$F`bl$\"1+DJMe-]AFay7$F`bl $\"#BF+Fht-F$6$7S7$$\"1++++++]SFBF+7$Ff[m$\"1+++vBSxXF/7$Ff[m$\"1,+D1d EN?#FB7 $Ff[m$\"1+]i+pb>EFB7$Ff[m$\"1+]i&fK.0$FB7$Ff[m$\"1+]iN9%e\\$FB7$Ff[m$ \"1**\\7B:#*RRFB7$Ff[m$\"1++]xCr'R%FB7$Ff[m$\"1++v)>a!*z%FB7$Ff[m$\"1+ +]#o&*>D&FB7$Ff[m$\"1++]()pz1dFB7$Ff[m$\"1++]P%f#>Fay7$Ff[m$\"1+++J/bn>Fay7$Ff[m$\"1+D1j=\">,#Fay7$Ff[m$\"1+v $RTrV0#Fay7$Ff[m$\"#@F+Fht-F$6$7S7$$\"1++++++]]FBF+7$F\\em$\"1++++&ech #F/7$F\\em$\"1+++v*G:*[F/7$F\\em$\"1+++]L)4X(F/7$F\\em$\"1+++X.u-5FB7$ F\\em$\"1++]Fy:f7FB7$F\\em$\"1++]d'*)o\\\"FB7$F\\em$\"1++](>ZIu\"FB7$F \\em$\"1++]xOi(*>FB7$F\\em$\"1****\\FPQ^AFB7$F\\em$\"1+++IrS7DFB7$F\\e m$\"1+++&o;Bu#FB7$F\\em$\"1+++!RS6+$FB7$F\\em$\"1+++]o-hKFB7$F\\em$\"1 +++5cZ6NFB7$F\\em$\"1++]xq!*QPFB7$F\\em$\"1+++!4X$4SFB7$F\\em$\"1+++g: WQUFB7$F\\em$\"1++]<_$\\]%FB7$F\\em$\"1+++gs#3u%FB7$F\\em$\"1++]<#Q'** \\FB7$F\\em$\"1++]_u3Y_FB7$F\\em$\"1+++v8B.bFB7$F\\em$\"1****\\n(p$RdF B7$F\\em$\"1+++Dp2%*fFB7$F\\em$\"1****\\xgkeiFB7$F\\em$\"1++]-V&*)['FB 7$F\\em$\"1*****\\\\$pPnFB7$F\\em$\"1+++?am%*pFB7$F\\em$\"1*****\\JigC (FB7$F\\em$\"1****\\P(>%F/7$Fb^n$\"1mmm\">K'*)\\F/7$Fb^n$\"1*****\\Kd,\" eF/7$Fb^n$\"1mmm\"fX(emF/7$Fb^n$\"1*****\\U7Y](F/7$Fb^n$\"1MLLL/pu$)F/ 7$Fb^n$\"1nmm;c0T\"*F/7$Fb^n$\"1+++I,Q+5FB7$Fb^n$\"1+++]*3q3\"FB7$Fb^n $\"1+++q=\\q6FB7$Fb^n$\"1nm;fBIY7FB7$Fb^n$\"1LLLj$[kL\"FB7$Fb^n$\"1LLL `Q\"GT\"FB7$Fb^n$\"1++]s]k,:FB7$Fb^n$\"1LLL`dF!e\"FB7$Fb^n$\"1++]sgam; FB7$Fb^n$\"1++]FB7$Fb^n $\"1nmmTc-)*>FB7$Fb^n$\"1mm;f`@'3#FB7$Fb^n$\"1++]nZ)H;#FB7$Fb^n$\"1mmm Jy*eC#FB7$Fb^n$\"1+++S^bJBFB7$Fb^n$\"1+++0TN:CFB7$Fb^n$\"1++]7RV'\\#FB 7$Fb^n$\"1+++:#fke#FB7$Fb^n$\"1LLL`4NnEFB7$Fb^n$\"1+++],s`FFB7$Fb^n$\" 1mm;zM)>$GFB7$Fb^n$\"1+++qfaY2aF/7$Fhgn$\"1mm;zXu9cF/7$Fhgn$\"1+ ++]y))GeF/7$Fhgn$\"1****\\i_QQgF/7$Fhgn$\"1***\\7y%3TiF/7$Fhgn$\"1**** \\P![hY'F/7$Fhgn$\"1LLL$Qx$omF/7$Fhgn$\"1+++v.I%)oF/7$Fhgn$\"1mm\"zpe* zqF/7$Fhgn$\"1+++D\\'QH(F/7$Fhgn$\"1KLe9S8&\\(F/7$Fhgn$\"1***\\i?=bq(F /7$Fhgn$\"1LLL3s?6zF/7$Fhgn$\"1++DJXaE\")F/7$Fhgn$\"1nmmm*RRL)F/7$Fhgn $\"1mm;a<.Y&)F/7$Fhgn$\"1LLe9tOc()F/7$Fhgn$\"1+++]Qk\\*)F/7$Fhgn$\"1LL $3dg6<*F/7$Fhgn$\"1mmmmxGp$*F/7$Fhgn$\"1++D\"oK0e*F/7$Fhgn$\"1++v=5s#y *F/7$Fhgn$\"\"\"F+Fht-F$6$7S7$$\"1,+++++]!)FBF+7$F^aoF[hn7$F^aoF^hn7$F ^aoFahn7$F^aoFdhn7$F^aoFghn7$F^aoFjhn7$F^aoF]in7$F^aoF`in7$F^aoFcin7$F ^aoFfin7$F^aoFiin7$F^aoF\\jn7$F^aoF_jn7$F^aoFbjn7$F^aoFejn7$F^aoFhjn7$ F^aoF[[o7$F^aoF^[o7$F^aoFa[o7$F^aoFd[o7$F^aoFg[o7$F^aoFj[o7$F^aoF]\\o7 $F^aoF`\\o7$F^aoFc\\o7$F^aoFf\\o7$F^aoFi\\o7$F^aoF\\]o7$F^aoF_]o7$F^ao Fb]o7$F^aoFe]o7$F^aoFh]o7$F^aoF[^o7$F^aoF^^o7$F^aoFa^o7$F^aoFd^o7$F^ao Fg^o7$F^aoFj^o7$F^aoF]_o7$F^aoF`_o7$F^aoFc_o7$F^aoFf_o7$F^aoFi_o7$F^ao F\\`o7$F^aoF_`o7$F^aoFb`o7$F^aoFe`o7$F^aoFh`oFht-F$6$7S7$$\"1,+++++]!* FBF+FcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFc doFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdo FcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFcdoFht-F$6$7S7$$ \"1++++++05FayF+FidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFi doFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFido FidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFidoFh t-F$6$7S7$$\"1++++++06FayF+F_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_e oF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF _eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_eoF_e oF_eoF_eoFht-F$6$7S7$$\"1++++++07FayF+FeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeo FeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFe eoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeoFeeo FeeoFeeoFeeoFeeoFeeoFht-F$6$7S7$$\"1++++++08FayF+F[foF[foF[foF[foF[foF [foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[f oF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF[foF [foF[foF[foF[foF[foF[foF[foF[foFht-F$6$7S7$$\"1++++++09FayF+FafoFafoFa foFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafo FafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFa foFafoFafoFafoFafoFafoFafoFafoFafoFafoFafoFht-F$6$7S7$$\"1++++++0:FayF +FgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoF gfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgf oFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFgfoFht-F$6$7S7$$!1++ +++++]F*F+7$F^go$\"1BdyA[Y9#eslF/7$F^go$\"1f#f8U,LD)F/7$F^go$\"1tw\\$HV:\")*F /7$F^go$\"1F0_\\7]U6FB7$F^go$\"1\"Gq^Ym$48FB7$F^go$\"1'>d>a'pv9FB7$F^g o$\"1(y(\\quyY;FB7$F^go$\"1tI[`W[(z\"FB7$F^go$\"1g6%QcLr'>FB7$F^go$\"1 Hfup#zu8#FB7$F^go$\"12pn@*Q;I#FB7$F^go$\"1ZfS7;r]CFB7$F^go$\"1A@IgM(zi #FB7$F^go$\"1\\]%p\\P\"yFFB7$F^go$\"1eq\"*pR\"G&HFB7$F^go$\"1%f'GYAV2J FB7$F^go$\"1!\\U6+tqF$FB7$F^go$\"1k#4j66'QMFB7$F^go$\"1.Z-$>frg$FB7$F^ go$\"1JYYU*Q>w$FB7$F^go$\"1SMD7)*))GRFB7$F^go$\"1coFX[I-TFB7$F^go$\"18 &389jKD%FB7$F^go$\"1j.mOAI;WFB7$F^go$\"1%yXbdPZe%FB7$F^go$\"1))e'Rp=& \\ZFB7$F^go$\"1)eC*)*G&*3\\FB7$F^go$\"1:Mu+u(f3&FB7$F^go$\"1A/6*yT]C&F B7$F^go$\"1#41`(p([T&FB7$F^go$\"1z7Z'*GxobFB7$F^go$\"1R'FB7$F^go$\"1jR!=U(4blFB7$F^go$\"1X,Qd\">>s'FB7$F^go$\"1!GfZ_ft )oFB7$F^go$\"1I)3//#QRqFB7$F^go$\"1/`2**ph8sFB7$F^go$\"1@mn%F/7$Ff`p$\"1+YWKj#pt )F/7$Ff`p$\"1LM\\([X3L\"FB7$Ff`p$\"1A\"4wkG5z\"FB7$Ff`p$\"1KN\\PY-\\AF B7$Ff`p$\"1QRoZbktEFB7$Ff`p$\"10J=0fJ8JFB7$Ff`p$\"1!pjyhB!oNFB7$Ff`p$ \"1CH9xIF@SFB7$Ff`p$\"1D>Xde\\([%FB7$Ff`p$\"1V%=h8X\")*[FB7$Ff`p$\"1>) z='*Q/O&FB7$Ff`p$\"1\\g]55jCeFB7$Ff`p$\"1+C0ig'>F'FB7$Ff`p$\"11.KW\"*= ymFB7$Ff`p$\"1mX()zwAhrFB7$Ff`p$\"1V!e(zYUqvFB7$Ff`p$\"17Mko!=k/)FB7$F f`p$\"1IGyrGvn%)FB7$Ff`p$\"1M2!)GR-I*)FB7$Ff`p$\"1\"3kEH:-P*FB7$Ff`p$ \"12!\\OG^-\"Fay7$Ff`p$\"1\"[@RZA12\"Fay7$Ff`p $\"1+SUq!yy6\"Fay7$Ff`p$\"1$z&41U,f6Fay7$Ff`p$\"1#f#efBW.7Fay7$Ff`p$\" 1.CV*)4M\\7Fay7$Ff`p$\"1Q`?WQC%H\"Fay7$Ff`p$\"1G*f]m*oP8Fay7$Ff`p$\"1s #4=%)GfQ\"Fay7$Ff`p$\"1=Ek()QFH9Fay7$Ff`p$\"1,X)e(Rbv9Fay7$Ff`p$\"1V_S 91\\<:Fay7$Ff`p$\"1-iAP%QLc\"Fay7$Ff`p$\"1Vj\\lwZ1;Fay7$Ff`p$\"1fA/*eq :l\"Fay7$Ff`p$\"16W/krl&p\"Fay7$Ff`p$\"1T0MA<\"=u\"Fay7$Ff`p$\"1!4zv/k iy\"Fay7$Ff`p$\"1*4F0(HsJ=Fay7$Ff`p$\"1%QT1Z0o(=Fay7$Ff`p$\"1v!fhgJ#=> Fay7$Ff`p$\"1X$*RK1rl>Fay7$Ff`p$\"1A-ZRk<3?Fay7$Ff`p$\"1fX*p%QX`?Fay7$ Ff`p$\"1tBE.,z'4#Fay7$Ff`p$\"1+++H3OV@FayF``p-F$6$7S7$$\"1++++++]>FBF+ 7$F\\jp$\"1a\"FB7$F\\jp$\"1*\\7q(pF8=FB7$ F\\jp$\"1\"zuF`w-W#FB7$F\\jp$\"1qk4pgHkIFB7$F\\jp$\"1xuKf?%Gk$FB7$F\\j p$\"1%4?:#H*=C%FB7$F\\jp$\"173:!=K9'[FB7$F\\jp$\"1f=_oX)*yaFB7$F\\jp$ \"19*)z1J@9hFB7$F\\jp$\"1Xa2\\FstmFB7$F\\jp$\"1m/y'3)f.tFB7$F\\jp$\"1% GRJvfg$zFB7$F\\jp$\"1LH(4w`ba)FB7$F\\jp$\"1Z`l5F.*4*FB7$F\\jp$\"1$*yX: FFay7$F\\jp$\"1U#yCHU/,#Fay7$F\\jp$\"1\"HsC'4e n?Fay7$F\\jp$\"1wi2@'[+8#Fay7$F\\jp$\"1bTg>e#))=#Fay7$F\\jp$\"1Ad%z#\\ E]AFay7$F\\jp$\"1\"='R')GL5BFay7$F\\jp$\"1d^3s%=KP#Fay7$F\\jp$\"1*yK_w %yLCFay7$F\\jp$\"1#zfxaAd\\#Fay7$F\\jp$\"1edy&F/7$Fbcq$\"1xM%Qk&G\"3\"FB7$Fbcq$\"1D bN@*fqk\"FB7$Fbcq$\"1=k.=Xe;AFB7$Fbcq$\"1%[CPE-My#FB7$Fbcq$\"1\\#*pr[ \"*3LFB7$Fbcq$\"18*y$QF0`QFB7$Fbcq$\"1%G'\\)*3!eT%FB7$Fbcq$\"1jr7=Vuw \\FB7$Fbcq$\"19#GRdVPb&FB7$Fbcq$\"1@]]K['>1'FB7$Fbcq$\"1&z%y(f,Tj'FB7$ Fbcq$\"1yDf6we3sFB7$Fbcq$\"1ChcN8@ixFB7$Fbcq$\"1_wi))FB7$Fbcq$\"1\\k!>8*=p$*FB7$Fbcq$\"1y')fL0Ge**FB7$Fbcq$\"1%HzT isz/\"Fay7$Fbcq$\"1\"4QG[$=06Fay7$Fbcq$\"133Z\\`()oFay7$Fbcq$\"1hT\\`M=))>Fay7$Fbcq$\"1q.fV1*R/#Fay7$Fbcq$ \"1f]`sBb)4#Fay7$Fbcq$\"1%GAiWtc:#Fay7$Fbcq$\"1GP=Hyo5AFay7$Fbcq$\"1(> l;)z%pE#Fay7$Fbcq$\"1\")*y;FUFK#Fay7$Fbcq$\"1\\3jH9,uBFay7$Fbcq$\"1(>l .nrFV#Fay7$Fbcq$\"19R;IvK&[#Fay7$Fbcq$\"1;1SQIOTDFay7$Fbcq$\"1$)Gvph* \\f#Fay7$Fbcq$\"1+++%HKEl#FayF``p-F$6$7S7$$\"1++++++]RFBF+7$Fh\\r$\"1z >#GEw*QRF/7$Fh\\r$\"1a<^Z%eiO(F/7$Fh\\r$\"1Q\\Cr&f?7\"FB7$Fh\\r$\"1q?n -#[+^\"FB7$Fh\\r$\"1J1'=z#>'*=FB7$Fh\\r$\"1%[')fD)>aAFB7$Fh\\r$\"1HIgx @*[i#FB7$Fh\\r$\"1'*>ihVE3IFB7$Fh\\r$\"1NZ=mqS!R$FB7$Fh\\r$\"1km84y[$y $FB7$Fh\\r$\"1Kq[TNrHTFB7$Fh\\r$\"1o=x7<[>XFB7$Fh\\r$\"1chh4.&3\"\\FB7 $Fh\\r$\"1JJ^rk+)G&FB7$Fh\\r$\"1@;s)R+0j&FB7$Fh\\r$\"1*Q%*[!)ex.'FB7$F h\\r$\"12e4&4gFQ'FB7$Fh\\r$\"1\"zMChySy'FB7$Fh\\r$\"1?c3,NJRrFB7$Fh\\r $\"1nmjvC1HvFB7$Fh\\r$\"1,I,nc>+zFB7$Fh\\r$\"1bdy+$y7\"Fay7$Fh\\r$\"1Y<*)H8]o6Fay7$Fh\\r$\"1i \"f7JY]?\"Fay7$Fh\\r$\"1m*Q#ze1W7Fay7$Fh\\r$\"1snF2MUz7Fay7$Fh\\r$\"1# 3pX_y!=8Fay7$Fh\\r$\"17uZx*\\WN\"Fay7$Fh\\r$\"11WsD'oCR\"Fay7$Fh\\r$\" 1h?c.)Q'H9Fay7$Fh\\r$\"1$H0]`_&o9Fay7$Fh\\r$\"19=R$3Jg]\"Fay7$Fh\\r$\" 1r]UL5FB7$F^fr$\"1Tg2\\!Q&G7FB7$F^fr$\"1i_GOi cI9FB7$F^fr$\"1[8]wS]R;FB7$F^fr$\"1FST]=xZ=FB7$F^fr$\"1(pY_&3+i?FB7$F^ fr$\"1!f5$zQp]AFB7$F^fr$\"1O`eKv6jCFB7$F^fr$\"1s'*y)GFB7$F^fr$\"1pF@;FioIFB7$F^fr$\"1]!)\\Z&y0H$FB7$F^fr$\"1.5f]UgyM FB7$F^fr$\"1AIYUGK(p$FB7$F^fr$\"1A*3n&e#4*QFB7$F^fr$\"1cPwZ!RL5%FB7$F^ fr$\"1x>*pj1cI%FB7$F^fr$\"1F\")3Z;l;XFB7$F^fr$\"1#Ryu%oX5ZFB7$F^fr$\"1 PCFQ>]>\\FB7$F^fr$\"1.]bO5kO^FB7$F^fr$\"1H&H*f7mD`FB7$F^fr$\"1\\Abz#3) HbFB7$F^fr$\"1kFiE@rSdFB7$F^fr$\"1D@'\\RSq%fFB7$F^fr$\"1]By.OnYhFB7$F^ fr$\"1mokXALojFB7$F^fr$\"1,&GSR-vc'FB7$F^fr$\"1dGWR!f,y'FB7$F^fr$\"1A/ =nv&G(pFB7$F^fr$\"1W'[j&z_$=(FB7$F^fr$\"1,'GYP_) >q'oy\\(*F*7$Fd_s$\"1n'=)GCIB=F/7$Fd_s$\"1z$R-?Jtx#F/7$Fd_s$\"1c(RO*Qo PPF/7$Fd_s$\"1aeK)>sMp%F/7$Fd_s$\"1n\"*=1.hzbF/7$Fd_s$\"1_#f(*)\\:(\\' F/7$Fd_s$\"1q'oS,\"3YuF/7$Fd_s$\"1pMe?Q'>R)F/7$Fd_s$\"1R(f3b?\\O*F/7$F d_s$\"1Zy%o8!>A5FB7$Fd_s$\"1aiPzem=6FB7$Fd_s$\"1#os0wPb@\"FB7$Fd_s$\"1 D&G!o<*)38FB7$Fd_s$\"1w@F]hm$R\"FB7$Fd_s$\"1yOJ!4rW\\\"FB7$Fd_s$\"1%eM r4m)z:FB7$Fd_s$\"1qD%4z+#z;FB7$Fd_s$\"1e(yl#)Grw\"FB7$Fd_s$\"1Bc^<)*fj =FB7$Fd_s$\"1$f*G9IYb>FB7$Fd_s$\"1f5/BEJ^?FB7$Fd_s$\"1zW)[[K$R@FB7$Fd_ s$\"16H&pztUB#FB7$Fd_s$\"1eBUP6*GL#FB7$Fd_s$\"1vW?Ayt=CFB7$Fd_s$\"1Vfh $fa9^#FB7$Fd_s$\"1+m#34Ssg#FB7$Fd_s$\"1pop7v%4q#FB7$Fd_s$\"19AqhUh\"z# FB7$Fd_s$\"1>^@GUG#*GFB7$Fd_s$\"11Ygq1u#)HFB7$Fd_s$\"1?kRIAKzIFB7$Fd_s $\"1v+Ty%Ro;$FB7$Fd_s$\"1TM#o%*=DE$FB7$Fd_s$\"1&3D,(ea_LFB7$Fd_s$\"1WT 0X*\\mW$FB7$Fd_s$\"1!=c*yQlQNFB7$Fd_s$\"1%zq7Zt\\j$FB7$Fd_s$\"1j[_'oSx s$FB7$Fd_s$\"1Ezyi&3E#QFB7$Fd_s$\"1ux8&)3p;RFB7$Fd_s$\"1w_OoF9.SFB7$Fd _s$\"15Sj#QEA5%FB7$Fd_s$\"11.U[z%3>%FB7$Fd_s$\"1!RV:-P`G%FB7$Fd_s$\"1Z i0@^xvVFB7$Fd_s$\"1+++5K'HZ%FBF``p-F$6$7S7$$\"1++++++]pFBF+7$Fjhs$\"1_ ?FWFX&z$F*7$Fjhs$\"1@w)G)e&y4(F*7$Fjhs$\"14B)*f`<\"3\"F/7$Fjhs$\"1X0;W p-b9F/7$Fjhs$\"1Z!=qn,r#=F/7$Fjhs$\"1J!eKei?<#F/7$Fjhs$\"1TSv?'\\#HDF/ 7$Fjhs$\"1`P*e'Hl)*GF/7$Fjhs$\"1\"o?)H;(oE$F/7$Fjhs$\"1Hl\\&yHck$F/7$F jhs$\"1x\"y^t%F/7$Fjhs$ \"1>U(ymG`4&F/7$Fjhs$\"1u))GFKMDaF/7$Fjhs$\"1wW7eCw]hF/7$Fjhs$\"1su+=zoF/7$Fjhs$\"1o*)HgysasF/ 7$Fjhs$\"1O&*)>;QBh(F/7$Fjhs$\"1x%>g*pY&)zF/7$Fjhs$\"1$zuQu:\"G$)F/7$F jhs$\"17Zn3(3xp)F/7$Fjhs$\"139*o#>h\"3*F/7$Fjhs$\"10=!p!3!eT*F/7$Fjhs$ \"1tlr_Rtw(*F/7$Fjhs$\"1'pn4Kh\\,\"FB7$Fjhs$\"1'f?;GS90\"FB7$Fjhs$\"1' [l\"FB7$Fjhs$\"1PYDo`!GB\"FB7$Fjhs$\"1#fa%=?0q7FB7$Fjhs$\"1BMf8#)40 8FB7$Fjhs$\"1Ach&eJc6X\"FB7$Fjhs$\"1ypk/p3)[\"FB7$Fjhs$\"1T3q3=rC:FB7$Fjhs $\"1`4\"=Hm$e:FB7$Fjhs$\"17;OF\"Qpf\"FB7$Fjhs$\"15j*oBP9j\"FB7$Fjhs$\" 1`a@$[?#o;FB7$Fjhs$\"1[$zLuEMq\"FB7$Fjhs$\"1+++y1ETF/7$F`bt$\"1TZS(\\4\\4#F/7$F`bt$\"1h$H[uFmA#F/7$F`bt$\"1j #)f72AVBF/7$F`bt$\"1e?Bx;9rCF/7$F`bt$\"1AkB'o_Hf#F/7$F`bt$\"1'\\=`p\\+ s#F/7$F`bt$\"1/\"[]Ukn$GF/7$F`bt$\"1!zt'GqliHF/7$F`bt$\"11#pB1CM4$F/7$ F`bt$\"1^Vh\\rD2KF/7$F`bt$\"1Q%)R-+?ILF/7$F`bt$\"1,7H**>@dMF/7$F`bt$\" 1vf+f%o9e$F/7$F`bt$\"1Q\"GH(Hp,PF/7$F`bt$\"12-#zz\"=NQF/7$F`bt$\"1A4LW u7bRF/7$F`bt$\"1l*fP-&>$3%F/7$F`bt$\"14wI&GV#*>%F/7$F`bt$\"1Z=a(o9hK%F /7$F`bt$\"1(oC(*y!\\XWF/7$F`bt$\"1'4mJ@t-d%F/7$F`bt$\"1M0BS3F#p%F/7$F` bt$\"1t`h\\5**>[F/7$F`bt$\"1K*f6z+I%\\F/7$F`bt$\"1!))z&=gzo]F/7$F`bt$ \"1]T+)f\\N>&F/7$F`bt$\"17&*3\"=&=3`F/7$F`bt$\"1_U#f\\q&RaF/7$F`bt$\"1 _\"ox3$3dbF/7$F`bt$\"1*[L0-wBo&F/7$F`bt$\"15z?WtH-eF/7$F`bt$\"1+++%=p6 $fF/F``p-F$6$7S7$$\"1************\\*)FBF+7$Ff[u$\"1_\\CpLF*7$Ff[u$\"1EGqy@()fPF*7$Ff[u$\"1]OV_o$R5%F*7$Ff [u$\"1B#eN=s7\\%F*7$Ff[u$\"11ow]z>!)[F*7$Ff[u$\"1cWm$****\\D&F*7$Ff[u$ \"1GHHbhN&f&F*7$Ff[u$\"1T>U^D2+gF*7$Ff[u$\"1$H!oD/#HM'F*7$Ff[u$\"1FzM3 StTnF*7$Ff[u$\"1m(p`g^Z4(F*7$Ff[u$\"1%*>efy1#[(F*7$Ff[u$\"1Cv!f_%)3&yF *7$Ff[u$\"1+#4e-1dB)F*7$Ff[u$\"1=+*4H#4*e)F*7$Ff[u$\"1&p8)=rEq*)F*7$Ff [u$\"1s()oP1?m$*F*7$Ff[u$\"1IvyH>'3r*F*7$Ff[u$\"1$)evc0J35F/7$Ff[u$\"1 vU%p'pwY5F/7$Ff[u$\"1R*3<1*Q%3\"F/7$Ff[u$\"1Ae)GP!z?6F/7$Ff[u$\"1$4%pA y?h6F/7$Ff[u$\"1Db)*)pCv>\"F/7$Ff[u$\"1j'3hr+jB\"F/7$Ff[u$\"1c:(3bP9F \"F/7$Ff[u$\"1_;7#R^)48F/7$Ff[u$\"1Tt2Be*fM\"F/7$Ff[u$\"1A/(z;xPQ\"F/7 $Ff[u$\"1CXSQ`r?9F/7$Ff[u$\"1ivfzhQf9F/7$Ff[u$\"1d\\!ezIm\\\"F/7$Ff[u$ \"1GhB&z=Z`\"F/7$Ff[u$\"1%z,oS\"\\s:F/7$Ff[u$\"1h*=R^+sg\"F/7$Ff[u$\"1 jffn5)pk\"F/7$Ff[u$\"1(4K!p7c#o\"F/7$Ff[u$\"1&)Q![=(\\?F[\\u7$F]eu$\"1jGlUu_QIF[ \\u7$F]eu$\"1#p7:C(>*3%F[\\u7$F]eu$\"1(e@-Rt[8&F[\\u7$F]eu$\"1()*Qf.\\ V5'F[\\u7$F]eu$\"1,%eG\"\\=3rF[\\u7$F]eu$\"1>0x7MNY\")F[\\u7$F]eu$\"1c p8$\\#>\"=*F[\\u7$F]eu$\"1e&=&o^cC5F*7$F]eu$\"1AR4#zA$=6F*7$F]eu$\"1Nq !)p;(QA\"F*7$F]eu$\"1lHNTR&)H8F*7$F]eu$\"1N[-)\\()>V\"F*7$F]eu$\"1%[_? bMZ_\"F*7$F]eu$\"1%[y*p(>]j\"F*7$F]eu$\"1dwp:eWGL$>F*7$F]eu$\"1_>Q\"\\j)Q?F*7$F]eu$\"1-#QH.m$R@F*7$F ]eu$\"1&>V;*)HUC#F*7$F]eu$\"1T#)))[w_SBF*7$F]eu$\"1^@X*y(RWCF*7$F]eu$ \"1qoKt'*G_DF*7$F]eu$\"1ouSv)4ik#F*7$F]eu$\"1GDomikZFF*7$F]eu$\"1&pFW) *RC&GF*7$F]eu$\"1'fVE>g\\&HF*7$F]eu$\"1tLelP:aIF*7$F]eu$\"1\"y\"H;8HkJ F*7$F]eu$\"1W(oT![DjKF*7$F]eu$\"1U4\"3X>*oLF*7$F]eu$\"1%>BbKnYY$F*7$F] eu$\"1C:\"H/X$pNF*7$F]eu$\"1p'[shQym$F*7$F]eu$\"1LW1#y#zqPF*7$F]eu$\"1 w)zkH\\9(QF*7$F]eu$\"1-ErQt#o(RF*7$F]eu$\"1@,'y\"*=$ySF*7$F]eu$\"1P?H; (3@=%F*7$F]eu$\"1)f!z#3R]G%F*7$F]eu$\"1\\%o'*R@'zVF*7$F]eu$\"1]%>%3M-) [%F*7$F]eu$\"1ugJd%z\\e%F*7$F]eu$\"1TixF[N)o%F*7$F]eu$\"1wmqZ#)H(y%F*7 $F]eu$\"1+++lki$*[F*F``p-F$6$7S7$$\"1++++++&4\"FayF+7$Fc^v$\"1&fG1tNCk #!#>7$Fc^v$\"1j9#3$\\gT\\Fh^v7$Fc^v$\"1T%Q!f6EFvFh^v7$Fc^v$\"1yY^zc+85 F[\\u7$Fc^v$\"12Ozz\"[?F\"F[\\u7$Fc^v$\"1S.IPP@7:F[\\u7$Fc^v$\"1il;h7* 3w\"F[\\u7$Fc^v$\"1$)e37R2=?F[\\u7$Fc^v$\"19?)\\xJWF#F[\\u7$Fc^v$\"1zu NZt7QDF[\\u7$Fc^v$\"1\"GQ5b!RqFF[\\u7$Fc^v$\"1*o#*p!R'=.$F[\\u7$Fc^v$ \"1\"\\ow*3T%H$F[\\u7$Fc^v$\"1r*e>dBua$F[\\u7$Fc^v$\"10@E>L=xPF[\\u7$F c^v$\"1wSk\"))*Q]SF[\\u7$Fc^v$\"1\"RXlnJ=G%F[\\u7$Fc^v$\"1@vX/N0^XF[\\ u7$Fc^v$\"1UH/S/O*y%F[\\u7$Fc^v$\"19w*z\"4#30&F[\\u7$Fc^v$\"1pu-BJz*H& F[\\u7$Fc^v$\"1'ow+]p&fbF[\\u7$Fc^v$\"1Tih2`7)z&F[\\u7$Fc^v$\"1Gy\"Q(* Ra0'F[\\u7$Fc^v$\"12=.&f&y**eo^)*F[\\u7$Fc^v$\"1Hk1A3J55F*7$Fc^v$\"1f$\\dQAg.\" F*7$Fc^v$\"10NqW5_h5F*7$Fc^v$\"1#4%4k;&\\3\"F*7$Fc^v$\"1Y(f!)z0=6\"F*7 $Fc^v$\"1#34cXCe8\"F*7$Fc^v$\"114=aLVh6F*7$Fc^v$\"1sftPV%f=\"F*7$Fc^v$ \"1+++PZG77F*F``p-F$6$7S7$$\"1++++++&>\"FayF+7$Fjgv$\"1U7DW6`+g!#?7$Fj gv$\"1d)>37c@7\"Fh^v7$Fjgv$\"1h&flb:$4l`)RV$Fh^v7$Fjgv$\"1Bm>w0p)*RFh^v7$Fj gv$\"1M#QK^4Fe%Fh^v7$Fjgv$\"1#*>9%\\b[;&Fh^v7$Fjgv$\"1e#Q.>kOw&Fh^v7$F jgv$\"1'H')e@&4\"H'Fh^v7$Fjgv$\"1VZ&*Gw&[)oFh^v7$Fjgv$\"1p(Qv3e5[(Fh^v 7$Fjgv$\"1)*Hk-xgb!)Fh^v7$Fjgv$\"1ZW*fz`td)Fh^v7$Fjgv$\"1**=Nl.w(>*Fh^ v7$Fjgv$\"1P8Y3hKB(*Fh^v7$Fjgv$\"1A^*4SoM.\"F[\\u7$Fjgv$\"1gzLMRe(3\"F [\\u7$Fjgv$\"1&pl7Ddp9\"F[\\u7$Fjgv$\"1vsa*o%\\.7F[\\u7$Fjgv$\"1^tAvb[ i7F[\\u7$Fjgv$\"1jJCBwl;8F[\\u7$Fjgv$\"1,xT&\\*3v8F[\\u7$Fjgv$\"1Z&G') y$yN9F[\\u7$Fjgv$\"17/H\\$=')[\"F[\\u7$Fjgv$\"18`DA6oX:F[\\u7$Fjgv$\"1 4WK]Dj/;F[\\u7$Fjgv$\"1mdj@]Ii;F[\\u7$Fjgv$\"1FF4.h5=F[\\u7$Fjgv$\"1qM=!fDz+#F[\\u7$Fjgv$\"1([#G*)GLj?F[\\u7$F jgv$\"14&G<[\\77#F[\\u7$Fjgv$\"1)yLq`ty<#F[\\u7$Fjgv$\"1\"4%HvO:PAF[\\ u7$Fjgv$\"1?.0%\\ZUH#F[\\u7$Fjgv$\"1B:(*pTj_BF[\\u7$Fjgv$\"1-ttut`5CF[ \\u7$Fjgv$\"1SGK1WujCF[\\u7$Fjgv$\"1-.$*ocsCDF[\\u7$Fjgv$\"1C$pr.o#zDF [\\u7$Fjgv$\"1GQtx:UPEF[\\u7$Fjgv$\"1P'yi$>3$p#F[\\u7$Fjgv$\"1+++$f'*G v#F[\\uF``p-F$6$7S7$$\"1++++++&H\"FayF+7$Faaw$\"1\"3eJk.yD\"F_hv7$Faaw $\"1$RIw'=@_BF_hv7$Faaw$\"1xDjw$))He$F_hv7$Faaw$\"1KEErK#>#[F_hv7$Faaw $\"1xZgC(p\\0'F_hv7$Faaw$\"11dyF;;)>(F_hv7$Faaw$\"1w-FF1(=Q)F_hv7$Faaw $\"1T/&y;kgg*F_hv7$Faaw$\"1E?UqJj#3\"Fh^v7$Faaw$\"1\"GKbo_\"37Fh^v7$Fa aw$\"1\"4jFM5(=8Fh^v7$Faaw$\"1&3&Hz?Fh^v7$F aaw$\"1=(*3T(e\"Q?Fh^v7$Faaw$\"1#3$pi(3j;#Fh^v7$Faaw$\"1tv*fCV(zAFh^v7 $Faaw$\"1^5e_))>/CFh^v7$Faaw$\"1AdB@+rADFh^v7$Faaw$\"1u))3$*ROYEFh^v7$ Faaw$\"1LC%)Rp\"*fFFh^v7$Faaw$\"11kQL\"*R#)GFh^v7$Faaw$\"1:_[4Pi4IFh^v 7$Faaw$\"1?)yp2t.7$Fh^v7$Faaw$\"1!)*\\zU&)*RKFh^v7$Faaw$\"1a,e!pcNO$Fh ^v7$Faaw$\"1zY3jqW%[$Fh^v7$Faaw$\"1[*\\Xr99g$Fh^v7$Faaw$\"1%[j!\\tGJPF h^v7$Faaw$\"1oE\\(o$)z%QFh^v7$Faaw$\"1C05P>esRFh^v7$Faaw$\"1wj=qk[&3%F h^v7$Faaw$\"1tV%o%4#*3UFh^v7$Faaw$\"1Ve!*pJ1DVFh^v7$Faaw$\"12b'yDlkW%F h^v7$Faaw$\"1mu7*)z:lXFh^v7$Faaw$\"17*)))y#=%*o%Fh^v7$Faaw$\"1W\\h1d44 [Fh^v7$Faaw$\"1'pe=;$[J\\Fh^v7$Faaw$\"1Q4=(=dG0&Fh^v7$Faaw$\"1YM.6tQk^ Fh^v7$Faaw$\"1Tk8>M@#H&Fh^v7$Faaw$\"1d)\\ShUlS&Fh^v7$Faaw$\"1+#>G_IF_hv7$Fgjw$\"18m\"p+*p'G$F_hv7$Fgjw$\"1n KhH@d*\\$F_hv7$Fgjw$\"1J)G?'))p_PF_hv7$Fgjw$\"1oH22/8nRF_hv7$Fgjw$\"1w .dq_c;UF_hv7$Fgjw$\"1vm()*QdtV%F_hv7$Fgjw$\"1.TUQ8gzYF_hv7$Fgjw$\"19\\ 9j[F5\\F_hv7$Fgjw$\"1k#\\p[e4:&F_hv7$Fgjw$\"1;v*o=\")>P&F_hv7$Fgjw$\"1 B8#z-%Q5cF_hv7$Fgjw$\"1Bs%ydtA#y'F_hv7$Fgjw$\"1'Q4 &)\\'*)4qF_hv7$Fgjw$\"1FH*\\I%oisF_hv7$Fgjw$\"1$G5BRD)*[(F_hv7$Fgjw$\" 1'*eJ&)pMKxF_hv7$Fgjw$\"1@^P%f2@&zF_hv7$Fgjw$\"1&Hgk)RO#>)F_hv7$Fgjw$ \"1$)RzFlU=%)F_hv7$Fgjw$\"1m\"[/-FZl)F_hv7$Fgjw$\"1N`\"Q%Rv&)))F_hv7$F gjw$\"1\\)>@soGh0O*F_hv7$Fgjw$\"1`%3lsz()f*F_hv7$ Fgjw$\"1kLC(fD]$)*F_hv7$Fgjw$\"1w,p,6@05Fh^v7$Fgjw$\"1L2fS:4I5Fh^v7$Fg jw$\"1d_u,\\M_5Fh^v7$Fgjw$\"1g%ohirg2\"Fh^v7$Fgjw$\"1%G.*G6y)4\"Fh^v7$ Fgjw$\"1+++\"f&=B6Fh^vF``p-F$6$7S7$$\"1++++++&\\\"FayF+7$F^dx$\"1Q&*) \\l3wW%!#A7$F^dx$\"1&ohW;\\uJ)Fcdx7$F^dx$\"1PZ8b*[pE\"F\\[x7$F^dx$\"1 \"*Hp3y.0*y\\a&F\\[x7$F^dx$\"1)zHt_O3(fF\\[x7$F^dx$\"1g)QM.cvN'F \\[x7$F^dx$\"1$)R^IHSh!)F\\[x7$F^dx$\"1^c4idF,&)F\\[x7$F^dx$\"1k\\?l EL?*)F\\[x7$F^dx$\"1#eV;euvN*F\\[x7$F^dx$\"1I$)3`\"*4f(*F\\[x7$F^dx$\" 1&pZZw>#>5F_hv7$F^dx$\"1OKIEl?k5F_hv7$F^dx$\"1NS[MwO.6F_hv7$F^dx$\"1:H 8F_hv7$F^dx$\"1CH]7;lg8F_hv7$F^dx$\"13\\) =q4ZS\"F_hv7$F^dx$\"1;Y=zGjW9F_hv7$F^dx$\"1'=mqXz#)[\"F_hv7$F^dx$\"1!R &R&[Z$H:F_hv7$F^dx$\"1\"*Q@3aFs:F_hv7$F^dx$\"1\"*)y#*HXUh\"F_hv7$F^dx$ \"1\")ol&)Q=e;F_hv7$F^dx$\"1\\Pcm>]+F_hv7$F^dx$\"1+Le`M'[&>F_hv7$F^dx$\"1*30-0>h*>F_hv7 $F^dx$\"1+++BQXS?F_hvF``p-%+AXESLABELSG6$%!GFd]y-%*AXESSTYLEG6#%$BOXG- %%VIEWG6$%(DEFAULTGF\\^y" 1 2 0 1 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 39 "We can make the fo llowing observations:" }}{PARA 0 "" 0 "" {TEXT -1 200 "1) for the chos en moderate number of recordings (event/no-event per second) we can fi nd that the simulation does well on the event rate, and on the expecte d average number of events per time interval." }}{PARA 0 "" 0 "" {TEXT -1 156 "2) a reasonable result is obtained for the entire distri bution; note that the Poissonian drawn for comparison is based on inpu t from the data run (variable " }{TEXT 19 4 "mymu" }{TEXT -1 368 " for the average number of events per time interval), and not the predeter mined rate used to generate the sequence. Substantial discrepancies oc cur for individual bins in the histogram. This should not cause too mu ch of an alarm: the Poisson distribution represents a limiting probabi lity distribution that should be reached in the limit of many performe d experiments." }}{PARA 0 "" 0 "" {TEXT -1 237 "3) We should explore f urther the consistency of the data with a Poisson distribution: either by considering larger data sets, or by performing a systematic test o n the hypothesis (chi-squared test). The latter is attempted further b elow." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 311 12 "Exercise P1:" }}{PARA 0 "" 0 "" {TEXT -1 65 "Explore the above data b y choosing different subinterval lengths " }{TEXT 19 4 "tint" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 312 12 "Exercise P2:" }}{PARA 0 "" 0 "" {TEXT -1 94 "Repeat the above step s for cases where the decay rate is increased and decreased respective ly." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 423 "We now consider an interesting additional check, namely a test of the hypothesis that the distribution generated from the independen t events is indeed compatible with a Poisson distribution. A chi-squar ed test of a hypothesis is a quantitative measure that essentially com pares the data distributions as done in the above graph by summing the squares of the discrepancies, and then turning the result into a usef ul measure." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 51 "In the graph above we have used the data to deduce " }{XPPEDIT 18 0 "mu;" "6#%#muG" }{TEXT -1 91 " value for the Poissonian, and then generated a distribution of expected values to observe " }{XPPEDIT 18 0 "nu = 0 .. infinity;" "6#/%#nuG;\"\"!%)infinityG" }{TEXT -1 486 " events per time interval. (A limitation of the Poissonian is perhaps \+ noticable here: our simulation, and any Geiger counter sequence has a \+ maximum number of observable events, namely when the counter is clicki ng all the time, while the Poissonian distributes the probability over an infinite interval; this is the price to be paid for the replacemen t of a binomial by a Poisson distribution; in practice this shouldn't \+ be a problem, as the predictions for high event counts are small)." }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 178 "We can c alculate the chi-squared for the results obtained above (make sure tha t the statements that produce the graph above were executed) from the \+ expected and observed results:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "chi_sq:=add((Nt*PD(mymu,nu)-C[nu])^2/(Nt*PD(mymu,nu)),nu=0..Em ax);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'chi_sqG$\"+>(4H]$!\"*" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 131 "Is this an acceptable value or no t? This depends on the number of degrees of freedom. Here we have a li ttle problem: should we use " }{TEXT 19 4 "Emax" }{TEXT -1 195 ", whic h includes data where we have no counts and small expected probabiliti es, or not? After all, there was some arbitrariness in choosing this v ariable, which acts as a cut-off for the diagram." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 103 "The criterion for accept ability is that chi-squared be much less than the number of degrees of freedom " }{TEXT 314 1 "d" }{TEXT -1 76 "; rejection of the hypothesi s is given when chi-squared is much bigger than " }{TEXT 313 1 "d" } {TEXT -1 42 ". The question is whether we can say that " }{TEXT 315 1 "d" }{TEXT -1 5 "=16 (" }{TEXT 19 7 "=Emax+1" }{TEXT -1 2 ")?" }} {PARA 0 "" 0 "" {TEXT -1 82 "Judging from the graph the right attitude should be to pick the 10 visible points:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 63 "chi_sq:=add((Nt*PD(mymu,nu)-C[nu])^2/(Nt*PD(mymu,nu )),nu=0..9);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'chi_sqG$\"+PzQQM!\" *" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 515 "The conclusion based on the limited data set would be that the data are not inconsistent with a P oisson distribution, and that further testing is desirable. Further te sting means to run with bigger data sets. An important issue when perf orming these tests is to choose reasonable time intervals (which depen d on the actual count rate). The chi-squared comparison has to have a \+ reasonable number of contributing data points (say six or more to acco unt for two degrees of freedom to be subtracted as explained below)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 86 "One sho uld also be a bit more precise when assessing the number of degrees of freedom " }{TEXT 316 1 "d" }{TEXT -1 275 ". We have simply taken the \+ number of non-negligible theoretical data points (we could have used t he number of data bins). It should be kept in mind that we really have one less independent value, as the data were used to infer the averag e number of decays per time interval (" }{TEXT 19 4 "mymu" }{TEXT -1 245 "). There can be other constraints (usually the number of data poi nts; in the case of a Gaussian distribution also its deviation) which \+ are used to determine a reduced number of degrees of freedom. A reduce d chi-squared is defined as chi-squared/" }{TEXT 317 1 "d" }{TEXT -1 217 ", and for this quantity the borderline between acceptance and non -acceptance of the hypothesis is the value of one. It is evident that \+ for acceptance we need to obtain much less than one, and vice versa fo r rejection." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 124 "A healthier attitude for dealing with the open-ended Poisson d istribution would be to call the bin neighbouring the largest " } {XPPEDIT 18 0 "nu;" "6#%#nuG" }{TEXT -1 100 " value that still receive d non-zero counts to be one bin for all events beyond the largest obse rved " }{XPPEDIT 18 0 "nu;" "6#%#nuG" }{TEXT -1 34 " value. In our cas e that would be " }{TEXT 19 4 "nu=9" }{TEXT -1 23 " (corresponding now to " }{XPPEDIT 18 0 "9 <= nu;" "6#1\"\"*%#nuG" }{TEXT -1 30 ", and re ceiving 0 counts) and " }{TEXT 318 1 "d" }{TEXT -1 132 " would be 10-2 =8, i.e., a borderline result [note that the precise value in the latt er statement can change for repeated data runs]." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "chi_sq/8;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$ \"+@\\)zH%!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 75 "Our own test wit h the code above produced the following interesting result:" }}{PARA 0 "" 0 "" {TEXT -1 171 "- for data runs with up to 3000 time ticks (at a decay rate of 1/10 per tick) consistency with a Poissonian was obse rved (reduced chi-squared observed between 0.4 and 1.4)" }}{PARA 0 "" 0 "" {TEXT -1 184 "- for data runs with 5000 time ticks and more the P oisson hypothesis was not confirmed (reduced chi-squared often around \+ 1.5 rather than falling further below one in a consistent way)." }} {PARA 0 "" 0 "" {TEXT -1 236 "This allows one to conclude that at the \+ level of random number sequences in the hundreds of thousands the pseu do-random number generator has problems with correlations, i.e., the r andom numbers are not truly independent from each other." }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 241 "We can also concl ude that it is hard to test a hypothesis, i.e., we have to run long si mulations, and we have to know that the tools are accurate enough. Sho rt simulations are often sufficient in order to obtain a qualitative u nderstanding. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 319 12 "Exercise P3:" }}{PARA 0 "" 0 "" {TEXT -1 95 "Explore the chi-squared test of the Poisson distribution hypothesis with differen t count rates." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "4 11 1 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }