Discussion :

[Jdevv]

Bonsoir,

Cela fait déjà 1 semaine que j'essaye de résoudre l'équation :

-(laplacien)u + V(x,y) * u = 1

Qui n'est d'autre qu'une equation elliptique (trouvable dans PDE toolbox) avec c=1 , a=V(x,y), f=1 (notation de Matlab). V(x,y) est une fonction que j'ai crée censée représentée un potentiel.

Le problème c'est que Matlab n'arrive pas à comprendre a=V(x,y) (même dans le cas trivial V(x,y)=1)

Quelqu'un sait comment lui faire comprendre que le coefficient a est une fonction de l'espace ?

Merci d'avance.

[Gooby]

Bonjour,

Peux tu nous montrer le code qui définit ta fonction V(x,y)?

[Jdevv]

Code :

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

1
2
3
4
5
6
7
function v = V(x,y)
     global M
     [i,j]=size(M);
     intx=int8(floor(i.*x))+1;
     inty=int8(floor(j.*x))+1
 
    v = M(intx,inty);

M étant une matrice deja crée par

Code :

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

global M;  M= rand(100,100);

L'erreur est :

Code :

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

 Expression evaluates to wrong size. Must be scalar or row vector. In a system case, pass first or second row; for example u(2,:)

[Gooby]

Comment fais tu appel à cette fonction? Peux tu nous montrer la ligne associée à cette erreur?