Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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) = 2*x(i)*y(j) + 3*sin(pi*x(i))*(sinh(pi*y(j))/sinh(pi));
- 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');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement