% unit square
%
load coordinates.dat; coordinates(:,1)=[];
load connectivity.dat; connectivity(:,1)=[];
% load dirichtlet.dat; dirichlet(:,1)=[];
% at the boundary nodes of the unit square the temperature is prescribed with zero
%
solution(:,1)=zeros(size(coordinates,1),1);
%
trisurf(connectivity,coordinates(:,1),coordinates(:,2), solution(:))
view(-60,45);
title('discretization of the unit square')
% 





