João Wesley de Abreu Silva

539605

Método de Euller
i x y fi yi+1 0.025
0 1 5 -12 4.7
1 1.025 4.7 -11.05 4.42375
2 1.05 4.42375 -10.17125 4.169469
3 1.075 4.169469 -9.358406 3.935509
4 1.1 3.935509 -8.606526 3.720345
5 1.125 3.720345 -7.911036 3.52257
6 1.15 3.52257 -7.267709 3.340877
7 1.175 3.340877 -6.67263 3.174061
8 1.2 3.174061 -6.122183 3.021006
9 1.225 3.021006 -5.613019 2.880681
10 1.25 2.880681

22.4
t p
 1950 2560
1955 2780
1960 3040
1965 3350
1970 3710
1975 4090
1980 4450
1985 4850
1990 5280
1995 5690
2000 6080

k estimado analiticamente
t p k k médio
1950 2560 0.016489
1955 2780 0.017881
1960 3040 0.019421 0.018551
1965 3350 0.020414
1970 3710
22.1
Método de Euller
i t y fi yi+1
0 0 1 -1.5 0.625
1 0.25 0.625 -0.927734 0.393066
2 0.5 0.393066 -0.540466 0.25795
3 0.75 0.25795 -0.278102 0.188424
4 1 0.188424 -0.094212 0.164871
5 1.25 0.164871 0.074707 0.183548
6 1.5 0.183548 0.344153 0.269586
7 1.75 0.269586 1.040434 0.529695
8 2 0.529695 3.443016 1.390449

Método de Euller
i t y fi yi+1
0 0 1 -1.5 0.25
1 0.5 0.25 -0.34375 0.078125
2 1 0.078125 -0.039063 0.058594
3 1.5 0.058594 0.109863 0.113525
4 2 0.113525 0.737915 0.482483
 22.9

1.5

0.5

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5

-0.5

-1

-1.5

Column E Column I

Método de Euller
 i t f'i fi yi yi+1 0.1
0 0 -4 0 1 1
1 0.1 -4 -0.4 1 0.96
2 0.2 -3.84 -0.784 0.96 0.8816
3 0.3 -3.5264 -1.13664 0.8816 0.767936
4 0.4 -3.071744 -1.443814 0.767936 0.623555
5 0.5 -2.494218 -1.693236 0.623555 0.454231
6 0.6 -1.816924 -1.874929 0.454231 0.266738
7 0.7 -1.066952 -1.981624 0.266738 0.068576
8 0.8 -0.274303 -2.009054 0.068576 -0.13233
9 0.9 0.529319 -1.956122 -0.13233 -0.327942
10 1 1.311768 -1.824945 -0.327942 -0.510436
11 1.1 2.041746 -1.620771 -0.510436 -0.672514
12 1.2 2.690054 -1.351765 -0.672514 -0.80769
13 1.3 3.23076 -1.028689 -0.80769 -0.910559
14 1.4 3.642236 -0.664466 -0.910559 -0.977006
15 1.5 3.908023 -0.273664 -0.977006 -1.004372
16 1.6 4.017488 0.128085 -1.004372 -0.991563
17 1.7 3.966254 0.524711 -0.991563 -0.939092
18 1.8 3.75637 0.900348 -0.939092 -0.849058
19 1.9 3.39623 1.239971 -0.849058 -0.725061
20 2 2.900242 1.529995 -0.725061 -0.572061
21 2.1 2.288244 1.758819 -0.572061 -0.396179
22 2.2 1.584717 1.917291 -0.396179 -0.20445
23 2.3 0.8178 1.999071 -0.20445 -0.004543
24 2.4 0.018172 2.000888 -0.004543 0.195546
25 2.5 -0.782183 1.92267 0.195546 0.387813
26 2.6 -1.551251 1.767545 0.387813 0.564567
27 2.7 -2.258269 1.541718 0.564567 0.718739
28 2.8 -2.874956 1.254222 0.718739 0.844161
29 2.9 -3.376645 0.916558 0.844161 0.935817
30 3 -3.743268 0.542231 0.935817 0.99004
31 3.1 -3.960161 0.146215 0.99004 1.004662
32 3.2 -4.018646 -0.25565 1.004662 0.979097
33 3.3 -3.916386 -0.647289 0.979097 0.914368
34 3.4 -3.657471 -1.013036 0.914368 0.813064
35 3.5 -3.252257 -1.338261 0.813064 0.679238
36 3.6 -2.716952 -1.609957 0.679238 0.518242
37 3.7 -2.07297 -1.817254 0.518242 0.336517
38 3.8 -1.346068 -1.95186 0.336517 0.141331
39 3.9 -0.565324 -2.008393 0.141331 -0.059508
40 4 0.238033 -1.984589 -0.059508 -0.257967

22.8
Método de Euller
i t f'i fi yi
0 0 0 0 0
1 0.25 0 0 0
2 0.5 0 0 0
3 0.75 0 0 0
4 1 0 0 0
5 1.25 0 0 0
6 1.5 0 0 0
7 1.75 0 0 0
8 2 0 0 0
9 2.25 0 0 0
10 2.5 0 0 0
11 2.75 0 0 0
12 3 0 0 0
13 3.25 0 0 0
14 3.5 0 0 0
15 3.75 0 0 0
16 4 0 0 0
17 4.25 0 0 0
18 4.5 0 0 0
19 4.75 0 0 0
20 5 0 0 0
21 5.25 0 0 0
22 5.5 0 0 0
23 5.75 0 0 0
24 6 0 0 0
25 6.25 0 0 0
26 6.5 0 0 0
27 6.75 0 0 0
28 7 0 0 0
29 7.25 0 0 0
30 7.5 0 0 0
31 7.75 0 0 0
32 8 0 0 0
33 8.25 0 0 0
34 8.5 0 0 0
35 8.75 0 0 0
36 9 0 0 0
37 9.25 0 0 0
38 9.5 0 0 0
39 9.75 0 0 0
40 10 0 0 0
22.20

E 1.3E+08 Métod
dy/dz 0 i z
z 0 0 0
L 30 1 2
y 0 2 4
I 0.05 3 6
f(z) 0 4 8
5 10
6 12
7 14
8 16
9 18
10 20
11 22
12 24
13 26
14 28
15 30
16 32
17 34
18 36
19 38
20 40
21 42
22 44
23 46
24 48
25 50
 Método de Heun
i x y fi
0 1 5 -12
1 1.025 4.711875 -11.08563
2 1.05 4.445752 -10.23726
3 1.075 4.200043 -9.450129
4 1.1 3.973274 -8.719823
5 1.125 3.764079 -8.042236
6 1.15 3.571187 -7.413562
7 1.175 3.393424 -6.830271
8 1.2 3.229695 -6.289085
9 1.225 3.078989 -5.786967
10 1.25 2.940365

Método de Range Kutta

i x yi fi yi+1/2 fi+1/2
 0 1950 2560 47.49115 2678.728 49.6937
1 1955 2808.817 52.10701 2939.084 54.52363
2 1960 3081.817 57.1715 3224.746 59.82301
3 1965 3381.352 62.72824 3538.172 65.63745
4 1970 3709.999 68.82505 3882.062 72.01703
5 1975 4070.589 75.51444 4259.375 79.01666
6 1980 4466.226 82.854 4673.361 86.69661
7 1985 4900.317 90.90692 5127.585 95.12301
8 1990 5376.599 99.74253 5625.956 104.3684
9 1995 5899.173 109.4369 6172.765 114.5124
10 2000 6472.538 120.0735 6772.722 125.6423
11 2005 7101.63 131.744 7430.99 137.854
12 2010 7791.867 144.5487 8153.239 151.2526
13 2015 8549.19 158.598 8945.685 165.9535
14 2020 9380.121 174.0128 9815.153 182.0832
15 2025 10291.81 190.9258 10769.13 199.7806
16 2030 11292.12 209.4827 11815.82 219.1981
17 2035 12389.64 229.8432 12964.25 240.5028
18 2040 13593.84 252.1826 14224.3 263.8783
19 2045 14915.09 276.6932 15606.82 289.5257
20 2050 16364.74 303.5862 17123.71 317.6659
0.25 Método do Ponto Médio
i t yi fi
0 0 1 -1.5
1 0.5 0.536133 -0.737183
2 1 0.346471 -0.173235
3 1.5 0.415156 0.778417
4 2 1.591802

M
0.5 i t yi
0 0 1
1 0.5 0.481096
2 1 0.278038
3 1.5 0.351348
4 2 1.474755
4 4.5
 cos(2t) i t f'i
1 0 0 -4
0.980066578 1 0.1 -4.054128
0.921060994 2 0.2 -3.947165
0.825335615 3 0.3 -3.686523
0.696706709 4 0.4 -3.285533
0.540302306 5 0.5 -2.762801
0.362357754 6 0.6 -2.14137
0.169967143 7 0.7 -1.447723
-0.02919952 8 0.8 -0.710668
-0.22720209 9 0.9 0.039842
-0.41614684 10 1 0.773919
-0.58850112 11 1.1 1.462913
-0.73739372 12 1.2 2.080521
-0.85688875 13 1.3 2.603781
-0.94222234 14 1.4 3.013907
-0.9899925 15 1.5 3.296951
-0.99829478 16 1.6 3.444257
-0.96679819 17 1.7 3.452699
-0.89675842 18 1.8 3.324694
-0.79096771 19 1.9 3.067996
-0.65364362 20 2 2.695285
-0.49026082 21 2.1 2.22357
-0.30733287 22 2.2 1.673429
-0.11215253 23 2.3 1.068131
0.087498983 24 2.4 0.432659
0.283662185 25 2.5 -0.207308
0.468516671 26 2.6 -0.82642
0.634692876 27 2.7 -1.400653
0.775565879 28 2.8 -1.908231
 0.885519517 29 2.9 -2.330439
0.960170287 30 3 -2.652303
0.996542097 31 3.1 -2.863106
0.993184919 32 3.2 -2.956725
0.950232592 33 3.3 -2.931787
0.86939749 34 3.4 -2.791624
0.753902254 35 3.5 -2.544048
0.608351315 36 3.6 -2.20096
0.438547328 37 3.7 -1.777792
0.251259843 38 3.8 -1.292833
0.053955421 39 3.9 -0.766448
-0.14550003 40 4 -0.220233

yi+1 0.25
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
 Método de Euller
yi f(z) f'i fi yi+1 2
0 0 0 0 0
0 50.00990395 0.003016 0.006032 0.012064
0.012064 68.08251897 0.00354 0.013113 0.038289
0.038289 73.1258232 0.00324 0.019593 0.077474
0.077474 72.20261163 0.002688 0.024969 0.127412
0.127412 68.45561587 0.002106 0.029182 0.185775
0.185775 63.43467729 0.001581 0.032344 0.250462
0.250462 57.95126413 0.001141 0.034626 0.319714
0.319714 52.44248181 0.000791 0.036207 0.392129
0.392129 47.14344186 0.000522 0.037252 0.466632
0.466632 42.1755421 0.000324 0.037901 0.542433
0.542433 37.59444452 0.000185 0.038271 0.618975
0.618975 33.4173547 9.254E-05 0.038456 0.695886
0.695886 29.63906832 3.648E-05 0.038529 0.772944
0.772944 26.24164489 8.074E-06 0.038545 0.850034
0.850034 23.20033427 0 0.038545 0.927124
0.927124 20.48723529 6.304E-06 0.038558 1.004239
1.004239 18.07355063 2.224E-05 0.038602 1.081443
1.081443 15.93095765 4.412E-05 0.03869 1.158823
1.158823 14.03241594 6.908E-05 0.038828 1.23648
1.23648 12.35261355 9.502E-05 0.039018 1.314517
1.314517 10.86818141 0.00012 0.039259 1.393035
1.393035 9.5577598 0.000144 0.039547 1.47213
1.47213 8.401972197 0.000165 0.039878 1.551887
1.551887 7.383342607 0.000184 0.040246 1.632379
1.632379 6.486180609 0.0002 0.040646 1.71367
do de Heun
yi+1 fi+1 yi+1,corrigido 0.025
4.7 -11.05 4.711875
4.434734 -10.2042 4.445752 Mét
4.189821 -9.419462 4.200043 i x
3.96379 -8.69137 3.973274 0 1
3.755279 -8.015837 3.764079 1 1.025
3.563023 -7.389068 3.571187 2 1.05
3.385848 -6.807545 3.393424 3 1.075
3.222667 -6.268 3.229695 4 1.1
3.072468 -5.767404 3.078989 5 1.125
2.934315 -5.302945 2.940365 6 1.15
7 1.175
8 1.2
9 1.225
10 1.25

e Range Kutta
yi+1/2,c fi+1/2,c yi+1 fi+1 yi+1,rk 5
2684.234 49.79585 2808.979 52.11002 2808.817
2945.126 54.63571 3081.995 57.17481 3081.817
3231.375 59.94598 3381.547 62.73186 3381.352
3545.445 65.77238 3710.214 68.82903 3709.999
3890.042 72.16506 4070.825 75.51881 4070.589
4268.131 79.17908 4466.485 82.85879 4466.226
4682.968 86.87482 4900.601 90.91217 4900.317
5138.125 95.31854 5376.91 99.7483 5376.599
5637.52 104.5829 5899.514 109.4432 5899.173
6185.454 114.7478 6472.912 120.0805 6472.538
6786.644 125.9006 7102.041 131.7516 7101.63
7446.265 138.1374 7792.317 144.5571 7791.867
8169.998 151.5635 8549.684 158.6072 8549.19
8964.074 166.2946 9380.663 174.0229 9380.121
9835.329 182.4575 10292.41 190.9368 10291.81
10791.27 200.1912 11292.77 209.4948 11292.12
11840.11 219.6487 12390.36 229.8564 12389.64
12990.9 240.9972 13594.63 252.1971 13593.84
14253.54 264.4207 14915.95 276.7092 14915.09
15638.9 290.1209 16365.69 303.6037 16364.74
17158.91 318.3189 17956.34 333.1122 17955.3
do Ponto Médio
yi+1/2 fi+1/2 yi+1 0.5
0.625 -0.927734 0.536133
0.351837 -0.379324 0.346471
0.303162 0.13737 0.415156
0.60976 2.353292 1.591802

Método de Range Kutta

fi yi+1/2 fi+1/2 yi+1/2,c fi+1/2,c yi+1 fi+1
-1.5 0.625 -0.927734 0.768066 -1.140099 0.429951 -0.591182
-0.661507 0.315719 -0.340385 0.396 -0.426937 0.481096 -0.240548
-0.139019 0.243283 0.110238 0.305597 0.138474 0.278038 0.521321
0.658778 0.516043 1.991602 0.849249 3.277568 0.351348 2.283763
 RK e Cos2t
1.5

0.5

0
-1.5 -1 -0.5 0 0.5 1 1.5

-0.5

-1

Column O

Método de Range Kutta

fi yi yi+1/2 f'i+1/2 fi+1/2 yi+1/2,c f'i+1/2,c
0 1 1 -4 -0.2 0.99 -3.96
-0.405413 1.013532 0.993261 -3.973045 -0.604065 0.983329 -3.933315
-0.800129 0.986791 0.946785 -3.787139 -0.989486 0.937317 -3.749268
-1.168782 0.921631 0.863192 -3.452767 -1.34142 0.85456 -3.418239
-1.497335 0.821383 0.746516 -2.986066 -1.646638 0.739051 -2.956205
-1.773615 0.6907 0.602019 -2.408078 -1.894019 0.595999 -2.383997
-1.987752 0.535342 0.435955 -1.743819 -2.074943 0.431595 -1.726381
-2.132524 0.361931 0.255304 -1.021218 -2.183585 0.252751 -1.011006
-2.203591 0.177667 0.067487 -0.26995 -2.217089 0.066813 -0.26725
-2.199607 -0.009961 -0.119941 0.479764 -2.175619 -0.118741 0.474966
-2.122215 -0.19348 -0.29959 1.198362 -2.062297 -0.296595 1.186378
-1.975924 -0.365728 -0.464524 1.858097 -1.883019 -0.459879 1.839516
-1.767871 -0.52013 -0.608524 2.434096 -1.646167 -0.602439 2.409755
-1.507493 -0.650945 -0.72632 2.90528 -1.362229 -0.719057 2.876227
-1.206103 -0.753477 -0.813782 3.255127 -1.043346 -0.805644 3.222576
-0.876408 -0.824238 -0.868058 3.472232 -0.702796 -0.859377 3.43751
-0.531982 -0.861064 -0.887663 3.550653 -0.354449 -0.878787 3.515147
-0.186712 -0.863175 -0.87251 3.490042 -0.01221 -0.863785 3.455141
0.145757 -0.831174 -0.823886 3.295543 0.310535 -0.815647 3.262587
0.452557 -0.766999 -0.744371 2.977485 0.601431 -0.736928 2.94771
0.722086 -0.673821 -0.637717 2.550868 0.849629 -0.63134 2.52536
0.944443 -0.555893 -0.50867 2.034681 1.046177 -0.503584 2.014335
1.111786 -0.418357 -0.362768 1.451072 1.184339 -0.35914 1.436562
1.218599 -0.267033 -0.206103 0.824411 1.259819 -0.204042 0.816167
1.261865 -0.108165 -0.045071 0.180286 1.270879 -0.044621 0.178483
1.241134 0.051827 0.113884 -0.455535 1.218357 0.112745 -0.450979
1.158492 0.206605 0.26453 -1.058118 1.105586 0.261884 -1.047537
1.018426 0.350163 0.401085 -1.604338 0.93821 0.397074 -1.588295
0.827603 0.477058 0.518438 -2.073752 0.723916 0.513254 -2.053014
 0.594559 0.58261 0.612338 -2.449351 0.472092 0.606214 -2.424857
0.329329 0.663076 0.679542 -2.718169 0.193421 0.672747 -2.690987
0.043019 0.715776 0.717927 -2.871709 -0.100567 0.710748 -2.842992
-0.252654 0.739181 0.726549 -2.906195 -0.397964 0.719283 -2.877133
-0.545833 0.732947 0.705655 -2.822621 -0.686964 0.698599 -2.794395
-0.824995 0.697906 0.656656 -2.626625 -0.956326 0.65009 -2.600359
-1.0794 0.636012 0.582042 -2.328169 -1.195808 0.576222 -2.304887
-1.299496 0.55024 0.485265 -1.941061 -1.396549 0.480413 -1.92165
-1.477275 0.444448 0.370584 -1.482337 -1.551392 0.366878 -1.467514
-1.606558 0.323208 0.24288 -0.971521 -1.655134 0.240452 -0.961806
-1.683203 0.191612 0.107452 -0.429807 -1.704694 0.106377 -0.425509
-1.705227 0.055058 -0.030203 0.120812 -1.699186 -0.029901 0.119604
 Método do Ponto Médio
yi fi yi+1/2 fi+1/2 yi+1 0.025
5 -12 4.85 -11.525 4.711875
4.711875 -11.08563 4.573305 -10.64491 4.445752
4.4457521484 -10.23726 4.317786 -9.828359 4.200043
4.2000431652 -9.450129 4.081917 -9.07075 3.973274
3.9732744242 -8.719823 3.864277 -8.36783 3.764079
3.7640786767 -8.042236 3.663551 -7.715652 3.571187
3.5711873723 -7.413562 3.478518 -7.110554 3.393424
3.3934235338 -6.830271 3.308045 -6.549135 3.229695
3.2296951475 -6.289085 3.151082 -6.028245 3.078989
3.078989029 -5.786967 3.006652 -5.544956 2.940365
2.9403651335
yi+1,rk 0.5
0.4810959895
0.2780377279
0.3513481034
1.4747548165
1.5
fi+1/2,c yi+1 f'i+1 fi+1 yi+1,rk
-0.198 0.9802 -3.9208 -0.39208 1.013532
-0.602078549 0.953324 -3.813297 -0.786742 0.986791
-0.987592712 0.888032 -3.552128 -1.155342 0.921631
-1.339693598 0.787661 -3.150646 -1.483846 0.821383
-1.645145186 0.656869 -2.627475 -1.760082 0.6907
-1.892814825 0.501419 -2.005675 -1.974182 0.535342
-2.074070996 0.327935 -1.311741 -2.118926 0.361931
-2.183074477 0.143623 -0.574493 -2.189973 0.177667
-2.21695353 -0.044028 0.176113 -2.18598 -0.009961
-2.175858491 -0.227546 0.910186 -2.108588 -0.19348
-2.062895995 -0.399769 1.599077 -1.962307 -0.365728
-1.883947816 -0.554123 2.216492 -1.754274 -0.52013
-1.647383756 -0.684869 2.739475 -1.493924 -0.650945
-1.363682022 -0.787313 3.149254 -1.192568 -0.753477
-1.044973892 -0.857974 3.431896 -0.862913 -0.824238
-0.704532153 -0.894691 3.578763 -0.518531 -0.861064
-0.356224612 -0.896687 3.586747 -0.173307 -0.863175
-0.013954945 -0.86457 3.458281 0.159116 -0.831174
0.3088868029 -0.800285 3.20114 0.465871 -0.766999
0.599942561 -0.707005 2.828019 0.735359 -0.673821
0.8483535769 -0.588986 2.355944 0.95768 -0.555893
1.0451593315 -0.451377 1.805506 1.124993 -0.418357
1.1836136286 -0.299996 1.199984 1.231784 -0.267033
1.2594070253 -0.141092 0.564368 1.275035 -0.108165
1.2707886771 0.018914 -0.075657 1.254299 0.051827
1.2185847812 0.173685 -0.694742 1.17166 0.206605
1.1061149376 0.317216 -1.268866 1.031605 0.350163
0.939011749 0.444064 -1.776258 0.840801 0.477058
0.7249526983 0.549553 -2.198212 0.607782 0.58261
0.4733166219 0.629941 -2.519766 0.342583 0.663076
0.194779832 0.682554 -2.730215 0.056308 0.715776
-0.099130976 0.705863 -2.823453 -0.239327 0.739181
-0.396510528 0.69953 -2.798121 -0.532466 0.732947
-0.685552378 0.664392 -2.657566 -0.811589 0.697906
-0.955012956 0.602405 -2.409619 -1.065957 0.636012
-1.194644216 0.516548 -2.066191 -1.286019 0.55024
-1.395578419 0.410682 -1.642729 -1.463769 0.444448
-1.550650817 0.289383 -1.157532 -1.593028 0.323208
-1.654648728 0.157743 -0.630974 -1.669656 0.191612
-1.704478631 0.021164 -0.084656 -1.691669 0.055058
-1.699246306 -0.114866 0.459465 -1.65928 -0.080965
 Método de Range Kutta
i x yi fi yi+1/2 fi+1/2
0 1 5 -12 4.7 -11.05
1 1.05 4.445214 -10.23564 4.189323 -9.417969
2 1.1 3.972348 -8.717044 3.754422 -8.013266
3 1.15 3.569991 -7.409974 3.384742 -6.804226
4 1.2 3.228322 -6.284967 3.071198 -5.763595
5 1.25 2.938888
de Range Kutta
yi+1/2,c fi+1/2,c yi+1 fi+1 yi+1,rk 0.05
4.72375 -11.12125 4.443938 -10.23181 4.445214
4.209765 -9.479295 3.971249 -8.713748 3.972348
3.772016 -8.066049 3.569046 -7.407137 3.569991
3.399886 -6.849657 3.227508 -6.282525 3.228322
3.084233 -5.802698 2.938188 -5.314563 2.938888
0.1