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 : Sélectionner tout - Visualiser dans une fenêtre à part
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