% Plot reconstructions
%
% Samuli Siltanen May 2012

% Plot parameters
lwidth = 2;
fsize  = 16;

% Create plot objects
load data/recon xvec recon
tBIE_recon = recon;


% Plot non-noisy tBIE, Texp and t0
figure(2)
clf
plot(linspace(0,1,200),sigma(linspace(0,1,200)),'k','linewidth',lwidth)
hold on
plot(xvec,recon,'m','linewidth',lwidth)
axis([0 1 .5 3])
set(gca,'ytick',[0:.5:3],'fontsize',fsize)
set(gca,'xtick',[0:.2:1],'fontsize',fsize)
box off
legend('conductivity','reconstruction')

% Save image to file
print -dpng recon.png

