Ministério da Educação

Universidade Tecnológica Federal do Paraná

Pró-reitora de Graduação e Educação Profissional

PROJETO 3 – RESOLUÇÃO NUMÉRICA DE

EDO DE SEGUNDA ORDEM

GUILHERME LUIZ GUERRA VICENTE

Professor: Alireza Mohebi Ashtiani

Londrina
2018
1. INTRODUÇÃO

A disciplina de Cálculo Numérico tem como objetivo principal buscar respostas

aproximadas de problemas de alta complexidade (polinômios, expressões algébricas,
sistemas de equações, etc.), já que pelos métodos tradicionais encontrar uma solução
exata para tais problemas é: ou extremamente difícil e requer mais base e estudo na
área de matemática, ou impossível encontrar devido a sua complexidade.
O presente relatório a seguir trata exatamente dessa questão. No caso, uma
aproximação de uma resolução de EDO de segunda ordem, realizando uma
comparação entre os métodos de Euler, Taylor de 2ªOrdem, Taylor de 4ªOrdem,
Runge-Kutta de 2ªOrdem e Runge-Kutta de 4ªOrdem
O resultado numérico de uma EDO se dá por uma tabela de pontos que é gerada
baseada nos métodos apresentados e a partir dela é possível deduzir a forma
equacionaria da resposta.

2. OBJETIVO(S)

 Encontrar a solução numérica para a EDO;

 Compreender os métodos;
 Analisar as diferenças entre os métodos;
 Compreender a utilização dos softwares MatLab e Octave;
3. PROCEDIMENTOS EXPERIMENTAIS

Para início do projeto foi escolhido uma EDO de Segunda Ordem a ser resolvida
utilizando os códigos aprendidos em sala de aula, com as devidas adaptações para
sua utilização na EDO escolhida. A EDO e seus respectivos valores iniciais são:

𝑦 (𝑥) + 𝑦 (𝑥) = 𝑥𝑒
𝑦 (0) = 1 , 𝑥 ∈ [0,3]
𝑦(0) = 1

Realizando uma redução de ordem na presente equação obtemos:

𝑧 = 𝑦′′
𝑧 = 𝑥𝑒 − 𝑧
𝑧(0) = 1
𝑦(0) = 1

Utilizando-se então do software online WolframAlpha a solução exata da EDO

encontrada foi a seguinte:

A solução exata da EDO foi encontrada para fins de comparação dos métodos
analisados.
4. RESULTADOS

A partir então do software Octave foram gerados os seguintes gráficos e tabelas

para diversos números de subintervalos.

4.1. NÚMERO DE SUBINTERVALOS = 4

VALORES DE Y PARA CADA MÉTODO
EULER TAYLOR 2 TAYLOR 4 RK2 RK4 EXATO
1 1 1 1 1 1
1,75 1,4688 1,6094 1,4688 1,609 1,6157
1,9375 2,2962 2,7017 2,4434 2,7043 2,7211
2,8775 5,2265 6,7716 5,8662 6,3844 6,4262
6,8939 14,737 18,1581 16,7319 17,8967 18,0019

ERROS DE CADA MÉTODO

EULER TAYLOR 2 TAYLOR 4 RK2 RK4
0 0 0 0 0
0,13433 0,14692 0,00629 0,14692 0,00667
0,78359 0,42492 0,01939 0,27766 0,01676
3,54867 1,19965 0,34547 0,55998 0,04179
11,10801 3,26493 0,15616 1,27 0,10519
4.2. NÚMERO DE SUBINTERVALOS = 10
VALORES DE Y PARA CADA MÉTODO
EULER TAYLOR 2 TAYLOR 4 RK2 RK4 EXATO
1 1 1 1 1 1
1,3 1,255 1,264 1,255 1,264 1,2641
1,51 1,4747 1,4939 1,4787 1,4939 1,494
1,6934 1,7204 1,8762 1,7336 1,7537 1,7539
1,9203 2,071 2,4017 2,1004 2,1252 2,1255
2,2782 2,6348 3,1842 2,6905 2,7207 2,7211
2,8874 3,5667 4,3388 3,6628 3,7002 3,7008
3,9189 5,0913 6,0612 5,2483 5,2959 5,2968
5,6209 7,5362 8,6298 7,7837 7,8457 7,847
8,3558 11,3789 12,4589 11,7593 11,842 11,8438
12,6512 17,313 18,1523 17,887 17,9992 18,0019
 ERROS DE CADA MÉTODO
EULER TAYLOR 2 TAYLOR 4 RK2 RK4
0 0 0 0 0
0,03594 0,00906 0,00006 0,00906 0,00006
0,01597 0,01936 0,00014 0,01534 0,00012
0,06046 0,03349 0,12233 0,02029 0,0002
0,20524 0,05448 0,27617 0,02505 0,00029
0,44284 0,08625 0,46311 0,03062 0,00042
0,81341 0,1341 0,63801 0,03803 0,00062
1,37792 0,20546 0,7644 0,04844 0,00089
2,22612 0,31079 0,78277 0,06331 0,0013
3,48808 0,46493 0,6151 0,08458 0,00188
5,35076 0,68894 0,15035 0,11488 0,00271
4.3. NÚMERO DE SUBINTERVALOS = 50
VALORES DE Y PARA CADA MÉTODO
EULER TAYLOR 2 TAYLOR 4 RK2 RK4 EXATO
1 1 1 1 1 1
1,06 1,0582 1,0583 1,0582 1,0583 1,0583
1,1164 1,1132 1,1134 1,1132 1,1134 1,1134
1,1696 1,1655 1,1724 1,1656 1,1657 1,1657
1,2202 1,2155 1,2349 1,2156 1,2158 1,2158
1,2685 1,2637 1,301 1,2638 1,2641 1,2641
1,3149 1,3105 1,3711 1,3106 1,3109 1,3109
1,3601 1,3563 1,4454 1,3564 1,3568 1,3568
1,4044 1,4017 1,5241 1,4019 1,4023 1,4023
1,4483 1,4472 1,6076 1,4474 1,4479 1,4479
1,4924 1,4932 1,6962 1,4935 1,494 1,494
1,5372 1,5404 1,7902 1,5408 1,5413 1,5413
1,5833 1,5893 1,89 1,5898 1,5903 1,5903
1,6311 1,6405 1,996 1,6411 1,6417 1,6417
1,6815 1,6946 2,1086 1,6953 1,696 1,696
1,7349 1,7524 2,2282 1,7533 1,7539 1,7539
1,7921 1,8146 2,3554 1,8156 1,8162 1,8162
1,8539 1,8819 2,4906 1,883 1,8837 1,8837
1,921 1,9551 2,6344 1,9564 1,9571 1,9571
1,9942 2,0352 2,7874 2,0367 2,0374 2,0374
2,0745 2,123 2,9502 2,1247 2,1255 2,1255
 2,1628 2,2197 3,1235 2,2216 2,2224 2,2224
2,2602 2,3262 3,3081 2,3284 2,3292 2,3292
2,3677 2,4437 3,5047 2,4462 2,447 2,447
2,4865 2,5736 3,7141 2,5764 2,5772 2,5772
2,618 2,7171 3,9373 2,7202 2,7211 2,7211
2,7635 2,8758 4,1752 2,8792 2,8801 2,8801
2,9244 3,0511 4,4289 3,0549 3,0558 3,0558
3,1024 3,2448 4,6994 3,249 3,2499 3,2499
3,2992 3,4586 4,988 3,4633 3,4643 3,4643
3,5166 3,6947 5,296 3,6998 3,7008 3,7008
3,7567 3,9549 5,6247 3,9606 3,9617 3,9617
4,0215 4,2418 5,9756 4,248 4,2491 4,2491
4,3135 4,5577 6,3503 4,5645 4,5657 4,5657
4,6352 4,9052 6,7504 4,9127 4,9139 4,9139
4,9891 5,2873 7,1779 5,2955 5,2968 5,2968
5,3783 5,707 7,6346 5,716 5,7173 5,7173
5,8059 6,1677 8,1227 6,1775 6,1789 6,1789
6,2752 6,6729 8,6444 6,6837 6,685 6,6851
6,79 7,2265 9,2022 7,2383 7,2397 7,2397
7,3541 7,8327 9,7986 7,8455 7,847 7,847
7,9719 8,496 10,4364 8,5099 8,5115 8,5115
8,6478 9,2212 11,1186 9,2364 9,238 9,238
9,3868 10,0135 11,8484 10,03 10,0317 10,0317
10,1942 10,8786 12,6294 10,8965 10,8983 10,8983
11,0758 11,8225 13,4651 11,8419 11,8438 11,8438
12,0376 12,8517 14,3595 12,8728 12,8748 12,8748
13,0864 13,9732 15,317 13,9961 13,9982 13,9982
14,2292 15,1946 16,342 15,2195 15,2217 15,2217
15,4738 16,5241 17,4395 16,551 16,5533 16,5533
16,8284 17,9703 18,6148 17,9995 18,0019 18,0019

ERROS DE CADA MÉTODO

EULER TAYLOR 2 TAYLOR 4 RK2 RK4
0 0 0 0 0
0,00173 0,00007 0 0,00007 0
0,00302 0,00014 0 0,00014 0
0,0039 0,00022 0,00668 0,0002 0
0,00435 0,00029 0,01904 0,00025 0
0,0044 0,00037 0,03695 0,0003 0
0,00404 0,00045 0,06018 0,00035 0
0,00327 0,00054 0,08853 0,00039 0
0,00208 0,00063 0,12178 0,00043 0
0,00046 0,00073 0,15973 0,00047 0
0,00159 0,00083 0,20217 0,0005 0
0,00409 0,00094 0,2489 0,00054 0
0,00706 0,00106 0,2997 0,00057 0
0,01052 0,00119 0,35436 0,0006 0
0,01449 0,00134 0,41265 0,00062 0
0,01901 0,00149 0,47434 0,00065 0
0,0241 0,00165 0,53919 0,00067 0
0,0298 0,00183 0,60693 0,0007 0
0,03614 0,00203 0,6773 0,00072 0
0,04319 0,00224 0,75001 0,00074 0
0,05097 0,00247 0,82474 0,00077 0
0,05955 0,00272 0,90116 0,00079 0
0,06898 0,00299 0,97892 0,00081 0
0,07934 0,00329 1,05763 0,00084 0
0,09068 0,00361 1,13687 0,00086 0
0,10309 0,00395 1,21619 0,00089 0
0,11664 0,00433 1,29509 0,00092 0
0,13143 0,00474 1,37304 0,00094 0
0,14755 0,00518 1,44947 0,00097 0
0,16511 0,00566 1,52374 0,001 0
0,18423 0,00617 1,59516 0,00104 0
0,20501 0,00673 1,663 0,00107 0
0,2276 0,00734 1,72643 0,00111 0
0,25214 0,00799 1,78458 0,00115 0
0,27878 0,0087 1,83649 0,00119 0
0,30768 0,00947 1,8811 0,00124 0
0,33903 0,01029 1,91729 0,00129 0
0,37301 0,01119 1,94383 0,00134 0
0,40984 0,01215 1,95937 0,00139 0
0,44972 0,01319 1,96246 0,00145 0
0,49291 0,01431 1,95153 0,00152 0
0,53966 0,01552 1,92485 0,00158 0
0,59024 0,01682 1,88058 0,00166 0
0,64496 0,01823 1,81671 0,00173 0
0,70412 0,01974 1,73104 0,00182 0
0,76807 0,02138 1,62123 0,0019 0
0,83718 0,02313 1,48472 0,002 0
0,91184 0,02502 1,31874 0,0021 0
0,99248 0,02706 1,12031 0,00221 0
1,07955 0,02925 0,88618 0,00232 0
1,17355 0,03161 0,61288 0,00245 0
4.4. NÚMERO DE SUBINTERVALOS = 100
VALORES DE Y PARA CADA MÉTODO
EULER TAYLOR 2 TAYLOR 4 RK2 RK4 EXATO
1 1 1 1 1 1
1,03 1,0295 1,0296 1,0295 1,0296 1,0296
1,0591 1,0583 1,0583 1,0583 1,0583 1,0583
1,0874 1,0862 1,0879 1,0862 1,0862 1,0862
1,1148 1,1133 1,1184 1,1133 1,1134 1,1134
1,1415 1,1398 1,1499 1,1398 1,1399 1,1399
1,1676 1,1657 1,1822 1,1657 1,1657 1,1657
1,193 1,191 1,2155 1,191 1,191 1,191
1,2179 1,2158 1,2498 1,2158 1,2158 1,2158
1,2422 1,2401 1,2851 1,2401 1,2401 1,2401
1,2661 1,264 1,3215 1,264 1,2641 1,2641
1,2896 1,2875 1,3589 1,2875 1,2876 1,2876
1,3127 1,3108 1,3974 1,3108 1,3109 1,3109
1,3356 1,3338 1,4371 1,3339 1,334 1,334
1,3583 1,3567 1,478 1,3567 1,3568 1,3568
1,3807 1,3795 1,5201 1,3795 1,3796 1,3796
1,4031 1,4022 1,5635 1,4022 1,4023 1,4023
1,4255 1,4249 1,6081 1,425 1,4251 1,4251
1,4479 1,4477 1,6541 1,4478 1,4479 1,4479
1,4703 1,4707 1,7015 1,4707 1,4708 1,4708
1,493 1,4938 1,7503 1,4939 1,494 1,494
1,5158 1,5173 1,8006 1,5174 1,5175 1,5175
1,539 1,5411 1,8524 1,5412 1,5413 1,5413
1,5625 1,5653 1,9057 1,5654 1,5656 1,5656
1,5865 1,5901 1,9607 1,5902 1,5903 1,5903
1,611 1,6154 2,0173 1,6155 1,6157 1,6157
1,636 1,6414 2,0757 1,6415 1,6417 1,6417
1,6618 1,6681 2,1358 1,6683 1,6684 1,6684
1,6883 1,6956 2,1978 1,6958 1,696 1,696
1,7157 1,7241 2,2616 1,7243 1,7244 1,7244
1,744 1,7535 2,3274 1,7538 1,7539 1,7539
1,7733 1,7841 2,3952 1,7843 1,7845 1,7845
1,8037 1,8158 2,4651 1,8161 1,8162 1,8162
1,8354 1,8488 2,5371 1,8491 1,8493 1,8493
1,8683 1,8832 2,6114 1,8835 1,8837 1,8837
1,9027 1,9191 2,6879 1,9194 1,9196 1,9196
1,9385 1,9566 2,7667 1,9569 1,9571 1,9571
1,976 1,9958 2,848 1,9962 1,9964 1,9964
2,0152 2,0368 2,9318 2,0372 2,0374 2,0374
2,0563 2,0798 3,0182 2,0802 2,0804 2,0804
2,0994 2,1249 3,1073 2,1253 2,1255 2,1255
2,1445 2,1721 3,1991 2,1726 2,1728 2,1728
2,1919 2,2217 3,2938 2,2222 2,2224 2,2224
 2,2417 2,2737 3,3914 2,2743 2,2745 2,2745
2,294 2,3284 3,492 2,329 2,3292 2,3292
2,3489 2,3858 3,5957 2,3864 2,3866 2,3866
2,4066 2,4462 3,7027 2,4468 2,447 2,447
2,4672 2,5096 3,8131 2,5103 2,5105 2,5105
2,531 2,5763 3,9269 2,577 2,5772 2,5772
2,5981 2,6464 4,0442 2,6471 2,6474 2,6474
2,6686 2,7201 4,1652 2,7209 2,7211 2,7211
2,7428 2,7976 4,29 2,7984 2,7986 2,7986
2,8208 2,879 4,4187 2,8799 2,8801 2,8801
2,9028 2,9646 4,5515 2,9655 2,9658 2,9658
2,989 3,0546 4,6885 3,0556 3,0558 3,0558
3,0797 3,1492 4,8297 3,1502 3,1505 3,1505
3,175 3,2486 4,9754 3,2497 3,2499 3,2499
3,2752 3,3531 5,1258 3,3542 3,3545 3,3545
3,3805 3,4628 5,2808 3,464 3,4643 3,4643
3,4911 3,5781 5,4408 3,5794 3,5796 3,5796
3,6073 3,6993 5,6058 3,7006 3,7008 3,7008
3,7294 3,8264 5,7761 3,8278 3,8281 3,8281
3,8577 3,96 5,9518 3,9614 3,9617 3,9617
3,9924 4,1001 6,1331 4,1017 4,1019 4,1019
4,1337 4,2473 6,3201 4,2489 4,2491 4,2491
4,2821 4,4017 6,5131 4,4033 4,4036 4,4036
4,4379 4,5636 6,7123 4,5654 4,5657 4,5657
4,6013 4,7335 6,9178 4,7354 4,7357 4,7357
4,7727 4,9117 7,1299 4,9137 4,9139 4,9139
4,9524 5,0985 7,3488 5,1006 5,1008 5,1008
5,1409 5,2944 7,5747 5,2965 5,2968 5,2968
5,3385 5,4996 7,8079 5,5018 5,5021 5,5021
5,5456 5,7147 8,0485 5,717 5,7173 5,7173
5,7626 5,94 8,2968 5,9424 5,9428 5,9428
5,99 6,176 8,5532 6,1786 6,1789 6,1789
6,2281 6,4232 8,8178 6,4258 6,4262 6,4262
6,4775 6,682 9,0909 6,6847 6,685 6,6851
6,7387 6,9528 9,3729 6,9557 6,9561 6,9561
7,012 7,2364 9,6639 7,2394 7,2397 7,2397
7,2981 7,533 9,9644 7,5362 7,5365 7,5365
7,5975 7,8434 10,2746 7,8467 7,847 7,847
7,9108 8,168 10,5949 8,1715 8,1718 8,1718
8,2384 8,5076 10,9255 8,5112 8,5115 8,5115
8,581 8,8626 11,2669 8,8663 8,8667 8,8667
8,9393 9,2337 11,6193 9,2376 9,238 9,238
9,3139 9,6217 11,9832 9,6257 9,6261 9,6261
9,7054 10,0271 12,359 10,0313 10,0317 10,0317
10,1146 10,4507 12,747 10,4552 10,4556 10,4556
10,5421 10,8933 13,1476 10,8979 10,8983 10,8983
10,9888 11,3556 13,5614 11,3604 11,3608 11,3608
11,4553 11,8384 13,9886 11,8434 11,8438 11,8438
11,9426 12,3425 14,4298 12,3478 12,3482 12,3482
12,4514 12,8689 14,8855 12,8743 12,8748 12,8748
12,9826 13,4184 15,3561 13,424 13,4245 13,4245
13,5371 13,9919 15,8421 13,9977 13,9982 13,9982
14,1158 14,5903 16,344 14,5965 14,5969 14,5969
14,7198 15,2148 16,8625 15,2212 15,2217 15,2217
15,35 15,8663 17,3979 15,8729 15,8734 15,8734
16,0075 16,5459 17,951 16,5528 16,5533 16,5533
16,6933 17,2547 18,5223 17,2619 17,2624 17,2624
17,4086 17,9939 19,1124 18,0014 18,0019 18,0019

ERROS DE CADA MÉTODO

EULER TAYLOR 2 TAYLOR 4 RK2 RK4
0 0 0 0 0
0,00044 0,00001 0 0,00001 0
0,00083 0,00002 0 0,00002 0
0,00116 0,00003 0,00173 0,00003 0
0,00144 0,00004 0,00507 0,00003 0
0,00167 0,00005 0,00999 0,00004 0
0,00185 0,00005 0,01647 0,00005 0
0,00197 0,00006 0,02447 0,00006 0
0,00205 0,00007 0,03398 0,00006 0
0,00207 0,00008 0,04496 0,00007 0
0,00205 0,00009 0,0574 0,00007 0
0,00197 0,0001 0,07126 0,00008 0
0,00184 0,00011 0,08653 0,00009 0
0,00166 0,00012 0,10317 0,00009 0
0,00143 0,00014 0,12117 0,0001 0
0,00114 0,00015 0,1405 0,0001 0
0,00081 0,00016 0,16114 0,00011 0
0,00042 0,00017 0,18306 0,00011 0
0,00002 0,00018 0,20624 0,00011 0
0,00052 0,0002 0,23065 0,00012 0
0,00107 0,00021 0,25627 0,00012 0
0,00168 0,00022 0,28308 0,00013 0
0,00235 0,00024 0,31105 0,00013 0
0,00308 0,00025 0,34015 0,00013 0
0,00386 0,00027 0,37037 0,00014 0
0,00471 0,00029 0,40167 0,00014 0
0,00562 0,0003 0,43404 0,00014 0
0,0066 0,00032 0,46743 0,00015 0
0,00764 0,00034 0,50183 0,00015 0
0,00875 0,00036 0,5372 0,00015 0
0,00993 0,00038 0,57352 0,00016 0
0,01118 0,0004 0,61076 0,00016 0
0,0125 0,00042 0,64889 0,00016 0
0,01391 0,00044 0,68787 0,00017 0
0,01539 0,00047 0,72768 0,00017 0
0,01695 0,00049 0,76827 0,00017 0
0,0186 0,00051 0,80962 0,00017 0
0,02034 0,00054 0,85168 0,00018 0
0,02217 0,00057 0,89442 0,00018 0
0,02409 0,0006 0,93781 0,00018 0
0,0261 0,00063 0,98179 0,00018 0
0,02822 0,00066 1,02632 0,00019 0
0,03044 0,00069 1,07137 0,00019 0
0,03277 0,00073 1,11688 0,00019 0
0,03522 0,00076 1,16281 0,00019 0
0,03777 0,0008 1,2091 0,0002 0
0,04045 0,00084 1,25571 0,0002 0
0,04326 0,00088 1,30257 0,0002 0
0,04619 0,00092 1,34964 0,00021 0
0,04926 0,00096 1,39684 0,00021 0
0,05246 0,001 1,44412 0,00021 0
0,05581 0,00105 1,49141 0,00021 0
0,05932 0,0011 1,53864 0,00022 0
0,06297 0,00115 1,58575 0,00022 0
0,06679 0,0012 1,63264 0,00022 0
0,07078 0,00126 1,67926 0,00023 0
0,07494 0,00132 1,7255 0,00023 0
0,07929 0,00138 1,77129 0,00023 0
0,08382 0,00144 1,81654 0,00024 0
0,08855 0,0015 1,86115 0,00024 0
0,09348 0,00157 1,90502 0,00024 0
0,09863 0,00164 1,94805 0,00025 0
0,10399 0,00171 1,99014 0,00025 0
0,10958 0,00179 2,03116 0,00026 0
0,11541 0,00187 2,071 0,00026 0
0,12148 0,00195 2,10954 0,00026 0
0,12781 0,00203 2,14664 0,00027 0
0,13441 0,00212 2,18218 0,00027 0
0,14128 0,00221 2,21601 0,00028 0
0,14843 0,00231 2,24798 0,00028 0
0,15589 0,00241 2,27793 0,00029 0
0,16365 0,00251 2,30572 0,00029 0
0,17173 0,00262 2,33116 0,0003 0
0,18014 0,00273 2,35409 0,0003 0
0,1889 0,00284 2,37432 0,00031 0
0,19802 0,00296 2,39165 0,00031 0
0,20751 0,00309 2,40589 0,00032 0
0,21739 0,00322 2,41682 0,00033 0
0,22767 0,00335 2,42424 0,00033 0
0,23837 0,00349 2,4279 0,00034 0
0,2495 0,00364 2,42758 0,00035 0
0,26107 0,00379 2,42303 0,00035 0
0,27312 0,00395 2,41398 0,00036 0
0,28565 0,00411 2,40016 0,00037 0
0,29868 0,00428 2,3813 0,00038 0
0,31223 0,00445 2,3571 0,00038 0
0,32632 0,00463 2,32725 0,00039 0
0,34097 0,00482 2,29143 0,0004 0
0,35621 0,00502 2,24932 0,00041 0
0,37205 0,00522 2,20056 0,00042 0
0,38852 0,00543 2,14479 0,00043 0
0,40564 0,00565 2,08164 0,00044 0
0,42343 0,00588 2,01071 0,00045 0
0,44193 0,00612 1,9316 0,00046 0
0,46115 0,00636 1,84387 0,00047 0
0,48113 0,00662 1,74709 0,00048 0
0,50189 0,00688 1,64079 0,00049 0
0,52346 0,00715 1,52449 0,0005 0
0,54587 0,00744 1,39769 0,00052 0
0,56916 0,00773 1,25987 0,00053 0
0,59335 0,00803 1,11048 0,00054 0
5. CONCLUSÃO

Com os resultados obtidos, é possível se notar que, para um número grande de

subintervalos todos os métodos, menos o método de Taylor de Ordem 4, se
aproximam da resposta exata da Equação Diferencial requerida. A possível explicação
para a divergência do método de Taylor de Ordem 4 se dá pelo possível erro de
derivação da EDO, já que o método exige que seja obtida a quarta derivada da Edo,
acarretando em possíveis erros humanos já que a obtenção dessa derivada pode ser
de alta complexidade. O comportamento do método demonstra insatisfatório a partir
de x=0, mas é perceptível a tendência de retornar a um erro mínimo para x>2,5, o que
pode significar que em um intervalo maior o erro volte a aumentar novamente,
mostrando que o método é ineficiente e/ou que há algum erro, como anteriormente
citado.
Ainda sobre um número grande de subintervalos, é perceptível a tendência a
gerar um erro maior conforme nos afastamos no ponto inicial no método de Euler,
visto que este método é o mais simples entre todos os métodos e devido a isso gera
erros maiores que os outros.
Quanto a um pequeno número de subintervalos o método de Euler demonstra-
se insatisfatório visto que conforme o valor de x aumenta o erro do método também
aumenta e de forma exponencial. O mesmo é possível notar para o método de Taylor
de Ordem 2, entretanto este possui uma precisão maior que o método de Euler, visto
que Taylor de Ordem 2 utiliza-se da segunda derivada da EDO.
Os métodos de Runge-Kutta, tanto de Ordem 2 quanto Ordem 4, demonstram-
se extremamente precisos. Obviamente, o método de Runge-Kutta de Ordem 4 é mais
preciso devido a sua melhor aproximação, visto que, desde os menores números de
subintervalos, os valores se aproximam e muito da resposta exata, basta notar que os
pontos que denotam o erro do método acompanham o eixo x (y≈0) desde n=4 e
chegando erro zero a partir de n=50. O método de Runge-Kutta de Ordem 2 segue a
mesma tendência de erro próximo de zero a partir de n=10, mas para números
maiores de subintervalos ambos os métodos demonstram-se muito precisos.