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).