Bonjour,
J'ai de la misère a finir le calcul de D par la méthode de la bissection parce que je ne suis pas capable de trouver la valeur de C(1) pour un rayon donné. C'est a dire, la boucle de la concentration en fonction du temps et ensuite la variation par rapport au rayon avec 5001 nœuds n'aboutissent pas a un résultat parce que ca prend beaucoup trop de temps et ne finit jamais, bien que ce n'est pas une boucle infinie. Si je fais pour n=100 noeuds la boucle marche, mais je trouve une valeur de 0 pour la concentration a t=5hrs et a r=1cm. Auriez vous une idée de la source de mon erreur?
Merci.