library(deSolve)
## Warning: package 'deSolve' was built under R version 3.4.4
metopob <- function (t, y, parms) {
f = y[1]
with(as.list(parms), {
dp = (pi*f*(1-f))-pe*f
return(list(dp))
})
}
prms <- c(pi=0.15, pe=0.15)
Initial.f=c(0.9)
out.M <- data.frame(ode(y = Initial.f, parms = prms, times = 1:1000, func = metopob))
out.M
## time X1
## 1 1 0.900000000
## 2 2 0.792952860
## 3 3 0.708661511
## 4 4 0.640569162
## 5 5 0.584415967
## 6 6 0.537313893
## 7 7 0.497238048
## 8 8 0.462725686
## 9 9 0.432693179
## 10 10 0.406321419
## 11 11 0.382979574
## 12 12 0.362173942
## 13 13 0.343512430
## 14 14 0.326679750
## 15 15 0.311419653
## 16 16 0.297521596
## 17 17 0.284811020
## 18 18 0.273141972
## 19 19 0.262391475
## 20 20 0.252455177
## 21 21 0.243243960
## 22 22 0.234681250
## 23 23 0.226700890
## 24 24 0.219245425
## 25 25 0.212264719
## 26 26 0.205714823
## 27 27 0.199557049
## 28 28 0.193757209
## 29 29 0.188284975
## 30 30 0.183113353
## 31 31 0.178218233
## 32 32 0.173578019
## 33 33 0.169173305
## 34 34 0.164986606
## 35 35 0.161002128
## 36 36 0.157205563
## 37 37 0.153583927
## 38 38 0.150125400
## 39 39 0.146819206
## 40 40 0.143655499
## 41 41 0.140625260
## 42 42 0.137720219
## 43 43 0.134932774
## 44 44 0.132255925
## 45 45 0.129683219
## 46 46 0.127208694
## 47 47 0.124826835
## 48 48 0.122532533
## 49 49 0.120321047
## 50 50 0.118187972
## 51 51 0.116129210
## 52 52 0.114140946
## 53 53 0.112219618
## 54 54 0.110361902
## 55 55 0.108564691
## 56 56 0.106825077
## 57 57 0.105140333
## 58 58 0.103507905
## 59 59 0.101925392
## 60 60 0.100390541
## 61 61 0.098901228
## 62 62 0.097455459
## 63 63 0.096051350
## 64 64 0.094687126
## 65 65 0.093361111
## 66 66 0.092071724
## 67 67 0.090817465
## 68 68 0.089596921
## 69 69 0.088408748
## 70 70 0.087251676
## 71 71 0.086124500
## 72 72 0.085026076
## 73 73 0.083955317
## 74 74 0.082911192
## 75 75 0.081892719
## 76 76 0.080898963
## 77 77 0.079929037
## 78 78 0.078982092
## 79 79 0.078057323
## 80 80 0.077153958
## 81 81 0.076271264
## 82 82 0.075408538
## 83 83 0.074565111
## 84 84 0.073740342
## 85 85 0.072933620
## 86 86 0.072144358
## 87 87 0.071371995
## 88 88 0.070615994
## 89 89 0.069875841
## 90 90 0.069151043
## 91 91 0.068441127
## 92 92 0.067745639
## 93 93 0.067064143
## 94 94 0.066396222
## 95 95 0.065741474
## 96 96 0.065099514
## 97 97 0.064469969
## 98 98 0.063852484
## 99 99 0.063246715
## 100 100 0.062652332
## 101 101 0.062069017
## 102 102 0.061496463
## 103 103 0.060934376
## 104 104 0.060382470
## 105 105 0.059840473
## 106 106 0.059308119
## 107 107 0.058785154
## 108 108 0.058271330
## 109 109 0.057766411
## 110 110 0.057270168
## 111 111 0.056782377
## 112 112 0.056302826
## 113 113 0.055831307
## 114 114 0.055367620
## 115 115 0.054911571
## 116 116 0.054462974
## 117 117 0.054021647
## 118 118 0.053587415
## 119 119 0.053160108
## 120 120 0.052739562
## 121 121 0.052325618
## 122 122 0.051918121
## 123 123 0.051516921
## 124 124 0.051121875
## 125 125 0.050732841
## 126 126 0.050349684
## 127 127 0.049972271
## 128 128 0.049600474
## 129 129 0.049234168
## 130 130 0.048873233
## 131 131 0.048517551
## 132 132 0.048167010
## 133 133 0.047821497
## 134 134 0.047480906
## 135 135 0.047145132
## 136 136 0.046814073
## 137 137 0.046487632
## 138 138 0.046165712
## 139 139 0.045848219
## 140 140 0.045535064
## 141 141 0.045226158
## 142 142 0.044921414
## 143 143 0.044620750
## 144 144 0.044324084
## 145 145 0.044031337
## 146 146 0.043742431
## 147 147 0.043457292
## 148 148 0.043175846
## 149 149 0.042898023
## 150 150 0.042623751
## 151 151 0.042352965
## 152 152 0.042085597
## 153 153 0.041821585
## 154 154 0.041560863
## 155 155 0.041303373
## 156 156 0.041049053
## 157 157 0.040797846
## 158 158 0.040549695
## 159 159 0.040304545
## 160 160 0.040062340
## 161 161 0.039823030
## 162 162 0.039586561
## 163 163 0.039352885
## 164 164 0.039121950
## 165 165 0.038893711
## 166 166 0.038668119
## 167 167 0.038445129
## 168 168 0.038224696
## 169 169 0.038006776
## 170 170 0.037791327
## 171 171 0.037578307
## 172 172 0.037367675
## 173 173 0.037159391
## 174 174 0.036953416
## 175 175 0.036749712
## 176 176 0.036548241
## 177 177 0.036348967
## 178 178 0.036151855
## 179 179 0.035956869
## 180 180 0.035763975
## 181 181 0.035573139
## 182 182 0.035384330
## 183 183 0.035197514
## 184 184 0.035012660
## 185 185 0.034829737
## 186 186 0.034648717
## 187 187 0.034469568
## 188 188 0.034292262
## 189 189 0.034116771
## 190 190 0.033943066
## 191 191 0.033771122
## 192 192 0.033600911
## 193 193 0.033432407
## 194 194 0.033265585
## 195 195 0.033100419
## 196 196 0.032936885
## 197 197 0.032774960
## 198 198 0.032614618
## 199 199 0.032455838
## 200 200 0.032298596
## 201 201 0.032142871
## 202 202 0.031988640
## 203 203 0.031835882
## 204 204 0.031684576
## 205 205 0.031534701
## 206 206 0.031386238
## 207 207 0.031239166
## 208 208 0.031093466
## 209 209 0.030949119
## 210 210 0.030806105
## 211 211 0.030664408
## 212 212 0.030524008
## 213 213 0.030384887
## 214 214 0.030247029
## 215 215 0.030110417
## 216 216 0.029975033
## 217 217 0.029840861
## 218 218 0.029707884
## 219 219 0.029576088
## 220 220 0.029445456
## 221 221 0.029315972
## 222 222 0.029187623
## 223 223 0.029060392
## 224 224 0.028934266
## 225 225 0.028809230
## 226 226 0.028685270
## 227 227 0.028562372
## 228 228 0.028440523
## 229 229 0.028319709
## 230 230 0.028199917
## 231 231 0.028081134
## 232 232 0.027963347
## 233 233 0.027846545
## 234 234 0.027730714
## 235 235 0.027615843
## 236 236 0.027501920
## 237 237 0.027388933
## 238 238 0.027276870
## 239 239 0.027165721
## 240 240 0.027055473
## 241 241 0.026946117
## 242 242 0.026837642
## 243 243 0.026730036
## 244 244 0.026623290
## 245 245 0.026517393
## 246 246 0.026412335
## 247 247 0.026308106
## 248 248 0.026204697
## 249 249 0.026102097
## 250 250 0.026000298
## 251 251 0.025899289
## 252 252 0.025799063
## 253 253 0.025699609
## 254 254 0.025600919
## 255 255 0.025502984
## 256 256 0.025405795
## 257 257 0.025309345
## 258 258 0.025213624
## 259 259 0.025118624
## 260 260 0.025024338
## 261 261 0.024930756
## 262 262 0.024837872
## 263 263 0.024745678
## 264 264 0.024654165
## 265 265 0.024563327
## 266 266 0.024473155
## 267 267 0.024383644
## 268 268 0.024294784
## 269 269 0.024206570
## 270 270 0.024118995
## 271 271 0.024032050
## 272 272 0.023945731
## 273 273 0.023860029
## 274 274 0.023774938
## 275 275 0.023690452
## 276 276 0.023606565
## 277 277 0.023523269
## 278 278 0.023440559
## 279 279 0.023358429
## 280 280 0.023276873
## 281 281 0.023195883
## 282 282 0.023115456
## 283 283 0.023035584
## 284 284 0.022956263
## 285 285 0.022877485
## 286 286 0.022799247
## 287 287 0.022721542
## 288 288 0.022644365
## 289 289 0.022567710
## 290 290 0.022491572
## 291 291 0.022415947
## 292 292 0.022340828
## 293 293 0.022266211
## 294 294 0.022192091
## 295 295 0.022118463
## 296 296 0.022045322
## 297 297 0.021972663
## 298 298 0.021900481
## 299 299 0.021828772
## 300 300 0.021757531
## 301 301 0.021686753
## 302 302 0.021616435
## 303 303 0.021546571
## 304 304 0.021477157
## 305 305 0.021408189
## 306 306 0.021339662
## 307 307 0.021271573
## 308 308 0.021203917
## 309 309 0.021136690
## 310 310 0.021069888
## 311 311 0.021003506
## 312 312 0.020937542
## 313 313 0.020871991
## 314 314 0.020806849
## 315 315 0.020742112
## 316 316 0.020677777
## 317 317 0.020613840
## 318 318 0.020550297
## 319 319 0.020487144
## 320 320 0.020424379
## 321 321 0.020361996
## 322 322 0.020299994
## 323 323 0.020238368
## 324 324 0.020177116
## 325 325 0.020116232
## 326 326 0.020055716
## 327 327 0.019995562
## 328 328 0.019935768
## 329 329 0.019876330
## 330 330 0.019817246
## 331 331 0.019758512
## 332 332 0.019700126
## 333 333 0.019642083
## 334 334 0.019584381
## 335 335 0.019527017
## 336 336 0.019469989
## 337 337 0.019413292
## 338 338 0.019356925
## 339 339 0.019300884
## 340 340 0.019245167
## 341 341 0.019189770
## 342 342 0.019134692
## 343 343 0.019079929
## 344 344 0.019025478
## 345 345 0.018971337
## 346 346 0.018917503
## 347 347 0.018863975
## 348 348 0.018810748
## 349 349 0.018757820
## 350 350 0.018705190
## 351 351 0.018652854
## 352 352 0.018600811
## 353 353 0.018549057
## 354 354 0.018497590
## 355 355 0.018446407
## 356 356 0.018395508
## 357 357 0.018344888
## 358 358 0.018294547
## 359 359 0.018244480
## 360 360 0.018194688
## 361 361 0.018145166
## 362 362 0.018095913
## 363 363 0.018046926
## 364 364 0.017998204
## 365 365 0.017949745
## 366 366 0.017901546
## 367 367 0.017853605
## 368 368 0.017805920
## 369 369 0.017758489
## 370 370 0.017711310
## 371 371 0.017664381
## 372 372 0.017617700
## 373 373 0.017571265
## 374 374 0.017525075
## 375 375 0.017479126
## 376 376 0.017433418
## 377 377 0.017387948
## 378 378 0.017342715
## 379 379 0.017297717
## 380 380 0.017252951
## 381 381 0.017208417
## 382 382 0.017164112
## 383 383 0.017120034
## 384 384 0.017076183
## 385 385 0.017032555
## 386 386 0.016989150
## 387 387 0.016945965
## 388 388 0.016902999
## 389 389 0.016860251
## 390 390 0.016817718
## 391 391 0.016775400
## 392 392 0.016733294
## 393 393 0.016691398
## 394 394 0.016649712
## 395 395 0.016608234
## 396 396 0.016566962
## 397 397 0.016525894
## 398 398 0.016485030
## 399 399 0.016444367
## 400 400 0.016403904
## 401 401 0.016363640
## 402 402 0.016323573
## 403 403 0.016283702
## 404 404 0.016244025
## 405 405 0.016204541
## 406 406 0.016165248
## 407 407 0.016126146
## 408 408 0.016087232
## 409 409 0.016048506
## 410 410 0.016009965
## 411 411 0.015971609
## 412 412 0.015933437
## 413 413 0.015895447
## 414 414 0.015857637
## 415 415 0.015820007
## 416 416 0.015782555
## 417 417 0.015745280
## 418 418 0.015708180
## 419 419 0.015671255
## 420 420 0.015634503
## 421 421 0.015597924
## 422 422 0.015561514
## 423 423 0.015525275
## 424 424 0.015489204
## 425 425 0.015453300
## 426 426 0.015417562
## 427 427 0.015381989
## 428 428 0.015346580
## 429 429 0.015311333
## 430 430 0.015276249
## 431 431 0.015241324
## 432 432 0.015206559
## 433 433 0.015171952
## 434 434 0.015137502
## 435 435 0.015103208
## 436 436 0.015069070
## 437 437 0.015035085
## 438 438 0.015001253
## 439 439 0.014967573
## 440 440 0.014934044
## 441 441 0.014900665
## 442 442 0.014867435
## 443 443 0.014834353
## 444 444 0.014801417
## 445 445 0.014768628
## 446 446 0.014735983
## 447 447 0.014703483
## 448 448 0.014671125
## 449 449 0.014638910
## 450 450 0.014606836
## 451 451 0.014574902
## 452 452 0.014543107
## 453 453 0.014511451
## 454 454 0.014479932
## 455 455 0.014448550
## 456 456 0.014417304
## 457 457 0.014386192
## 458 458 0.014355214
## 459 459 0.014324370
## 460 460 0.014293658
## 461 461 0.014263077
## 462 462 0.014232627
## 463 463 0.014202307
## 464 464 0.014172115
## 465 465 0.014142052
## 466 466 0.014112116
## 467 467 0.014082306
## 468 468 0.014052622
## 469 469 0.014023063
## 470 470 0.013993628
## 471 471 0.013964316
## 472 472 0.013935127
## 473 473 0.013906059
## 474 474 0.013877113
## 475 475 0.013848287
## 476 476 0.013819580
## 477 477 0.013790992
## 478 478 0.013762523
## 479 479 0.013734170
## 480 480 0.013705934
## 481 481 0.013677814
## 482 482 0.013649809
## 483 483 0.013621919
## 484 484 0.013594142
## 485 485 0.013566478
## 486 486 0.013538927
## 487 487 0.013511487
## 488 488 0.013484159
## 489 489 0.013456940
## 490 490 0.013429832
## 491 491 0.013402832
## 492 492 0.013375941
## 493 493 0.013349157
## 494 494 0.013322480
## 495 495 0.013295910
## 496 496 0.013269446
## 497 497 0.013243086
## 498 498 0.013216832
## 499 499 0.013190681
## 500 500 0.013164633
## 501 501 0.013138688
## 502 502 0.013112846
## 503 503 0.013087104
## 504 504 0.013061464
## 505 505 0.013035923
## 506 506 0.013010483
## 507 507 0.012985141
## 508 508 0.012959899
## 509 509 0.012934754
## 510 510 0.012909706
## 511 511 0.012884755
## 512 512 0.012859901
## 513 513 0.012835142
## 514 514 0.012810478
## 515 515 0.012785909
## 516 516 0.012761434
## 517 517 0.012737053
## 518 518 0.012712764
## 519 519 0.012688568
## 520 520 0.012664464
## 521 521 0.012640452
## 522 522 0.012616530
## 523 523 0.012592698
## 524 524 0.012568957
## 525 525 0.012545305
## 526 526 0.012521741
## 527 527 0.012498266
## 528 528 0.012474879
## 529 529 0.012451579
## 530 530 0.012428366
## 531 531 0.012405240
## 532 532 0.012382199
## 533 533 0.012359244
## 534 534 0.012336374
## 535 535 0.012313588
## 536 536 0.012290886
## 537 537 0.012268268
## 538 538 0.012245733
## 539 539 0.012223281
## 540 540 0.012200910
## 541 541 0.012178622
## 542 542 0.012156414
## 543 543 0.012134288
## 544 544 0.012112242
## 545 545 0.012090276
## 546 546 0.012068389
## 547 547 0.012046582
## 548 548 0.012024853
## 549 549 0.012003203
## 550 550 0.011981630
## 551 551 0.011960135
## 552 552 0.011938716
## 553 553 0.011917375
## 554 554 0.011896109
## 555 555 0.011874919
## 556 556 0.011853805
## 557 557 0.011832765
## 558 558 0.011811801
## 559 559 0.011790910
## 560 560 0.011770093
## 561 561 0.011749349
## 562 562 0.011728678
## 563 563 0.011708080
## 564 564 0.011687555
## 565 565 0.011667101
## 566 566 0.011646718
## 567 567 0.011626407
## 568 568 0.011606166
## 569 569 0.011585996
## 570 570 0.011565895
## 571 571 0.011545865
## 572 572 0.011525903
## 573 573 0.011506010
## 574 574 0.011486186
## 575 575 0.011466431
## 576 576 0.011446743
## 577 577 0.011427122
## 578 578 0.011407569
## 579 579 0.011388082
## 580 580 0.011368662
## 581 581 0.011349308
## 582 582 0.011330020
## 583 583 0.011310797
## 584 584 0.011291640
## 585 585 0.011272547
## 586 586 0.011253518
## 587 587 0.011234554
## 588 588 0.011215654
## 589 589 0.011196817
## 590 590 0.011178043
## 591 591 0.011159332
## 592 592 0.011140684
## 593 593 0.011122098
## 594 594 0.011103573
## 595 595 0.011085111
## 596 596 0.011066709
## 597 597 0.011048369
## 598 598 0.011030089
## 599 599 0.011011870
## 600 600 0.010993711
## 601 601 0.010975611
## 602 602 0.010957571
## 603 603 0.010939591
## 604 604 0.010921669
## 605 605 0.010903806
## 606 606 0.010886001
## 607 607 0.010868254
## 608 608 0.010850565
## 609 609 0.010832934
## 610 610 0.010815359
## 611 611 0.010797842
## 612 612 0.010780381
## 613 613 0.010762977
## 614 614 0.010745629
## 615 615 0.010728336
## 616 616 0.010711099
## 617 617 0.010693918
## 618 618 0.010676791
## 619 619 0.010659720
## 620 620 0.010642702
## 621 621 0.010625739
## 622 622 0.010608830
## 623 623 0.010591975
## 624 624 0.010575173
## 625 625 0.010558425
## 626 626 0.010541729
## 627 627 0.010525086
## 628 628 0.010508496
## 629 629 0.010491958
## 630 630 0.010475471
## 631 631 0.010459037
## 632 632 0.010442654
## 633 633 0.010426322
## 634 634 0.010410041
## 635 635 0.010393811
## 636 636 0.010377632
## 637 637 0.010361503
## 638 638 0.010345424
## 639 639 0.010329394
## 640 640 0.010313415
## 641 641 0.010297484
## 642 642 0.010281603
## 643 643 0.010265771
## 644 644 0.010249987
## 645 645 0.010234252
## 646 646 0.010218565
## 647 647 0.010202926
## 648 648 0.010187335
## 649 649 0.010171792
## 650 650 0.010156295
## 651 651 0.010140846
## 652 652 0.010125444
## 653 653 0.010110089
## 654 654 0.010094780
## 655 655 0.010079518
## 656 656 0.010064301
## 657 657 0.010049130
## 658 658 0.010034005
## 659 659 0.010018926
## 660 660 0.010003892
## 661 661 0.009988903
## 662 662 0.009973958
## 663 663 0.009959059
## 664 664 0.009944203
## 665 665 0.009929392
## 666 666 0.009914625
## 667 667 0.009899902
## 668 668 0.009885223
## 669 669 0.009870587
## 670 670 0.009855994
## 671 671 0.009841445
## 672 672 0.009826938
## 673 673 0.009812474
## 674 674 0.009798053
## 675 675 0.009783673
## 676 676 0.009769336
## 677 677 0.009755041
## 678 678 0.009740788
## 679 679 0.009726576
## 680 680 0.009712406
## 681 681 0.009698277
## 682 682 0.009684189
## 683 683 0.009670142
## 684 684 0.009656136
## 685 685 0.009642170
## 686 686 0.009628244
## 687 687 0.009614359
## 688 688 0.009600513
## 689 689 0.009586708
## 690 690 0.009572942
## 691 691 0.009559215
## 692 692 0.009545528
## 693 693 0.009531880
## 694 694 0.009518271
## 695 695 0.009504701
## 696 696 0.009491169
## 697 697 0.009477676
## 698 698 0.009464221
## 699 699 0.009450805
## 700 700 0.009437426
## 701 701 0.009424085
## 702 702 0.009410782
## 703 703 0.009397516
## 704 704 0.009384288
## 705 705 0.009371097
## 706 706 0.009357942
## 707 707 0.009344825
## 708 708 0.009331745
## 709 709 0.009318701
## 710 710 0.009305693
## 711 711 0.009292722
## 712 712 0.009279787
## 713 713 0.009266887
## 714 714 0.009254024
## 715 715 0.009241196
## 716 716 0.009228404
## 717 717 0.009215647
## 718 718 0.009202926
## 719 719 0.009190239
## 720 720 0.009177587
## 721 721 0.009164971
## 722 722 0.009152388
## 723 723 0.009139841
## 724 724 0.009127327
## 725 725 0.009114848
## 726 726 0.009102403
## 727 727 0.009089992
## 728 728 0.009077615
## 729 729 0.009065271
## 730 730 0.009052961
## 731 731 0.009040684
## 732 732 0.009028441
## 733 733 0.009016230
## 734 734 0.009004053
## 735 735 0.008991908
## 736 736 0.008979797
## 737 737 0.008967717
## 738 738 0.008955671
## 739 739 0.008943656
## 740 740 0.008931674
## 741 741 0.008919724
## 742 742 0.008907805
## 743 743 0.008895919
## 744 744 0.008884064
## 745 745 0.008872241
## 746 746 0.008860449
## 747 747 0.008848688
## 748 748 0.008836959
## 749 749 0.008825261
## 750 750 0.008813594
## 751 751 0.008801957
## 752 752 0.008790351
## 753 753 0.008778776
## 754 754 0.008767231
## 755 755 0.008755717
## 756 756 0.008744232
## 757 757 0.008732778
## 758 758 0.008721354
## 759 759 0.008709959
## 760 760 0.008698595
## 761 761 0.008687260
## 762 762 0.008675954
## 763 763 0.008664678
## 764 764 0.008653431
## 765 765 0.008642213
## 766 766 0.008631025
## 767 767 0.008619865
## 768 768 0.008608734
## 769 769 0.008597632
## 770 770 0.008586558
## 771 771 0.008575513
## 772 772 0.008564496
## 773 773 0.008553508
## 774 774 0.008542548
## 775 775 0.008531615
## 776 776 0.008520711
## 777 777 0.008509835
## 778 778 0.008498986
## 779 779 0.008488165
## 780 780 0.008477371
## 781 781 0.008466605
## 782 782 0.008455866
## 783 783 0.008445154
## 784 784 0.008434470
## 785 785 0.008423812
## 786 786 0.008413182
## 787 787 0.008402578
## 788 788 0.008392001
## 789 789 0.008381450
## 790 790 0.008370926
## 791 791 0.008360428
## 792 792 0.008349957
## 793 793 0.008339512
## 794 794 0.008329093
## 795 795 0.008318700
## 796 796 0.008308332
## 797 797 0.008297991
## 798 798 0.008287675
## 799 799 0.008277385
## 800 800 0.008267121
## 801 801 0.008256882
## 802 802 0.008246668
## 803 803 0.008236479
## 804 804 0.008226316
## 805 805 0.008216178
## 806 806 0.008206064
## 807 807 0.008195976
## 808 808 0.008185912
## 809 809 0.008175873
## 810 810 0.008165859
## 811 811 0.008155869
## 812 812 0.008145903
## 813 813 0.008135962
## 814 814 0.008126045
## 815 815 0.008116152
## 816 816 0.008106283
## 817 817 0.008096438
## 818 818 0.008086618
## 819 819 0.008076820
## 820 820 0.008067047
## 821 821 0.008057297
## 822 822 0.008047571
## 823 823 0.008037868
## 824 824 0.008028189
## 825 825 0.008018533
## 826 826 0.008008900
## 827 827 0.007999290
## 828 828 0.007989703
## 829 829 0.007980139
## 830 830 0.007970598
## 831 831 0.007961080
## 832 832 0.007951585
## 833 833 0.007942112
## 834 834 0.007932661
## 835 835 0.007923234
## 836 836 0.007913828
## 837 837 0.007904445
## 838 838 0.007895084
## 839 839 0.007885745
## 840 840 0.007876428
## 841 841 0.007867134
## 842 842 0.007857861
## 843 843 0.007848610
## 844 844 0.007839381
## 845 845 0.007830173
## 846 846 0.007820987
## 847 847 0.007811823
## 848 848 0.007802680
## 849 849 0.007793558
## 850 850 0.007784458
## 851 851 0.007775379
## 852 852 0.007766321
## 853 853 0.007757284
## 854 854 0.007748268
## 855 855 0.007739273
## 856 856 0.007730299
## 857 857 0.007721346
## 858 858 0.007712413
## 859 859 0.007703502
## 860 860 0.007694610
## 861 861 0.007685739
## 862 862 0.007676889
## 863 863 0.007668059
## 864 864 0.007659249
## 865 865 0.007650460
## 866 866 0.007641690
## 867 867 0.007632941
## 868 868 0.007624212
## 869 869 0.007615503
## 870 870 0.007606813
## 871 871 0.007598143
## 872 872 0.007589494
## 873 873 0.007580863
## 874 874 0.007572253
## 875 875 0.007563662
## 876 876 0.007555090
## 877 877 0.007546538
## 878 878 0.007538005
## 879 879 0.007529491
## 880 880 0.007520997
## 881 881 0.007512522
## 882 882 0.007504065
## 883 883 0.007495628
## 884 884 0.007487210
## 885 885 0.007478811
## 886 886 0.007470430
## 887 887 0.007462069
## 888 888 0.007453726
## 889 889 0.007445401
## 890 890 0.007437095
## 891 891 0.007428808
## 892 892 0.007420539
## 893 893 0.007412289
## 894 894 0.007404057
## 895 895 0.007395843
## 896 896 0.007387647
## 897 897 0.007379469
## 898 898 0.007371310
## 899 899 0.007363169
## 900 900 0.007355045
## 901 901 0.007346939
## 902 902 0.007338852
## 903 903 0.007330782
## 904 904 0.007322730
## 905 905 0.007314695
## 906 906 0.007306678
## 907 907 0.007298679
## 908 908 0.007290697
## 909 909 0.007282733
## 910 910 0.007274785
## 911 911 0.007266856
## 912 912 0.007258943
## 913 913 0.007251048
## 914 914 0.007243170
## 915 915 0.007235309
## 916 916 0.007227465
## 917 917 0.007219638
## 918 918 0.007211828
## 919 919 0.007204035
## 920 920 0.007196259
## 921 921 0.007188499
## 922 922 0.007180756
## 923 923 0.007173030
## 924 924 0.007165321
## 925 925 0.007157628
## 926 926 0.007149951
## 927 927 0.007142291
## 928 928 0.007134647
## 929 929 0.007127020
## 930 930 0.007119409
## 931 931 0.007111814
## 932 932 0.007104236
## 933 933 0.007096673
## 934 934 0.007089127
## 935 935 0.007081596
## 936 936 0.007074082
## 937 937 0.007066584
## 938 938 0.007059101
## 939 939 0.007051634
## 940 940 0.007044183
## 941 941 0.007036748
## 942 942 0.007029329
## 943 943 0.007021925
## 944 944 0.007014536
## 945 945 0.007007164
## 946 946 0.006999806
## 947 947 0.006992464
## 948 948 0.006985138
## 949 949 0.006977827
## 950 950 0.006970531
## 951 951 0.006963250
## 952 952 0.006955985
## 953 953 0.006948734
## 954 954 0.006941499
## 955 955 0.006934279
## 956 956 0.006927074
## 957 957 0.006919884
## 958 958 0.006912708
## 959 959 0.006905548
## 960 960 0.006898403
## 961 961 0.006891272
## 962 962 0.006884156
## 963 963 0.006877054
## 964 964 0.006869967
## 965 965 0.006862895
## 966 966 0.006855838
## 967 967 0.006848794
## 968 968 0.006841766
## 969 969 0.006834752
## 970 970 0.006827752
## 971 971 0.006820766
## 972 972 0.006813795
## 973 973 0.006806838
## 974 974 0.006799895
## 975 975 0.006792966
## 976 976 0.006786052
## 977 977 0.006779151
## 978 978 0.006772264
## 979 979 0.006765392
## 980 980 0.006758533
## 981 981 0.006751689
## 982 982 0.006744858
## 983 983 0.006738041
## 984 984 0.006731237
## 985 985 0.006724448
## 986 986 0.006717672
## 987 987 0.006710910
## 988 988 0.006704161
## 989 989 0.006697426
## 990 990 0.006690704
## 991 991 0.006683996
## 992 992 0.006677301
## 993 993 0.006670620
## 994 994 0.006663952
## 995 995 0.006657298
## 996 996 0.006650656
## 997 997 0.006644028
## 998 998 0.006637413
## 999 999 0.006630812
## 1000 1000 0.006624223
plot(out.M[,2]~ out.M[,1], type = "l", ylim=c(0,1), ylab = "f", xlab = "tiempo")
