Bonjour à tous. Je ne sais pas si je suis sur le bon forum mais je me
lance. (Désolé si je suis dans l'erreur).
Je travail présentement sur un programme utilisant la méthode
Runge-Kutta-Fehlberg (RKF45). Je trouve plein de livre ou site web me
donnant des infos contradictoires.
Voici mes questions :
Je calcule mon cofficient d'ajustement pour mon pas h de la façon suivante :
E = abs(zi1 - yi1);
if(E == 0.0)
B = 1.0;
else
B = pow(e/(2*E),0.2);
1) est-ce bien ^0.2 ou ^0.25 ?
2) je vois souvent 2*E cependant mes notes de cours me donne pow(e/(E),0.2);
Par la suite si mon cofficient B est trop grand (par exemple si j'ajuste
mon pas h avec 170) il se peut que cela soit trop. Donc plusieurs font une
vérification du genre :
si (B > seuil) alors
B = seuil
ainsi je préviens d'avoir un pas h trop grand pour le prochain calcul.
Dois-je faire cette vérification ?