import matplotlib.pyplot as plt
from matplotlib import cm
import numpy as np

#Fonctions

#----------------A compléter

#Initialisation des paramètres
Nx = 100
Nz = 50
Vnuage = 100E6
Vsol = 0
Hnuage = 500
L = Hnuage/(Nz-1)*(Nx-1)
dx = L/(Nx-1)
dz = Hnuage/(Nz-1)

#Initialisation du tableau V[j][i]
# i repère les colonnes (abscisse) et j les lignes (ordonnée)

#-------------------A compléter
    
    
#Paramères des itérations
epsilon=Vnuage/1E4
ecartmax=10*epsilon
Nmax=10000
Niteration=0

#Remplissage du tableau

#-------------------A compléter


#Affichage du tableau
#print('Nombre itérations : ',Niteration)
#plt.subplot(211)
plt.imshow(V,cmap=cm.gist_gray)
plt.colorbar()
plt.axis('off')
plt.title('Potentiel életrostatique (V)')
plt.show()



#Calcul du champ électrique

#------------------A compléter

"""
plt.subplot(212)
plt.imshow(E,cmap=cm.gist_rainbow)
plt.colorbar()
plt.axis('off')
plt.title('Champ életrostatique (V/m)')
plt.show()
"""