function []=Plotinvsol(rho,g,H);

%Plotinvsol Plots the solution of 2D EIT problem with piecewise constant basis
% Function []=Plotinvsol(rho,g,H);
% plots the solution of 2D EIT problem.
%
% INPUT
%
% rho = the solution, constant on each element
% g = node coordinates (x,y) (N x 2 matrix)
% H = connectivity matrix 


[nH,mH]=size(H);

axis('off'),axis('square')
set(gcf,'defaultpatchedgecolor','none');

for ii=1:nH
  Hii=g(H(ii,:),:);
   patch(Hii(:,1),Hii(:,2),rho(ii));		
end