Solução clássica da EDP de Laplace. Código em Python--------------------------------------import numpy as npimport matplotlib.pyplot as pltmaxIter 500lenX lenY 21 delta 1Ttop 50Tbottom 0Tleft 10Tright 10Tguess 30colorinterpolation 50colourMap plt.cm.jet X, Y np.meshgrid(np.arange(0, lenX), np.arange(0, lenY))T np.empty((lenX, lenY))T.fill(Tguess)T[(lenY-1), ] TtopT[1, ] TbottomT[, (lenX-1)] TrightT[, 1] Tleftprint("um momento")for iteration in range(0, maxIter) for i in range(1, lenX-1, delta) for j in range(1, lenY-1, delta) T[i, j] 0.25 (T[i1][j] T[i-1][j] T[i][j1] T[i][j-1])print("Iteração Excedida")plt.title("Variacao da Temperatura")plt.contourf(X, Y, T, colorinterpolation, cmapcolourMap)plt.colorbar()plt.show()print("")---------------------------------------------------------------------
Compartilhar