etude numerique de diffusion
Bonjour, je suis bloquer sur un probleme de physique
Écrire un script Python pour tracer la condition initiale c(x,0) comme un vecteurde 201 éléments dans le domaine x ∈[−1,1]. avec
cr(xr,tr)= c0 /2 erfcxr 2√Dtr
Tracer la solution analytique au temps t = 0.01 sur le même graphique. On pourra utiliser la fonction Python erfc que l’on pourra charger grâce à la commande :
1 from scipy.special import erfc
Je n'ai fais que le debut, ouis je bloque
Code:
1
2
3
4
5
6
| import os, glob, matplotlib, time
import numpy as np
from scipy.special import erfc
from matplotlib import use
use('agg')
import matplotlib.pyplot as plt |
J'avais oublier la deuxieme ligne
merci d avance
bien à vous
mh
