Olá Pessoal,
A Dani já compartilhou, segue anexo o meu. Espero o de vocês.
Valeu,
Olá a todos,
O código que compartilho possui convergência pela norma do infinito, já a norma_2 não está convergindo como o esperado.
Seguem as seguintes especificações do código (que está em anexo):
u_exata: val = b*(cos(x+y))
forçante: val = (0.5)*b*(cos(x+y))
b = e^(-t)
comprimento = 1.0
largura = 1.0
mhi = 0.75
tempo_final = 0.5
tolerância: dif_relativa > (dt/10.0)
Pessoal, desculpem a demora. Segue meu código.
u_exata: mi * exp(-0.1 * t) * cos(x + y), mi = 1;
a = 0; b = 4; c = 0; d = 4;
Nx -> começa com 8; Ny -> começa com 8 (dobra a cada caso);
tMax = 0.5 (tempo final);
tolerância: diferença em relação à equação exata = dt / 10).