HP 40gs HP 39gs_40gs_Mastering The Graphing Calculator_English_E_F2224-90010.p - Page 121

Simulation of a normal die, Simulating Random Variables

Page 121 highlights

Simulation of a normal die Similarly the expression INT(RANDOM*6+1) will simulate one roll of the die. This means that MAKELIST(INT(RANDOM*6+1),X,1,500,1) will simulate 500 rolls of a normal die. We therefore need only store the resulting list into a Statistics aplet column to analyze and graph it. This is shown in the series of screen shots to the right. To simulate the adding of two six sided dice, use the similar expression MAKELIST(INT(RANDOM*6+1)+INT(RANDOM*6+1),X,1,500,1) Simulating Random Variables It is occasionally handy to be able to simulate a set of observations on random variables of various sorts. This can also be done using the MATH menu function MAKELIST. See the previous page for the syntax of MAKELIST. It is also covered in more detail on page 190) Example 1: Simulate 100 observations on a Uniform[0,1] random variable to be stored in C2. In the HOME view type: MAKELIST(RANDOM,X,1,100,1) C2 Note: 1. The symbol is one of the screen keys in the HOME view, appearing as and read as 'store'. 2. The X is simply a dummy variable here to count off the values. Example 2: Simulate 50 observations on a discrete uniform[3,7] random variable In HOME type: MAKELIST(INT(5*RANDOM+3),X,1,50,1) C2 Example 3: Simulate 50 observations on a Binomial random variable with n = 20 and p = 0.75. In HOME type: MAKELIST(Σ(I=1,20,RANDOM

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
  • 115
  • 116
  • 117
  • 118
  • 119
  • 120
  • 121
  • 122
  • 123
  • 124
  • 125
  • 126
  • 127
  • 128
  • 129
  • 130
  • 131
  • 132
  • 133
  • 134
  • 135
  • 136
  • 137
  • 138
  • 139
  • 140
  • 141
  • 142
  • 143
  • 144
  • 145
  • 146
  • 147
  • 148
  • 149
  • 150
  • 151
  • 152
  • 153
  • 154
  • 155
  • 156
  • 157
  • 158
  • 159
  • 160
  • 161
  • 162
  • 163
  • 164
  • 165
  • 166
  • 167
  • 168
  • 169
  • 170
  • 171
  • 172
  • 173
  • 174
  • 175
  • 176
  • 177
  • 178
  • 179
  • 180
  • 181
  • 182
  • 183
  • 184
  • 185
  • 186
  • 187
  • 188
  • 189
  • 190
  • 191
  • 192
  • 193
  • 194
  • 195
  • 196
  • 197
  • 198
  • 199
  • 200
  • 201
  • 202
  • 203
  • 204
  • 205
  • 206
  • 207
  • 208
  • 209
  • 210
  • 211
  • 212
  • 213
  • 214
  • 215
  • 216
  • 217
  • 218
  • 219
  • 220
  • 221
  • 222
  • 223
  • 224
  • 225
  • 226
  • 227
  • 228
  • 229
  • 230
  • 231
  • 232
  • 233
  • 234
  • 235
  • 236
  • 237
  • 238
  • 239
  • 240
  • 241
  • 242
  • 243
  • 244
  • 245
  • 246
  • 247
  • 248
  • 249
  • 250
  • 251
  • 252
  • 253
  • 254
  • 255
  • 256
  • 257
  • 258
  • 259
  • 260
  • 261
  • 262
  • 263
  • 264
  • 265
  • 266
  • 267
  • 268
  • 269
  • 270
  • 271
  • 272
  • 273
  • 274
  • 275
  • 276
  • 277
  • 278
  • 279
  • 280
  • 281
  • 282
  • 283
  • 284
  • 285
  • 286
  • 287
  • 288
  • 289
  • 290
  • 291
  • 292
  • 293
  • 294
  • 295
  • 296
  • 297
  • 298
  • 299
  • 300
  • 301
  • 302
  • 303
  • 304
  • 305
  • 306
  • 307
  • 308
  • 309
  • 310
  • 311
  • 312
  • 313
  • 314
  • 315
  • 316
  • 317
  • 318
  • 319
  • 320
  • 321
  • 322
  • 323
  • 324
  • 325
  • 326
  • 327
  • 328
  • 329
  • 330
  • 331
  • 332
  • 333
  • 334
  • 335
  • 336
  • 337
  • 338
  • 339
  • 340
  • 341
  • 342
  • 343
  • 344
  • 345
  • 346
  • 347
  • 348
  • 349
  • 350
  • 351
  • 352
  • 353
  • 354
  • 355
  • 356
  • 357
  • 358
  • 359
  • 360
  • 361
  • 362
  • 363
  • 364
  • 365
  • 366

Simulation of a normal die
Similarly the expression
INT(RANDOM*6+1)
will simulate one roll of the
die. This means that
MAKELIST(INT(RANDOM*6+1),X,1,500,1)
will
simulate 500 rolls of a normal die. We therefore need only store the
resulting list into a Statistics aplet column to analyze and graph it.
This
is shown in the series of screen shots to the right.
To simulate the adding of two six sided dice, use the similar expression
MAKELIST(INT(RANDOM*6+1)+INT(RANDOM*6+1),X,1,500,1)
Simulating Random Variables
It is occasionally handy to be able to simulate a set of observations on random variables of various sorts.
This
can also be done using the
MATH
menu function
MAKELIST
.
See the previous page for the syntax of
MAKELIST
. It is also covered in more detail on page 190)
Example 1:
Simulate 100 observations on a Uniform[0,1] random variable to be stored in
C2
.
In the
HOME
view type:
MAKELIST(RANDOM,X,1,100,1)
C2
symbol is one of the screen keys in the
HOME
view,
appearing as
Note: 1. The
and read as ‘store’.
2. The
X
is simply a dummy variable here to count off the values.
Example 2:
Simulate 50 observations on a discrete uniform[3,7] random variable
In
HOME
type:
MAKELIST(INT(5*RANDOM+3),X,1,50,1)
C2
Example 3:
Simulate 50 observations on a Binomial random variable with n = 20 and p = 0.75.
In
HOME
type:
MAKELIST(
Σ
(I=1,20,RANDOM<0.75),X,1,50,1)
C2
Note: The
Σ
and the = are both on the keyboard.
This will be a relatively slow calculation because it
involves evaluating 1000 random numbers.
121