Problème de code: Erreur 9999

**johannj** · 17/12/2009, 12h08

Bonjour le forum,

J'ai un problème avec le code suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
x=linspace(100,1500,3000);

cs1=input("boron concentration:");

xbasc();
Ea=-3.66;q=1.60217733e-19;k=1.380658e-23;
kbis=k/q;

gap=1.166-(4.73e-4*((273.15+x).^2)./(636+x+273.15));

Nc=2.86e19.*((273.15+x)./300).^(1.58);
Nv=3.1e19.*((273.15+x)./300).^(1.85);
A=sqrt(Nc.*Nv);
ni=A.*exp(-gap./((2*kbis).*(x+273.15)));

ne=4.65e21*exp(-0.39./(kbis*(x+273.15)));

D0=3.85*exp((Ea/kbis)./(x+273.15));
stacksize("max");
Di=4.44*exp(-4./(kbis.*(273.15+x)));

Dtail1=Di.*(cs1^3)./(ni.*ne.^2);


plot2d1(x, D0,logflag="nl",frameflag=1,rect=[700,1e-18,1200,1e-10])

plot2d1([x],[Dtail1],style=[2],logflag="nl", frameflag=1,rect=[700,1e-18,1200,1e-10])

La ligne en rouge correspond à l'erreur: celle-ci semble venir du signe "/". Quoique j'écrive, lors de l'exécution de ce code, une erreur 9999 (inconsistent element-wise operation) apparait pour cette ligne. J'ai beau me battre, là je ne comprends pas!!

Quelqu'un a-t-il une idée SVP?

Merci de votre aide,
Johann

**Jerome Briot** · 17/12/2009, 14h20

Quelles sont les dimensions des variables Di, cs1, ni et ne ?

**johannj** · 17/12/2009, 15h52

Bonjour Dut,

Et bien je dois dire que je ne sais pas

Ce qui est sur, c'est que cs1 est un scalaire.
Di, ne et ni dépendent tous de x, donc j'imagine qu'ils ont la même dimension que x, donc un vecteur colonne de "3000 points".

Merci de ton aide,
Johann