Questão 9

close all, clear all, clc;
%malha computacional bidimensional
N=1000; a=0; b=1;
h=(b-a)/(N-1);
x=a:h:b;
y=x;
for i=1:N
    for j=1:N
        V(i,j) = 10*exp(-pi*y(j)/b)*sin(pi*x(i)/b);
    end
end
[X,Y] = meshgrid(x,y);
mesh(X,Y,V);

colorbar
label(1)=xlabel('x');
label(2)=ylabel('y');
label(3)=zlabel('V(x,y)');
label(4)=title('Soluão exata');