| Autor: |
Moawwad El-Mikkawy, M. M. M Eisa |
| Rok vydání: |
2003 |
| Předmět: |
|
| Zdroj: |
Applied Mathematics and Computation. 143:259-267 |
| ISSN: |
0096-3003 |
| DOI: |
10.1016/s0096-3003(02)00358-2 |
| Popis: |
The well-known Dormand-Prince embedded RK 5(4) 7FM algorithm [J. Comput. Appl. Math. 6 (1980) 19] is of the FSAL type and uses seven stages per step. This algorithm has been recommended by Shampine [Math. Comput. 46 (1986) 135] as a candidate for an efficient production RK code. In fact the new MATLAB function ode45 is based on this algorithm. Later on another efficient RK 5(4) 7 embedded algorithm, also having orders 5 and 4, is developed by Sharp and Smart [SIAM J. Sci. Comput. 14 (1993) 338]. The last algorithm uses seven stages per step and it is of the non-FSAL type. The current paper shows that by using seven stages per step a general four-parameter, non-FSAL embedded RK algorithm having orders 6 and 4 may be designed. A special algorithm, called RK 6(4) 7 new is obtained by using suitable choices for the free parameters. This new algorithm together with the RK 5(4) 7FM in [J. Comput. Appl. Math. 6 (1980) 19] and the RK 5(4) 7 in [SIAM J. Sci. Comput. 14 (1993) 338] are applied to some test problems, which have known exact solutions. It is found that the new algorithm is competitive comparing with the algorithms in [J. Comput. Appl. Math. 6 (1980) 19] and in [SIAM J. Sci. Comput. 14 (1993) 338]. |
| Databáze: |
OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect