Discussion :

[furtivesquirrel]

Bonjour,
Lors de la construction d'un intégrateur pour simuler l'évolution temporelle d'un glacier monodimensionnel, nous devons implémenter une certaine équation contenant un terme de source. Cependant, nous rencontrons un problème pour ce qui est d'implémenter ce terme de source. Voici la définition de celui-ci :


Ainsi que son implémentation dans notre code :



Dans notre code, seule l'implémentation de S0 est réalisée à cet endroit.
Pour implémenter la partie concernant S(x,t,h), nous avons utiliser l'option nonnegative que nous avons ajoutée à notre intégrateur ode45.



Apparemment, nous avons un problème avec ce terme de source mais nous ne voyons pas d'autre façon d'implémenter celui-ci. Quelqu'un pourrait-il nous aider ?
Si besoin, je peux envoyer l'entièreté de notre énoncé et de notre code par mail.

Merci beaucoup,
Terry