Résolution itérative d'un Laplacien 2D
Bonjour tout le monde.
j'ai un problème de syntaxe dans mon algo et je m'arrache les cheveux dessus !
Méthode itérative de résolution d'un laplacien 2D, définition de la fonction d'itération :
Code:
1
2
3
4
5
6
7
8
9
10
| def iteration (u,x,y):
n=len(x)
dx2=[(x[-1]-x[0])/(n-1)]**2
dy2=[(y[-1]-y[0])/(n-1)]**2
l=u.shape[0],m=u.shape[1]
u1=zeros((n,n))
for i in range(l):
for j in range(m):
u1[i,j]=(dx2*dy2+u[i+1,j]+u[i-1,j]*dy2+(u[i,j+1]+u[i,j-1]*dx2)/(dx2+dy2)
return u1 |
Le terminal me dit que le problème de syntaxe vient de la ligne 10 return u1
Si vous voyez quelque chose, je suis preneur, moi je sèche ...
Merci bien,
Cordialement.
